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ABSTRACT 


Thia  report  describes  the  results  of  a  study  of  the  dynamic 
stability  and  response  during  the  feathering  and  folding 
of  the  blades  of  a  folding-proprotor  VTOL.  This  study  involved 
the  development  of  analytical  methods  of  predicting  the 
dynamic  characteristics  during  feathering  and  folding#  corre¬ 
lation  of  the  theoretical  methods  with  experimental  data,  a 
dynamic  analysis  of  a  representative  folding-proprotor  VTOL 
design,  and  a  parametric  study  to  identify  design  factors 
which  can  be  used  to  control  dynamic  characteristics  during 
feathering  and  folding.  Correlation  of  the  theory  with 
measured  dynamic  stability  and  response  characteristics 
during  feathering  and  folding  is  good,  and  indicates  the 
analytical  methods  can  be  used  with  confidence.  The  results 
of  the  dynamic  analysis  indicate  satisfactory  stability  and 
response  characteristics  can  be  achieved  in  a  66,000-pound 
gross-weight  class,  folding-proprotor  VTOL.  The  results  of 
the  parametric  study  are  summarised  in  terms  of  design  guide¬ 
lines.  This  volume  contains  the  development  of  the  theory, 
correlation  of  theory  with  experimental  data,  and  the  dynamic 
analysis  and  parametric  study.  Volume  II  is  a  guide  to  the 
computer  programs  and  contains  FORTRAN  listings. 
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SECTION  I 


INTRODUCTION 


The  folding- proprotor  VTOL  concept  combines  the  low-speed 
flight  characteristics  of  the  helicopter  with  the  high  subsonic 
cruise  characteristics  of  the  conventional  turbojet  airplane. 
Because  of  this  highly  desirable  combination  of  flight 
characteristics,  the  folding-proprotor  concept  is  currently 
being  investigated  for  both  military  and  commercial  applica¬ 
tions  . 

The  folding-proprotor  VTOL  (hereafter  referred  to  as  FPR) 
principle  of  operation  is  illustrated  in  Figure  1.  Takeoff  is 
accomplished  in  helicopter  mode  with  the  nacelles  vertical. 
After  takeoff,  the  FPR  is  accelerated  to  an  airspeed  at  which 
wing  lift  can  support  the  aircraft.  The  nacelles  are  then 
mechanically  tilted  90-degrees  so  that  the  proprotors  act  as 
propellers.  During  vertical  and  low-speed  flight,  control  of 
the  aircraft  is  by  use  of  the  usual  helicopter  collective  and 
cyclic  pitch  controls.  These  controls  are  washed-out  as  the 
nacelles  are  tilted  forward,  and  at  90  degrees  nacelle  tilt, 
control  is  by  conventional  fixed-wing  control.  The  nacelle- 
tilting  process  is  operationally  flexible;  it  can  be 
accomplished  over  a  wide  range  of  airspeeds  and  maneuver  con¬ 
ditions,  and  with  or  without  power. 

After  the  nacelles  are  tilted,  flight  can  be  continued  with 
the  proprotors  acting  as  propellers  or  the  proprotors  can  be 
feathered  and  folded,  and  fan  engine  thrust  used  for  propul¬ 
sion.  Proprotor  feathering  and  folding  is  accomplished  by 
transferring  power  from  the  proprotors  to  high  bypass  ratio 
fans  (compound  engines  capable  of  delivering  shaft  power  or 
fan  thrust  are  employed)  until  the  proprotors  are  windmilling. 
The  proprotors  are  then  declutched  from  the  engines  and  the 
blades  feathered.  When  stopped,  the  proprotors  are  indexed 
and  the  blades  folded.  With  the  blades  folded,  the  FPR  is 
operated  as  a  conventional  jet  aircraft.  For  landing,  the 
tilt/feather/fold  sequence  described  above  is  reversed. 

Preliminary  investigations  were  made  in  1967-68  to  determine 
feasibility  of  the  FPR  concept.  ^->2  These  concentrated  on  the 
feather/fold  process  since  helicopter,  nacelle-tilting,  and 
propeller  modes  of  flight  were  found  to  be  feasible  during 
the  XV-3  flight  tests,-'  and  since  technical  considerations 
after  the  btades  are  folded  are  expected  to  be  similar  to  those 
of  conventional  fixed  wing  aircraft. 

The  preliminary  investigations  included  small-scale  model  tests 
as  well  as  analytical  studies  and  showed  the  feather/fold  pro¬ 
cedure  to  be  relatively  simple  and  straightforward.  However, 
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FAN  PROPULSION  - 
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ROTORS  FOLDED 
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FIGURE  1.  Folding-Proprotor  Principle  of  Operation 
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they  also  identified  several  areas  requiring  further  investiga¬ 
tion.  Among  these  were  the  dynamic  and  aeroelastic  character¬ 
istics  during  proprotor  feathering  and  folding. 

A.  FOLDING  PROPROTOR  DYNAMIC  CONSIDERATIONS 

The  analytical  and  experimental  investigations  of  the  dynamics 
of  the  tilting-proprotor  VTOL4#5'*  and  the  results  of  the  pre¬ 
liminary  investigations  of  folding-proprotor  dynamics,  have 
provided  a  good  background  for  an  understanding  of  FPR  dynamic 
considerations.  Some  of  these  considerations  are  discussed 
below. 

1.  Dynamic  Stability 

Proprotor  instability,  similar  in  nature  to  propeller  whirl 
flutter,  was  encountered  during  tests  of  the  XV-3  converti- 
plane.  Subsequent  investigations  have  led  to  an  understanding 
of  proprotor  instability  and  have  provided  design  approaches 
which  can  provide  an  adequate  margin  of  stability  for  tilting- 
proprotor  VTOLs. 

The  airloads  generated  by  precession  are  the  cause  of  proprotor 
instability.  However,  the  proprotor' s  flapping  freedom  causes 
a  basic  difference  in  the  manner  in  which  the  precession  gen¬ 
erated  aerodynamic  forces  act  or.  the  pylon. 

In  the  case  of  the  propeller  the  rigid  attachment  to  the  hub 
processes  the  propeller  in  response  to  nacelle  pitching  (or 
yawing) .  If  the  nacelle  pitches  up  the  angle  of  attack  gen¬ 
erated  airloads  cause  a  yawing  moment  at  the  hub.  This  yawing 
moment  acts  as  negative  damping  on  the  nacelle  backward  whirl 
mode  and  at  some  airspeed  will  exceed  the  combined  nacelle 
structural  damping  and  aerodynamic  damping  caused  by  trans¬ 
verse  velocities  over  the  propeller  disc.  Flutter  in  the 
backward  whirl  mode  results. 

In  the  case  of  the  proprotor  airloads  are  required  to  precess 
the  proprotor  because  of  the  blade  flapping  freedom.  (In 
a  vacuum  the  proprotor  would  not  precess  in  response  to  pylon 
pitching  or  yawing.)  In  order  to  generate  the  precessional 
airloads  a  change  in  blade  angle  of  attack  is  required.  This 
means  the  rotor  must  develop  a  flapping  angle  with  respect  to 
its  control  plane  (swashplate) .  Consequently,  a  proprotor 
must  lag  the  pylon  motion.  This  leads  to  frequency  dependent 
aerodynamic  flapping  moments  and  hub  inplane  shear  forces. 

In  contrast,  a  propeller's  aerodynamics  forces  and  moments 
are  essentially  independent  of  the  frequency  of  the  nacelle 
motions . 

From  the  standpoint  of  proprotor  stability,  the  normal  shear 
force  is  the  most  significant.  It  can  be  conveniently  expressed 
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as  H~  ,  the  component  of  the  normal  force  in-phase  with  pylon 
pitch  and  Hq,  the  component  in-phase  with  pylon  pitch  rate. 

Ha  acts  as  3  negative  spring  reducing  the  pylon  natural  fre¬ 
quency.  Depending  on  the  pylon  pitching  frequency  H„  can  act 
as  either  positive  or  nagative  damping  and  thus  be  stabilizing 
or  destabilizing.  Small  side  forces  are  also  generated  but 
are  not  of  consequence  in  proprotor  instability. 

In  the  case  of  the  FPR  aircraft,  the  proprotor  stability 
requirements  are  somewhat  relaxed  since  the  blades  are  folded 
for  operation  at  high  speed.  Based  on  tilt-proprotor  design 
experience,  the  wing  stiffnesses  resulting  from  strength 
requirements  should  provide  an  adequate  proprotor  stability 
margin.  However,  the  provision  of  feather/fold  capability  may 
require  design  features  which  have  a  destabilizing  influence 
on  proprotor  stability.  For  example,  a  longer  mast  may  be 
required  to  accommodate  blade  folding.  Another  consideration 
is  that  during  feathering  and  unfeathering  proprotor  natural 
frequencies  can  become  coincident  with  wing/pylon  natural 
frequencies.  Any  resulting  adverse  dynamic  or  aerodynamic 
coupling  could  then  cause  instability.  An  example  of  this 
was  encountered  during  the  dynamic  model  test  conducted  under 
this  program  when  coincidence  of  the  blade  flapping  and  wing- 
beam  bending  mode  was  found  to  cause  the  blade  flapping  mode 
to  become  unstable  thereby  requiring  a  flapping  lockout  to 
insure  stability. 

Wing  and  blade  flutter  and  divergence  become  considerations 
when  the  proprotor  is  stopped.  Wing  divergence  is  of  concern 
since  the  proprotor  is  located  well  ahead  of  the  wing  torsional 
axis.  Wing  flutter  does  not  appear  to  be  a  serious  problem 
because  of  the  forward  location  of  the  pylon  center  of  gravity. 
However,  coupled  blade/wing  flutter  could  be  a  problem.  As  the 
blades  are  folded,  the  blade  natural  frequencies  change  and  can 
become  coincident  with  wing/pylon  frequencies.  Also,  at  large 
fold  angles,  when  the  blades  are  in  close  proximity  to  the 
nacelle,  aerodynamic  buffetting  could  occur.  Blade  divergence 
does  not  appear  to  be  a  problem  since  the  blades  are  not  sub¬ 
ject  to  large  forward  sweep  angles  or  reversed  flow,  in  con¬ 
trast  to  the  case  of  the  edgewise  stoppable-rotor  concept. 

2 .  Blade  Flapping  and  Loads 

The  XV-3  flight  tests  showed  that  provision  of  adequate  flap¬ 
ping  clearance  between  the  proprotor  blades  and  the  wing  lead¬ 
ing  edge  in  the  airplane  mode  is  required  for  tilting-proprotor 
aircraft.  In  general,  the  required  clearance  is  proportional 
to  the  square  of  the  proprotor  advance  ratio,  V/ftR,  and  can 
be  reduced  if  the  blade  flapping  natural  frequency  is  detuned 
from  one-per-rev  by  using  flapping  restraint  and/or  pitch-flap 
coupling. 
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For  the  FPR,  the  large  advance  ratios  occurring  during  blade 
feathering  will  cause  excessive  flapping  unless  flapping  is 
restrained.  An  important  result  of  this  study  is  the  finding 
that  rapid  rates  of  feathering  of  the  blades  increases  the 
flapping  caused  by  shaft  angle  of  attack. 

However,  restraining  flapping  or  the  use  of  very  stiff  hinge¬ 
less  blades,  increases  one-per-rev  loads  in  the  rotor,  since 
flapping  acts  to  alleviate  one-per-rev  airloads.  While  only  a 
limited  number  of  cycles  is  incurred  during  a  stop/start  cycle, 
these  one-per-rev  loads  may  limit  the  maneuver  capability  in 
helicopter  and  conversion  modes. 

Another  blade  load  consideration  is  maneuvers  and  gust  en¬ 
counters  when  the  proprotor  is  at  low  rpm  or  is  stopped  but  not 
folded.  Blade  stresses  under  these  conditions  could  be  greater 
than  those  at  normal  operating  rpm  at  the  same  forward  speed. 

Blade  load  build-up  due  to  traversing  blade  natural  frequency 
resonances  during  feathering  were  not  a  problem  during  pre¬ 
liminary  model  tests.  However,  the  model  proprotors  had  the 
lowest  inplane  natural  frequency  of  the  blades  designed  above 
operating  rpm.  In  rotor  types  where  the  lowest  inplane 
natural  frequency  is  below  operating  rpm,  high  blade  loads 
might  occur  when  the  blades  are  feathered  and  the  mode  transits 
one-per-rev  resonance. 

3.  Dynamic  Response 

Initially  it  was  thought  that  the  dynamic  response  during 
proprotor  feathering  would  be  negligible.  This  was  based  on 
the  assumption  that  the  axial  flow  condition  of  the  proprotor 
would  not  induce  significant  oscillatory  excitation. 

However,  during  wind  tunnel  tests  of  small-scale  dynamic 
models,  relatively  high  vibrations  of  the  wing  and  pylon  were 
recorded.  These  occurred  at  the  blade  passage  frequency  with 
peak  vibrations  occurring  when  the  proprotor  angular  velocity 
was  one-third  of  the  wing/pylon  natural  frequency  (a  three- 
bladed  proprotor  was  employed  on  the  model) .  The  source  of 
this  vibration  was  determined  to  be  wing-blade  aerodynamic 
interference,  principally  the  change  in  the  blade  airloads  as 
it  passed  immediately  in  front  of  the  wing. 

The  principal  risk,  insofar  as  this  blade  passage  frequency 
vibration  is  concerned,  appears  to  be  the  acceptability  of 
the  level  of  crew  station  and  cabin  vibrations  during  feather¬ 
ing  or  unfeathering.  Of  course,  the  design  of  the  nacelle 
must  take  into  account  these  vibrations. 

In  the  preliminary  tests,  only  very  slow  rates  of  feathering 
were  possible  because  of  mechanical  limitations.  A  fast  rate 
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of  feathering  appears  to  be  desirable  since  it  will  reduce 
resonance  buildup.  This  was  investigated  in  this  study  and  is 
discussed  in  Sections  III  and  V. 

B.  PROGRAM  APPROACH 

As  stated  earlier,  the  objective  of  this  program  of  study  was 
to  develop  analytical  methods  of  predicting  dynamic  behavior 
during  proprotor  feathering  and  folding.  To  accomplish  this, 
the  study  was  conducted  in  three  phases. 

In  Phase  I,  three  digital  computer  programs  were  developed;  one 
treating  the  feathering  and  stopping  of  the  proprotor,  and  the 
other  two  blade  folding.  The  programs  are  based  on  existing 
computer  programs,  which  were  extensively  modified  to  simulate 
the  FPR  configuration. 

Phase  II  involved  the  testing  of  a  small-scale,  aeroelastically- 
scaled  model.  Data  were  obtained  on  dynamic  stability  at  low 
rpm  and  on  the  dynamic  response  during  feathering.  The  effect 
of  feathering  rate  was  investigated.  The  accuracy  of  the  com¬ 
puter  programs  developed  in  Phase  I  was  then  evaluated  by  com¬ 
paring  predicted  dynamic  characteristics  with  those  recorded 
during  the  model  test  and  with  the  results  of  other  related 
tests. 

The  dynamic  characteristics  of  a  66,000-pound  gross  weight 
FPR  design  was  investigated  in  Phase  III.  This  investigation 
provided  the  basis  for  a  parametric  study  to  identify  design 
parameters  having  significant  influence  on  dynamic  characteris¬ 
tics  during  blade  feathering  and  folding. 

,  . 

This  Volume  contains  a  description  of  the  analytical  methods 
employed,  the  results  of  the  model  test  and  program  correla¬ 
tion,  and  the  parametric  study.  Volume  II  contains  the  infor¬ 
mation  required  for  utilization  of  the  computer  programs  and 
FORTRAN  listings. 
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SECTION  II 


DEVELOPMENT  OF  THEORY 


In  order  to  simplify  the  development  of  a  dynamic  theory  the 
transition  mode  was  considered  to  be  separated  into  two  distinct 
phases:  (1)  a  blade  feathering  or  unfeathering  phase,  in 

which  the  proprotor  is  stopped  or  started  and,  (2)  a  blade 
folding  or  unfolding  phase.  These  overlap  only  at  the  con¬ 
dition  of  zero  rpm  and  zero  fold  angle.  By  developing 
separate  analyses  for  each  phase,  the  most  applicable  dynamic 
and  aerodynamic  theory  was  used. 

The  theory  treating  the  blade  feathering/unfeathering  phase 
is  based  on  a  Bell  Helicopter  Company  tilt-proprotor  aero- 
elastic  analysis,  called  DYN5.  The  folding -proprotor  version 
of  that  analysis  is  referred  to  in  this  report  as  the  Feather- 
Stop  Analysis  and  its  digital  computer  mechanization  is  named 
program  ARAP06.  The  theory  treating  the  blade  folding  and 
unfolding  phase  is  based  on  a  fixed  wing  flutter  analysis 
developed  for  the  U.S.  Air  Force  by  the  Aerospace  Corporation. 8 
It  is  referred  to  in  this  report  as  the  Folding  Dynamics 
Analysis  and  its  digital  computer  mechanizations  are  named 
DFAL17  and  DFAL18  (two  separate  programs). 

This  section  is  intended  to  provide  the  reader  with  an  under¬ 
standing  of  the  theory  and  assumptions  used  in  the  analyses. 

No  attempt  is  made  to  give  the  complete  equations  of  motion 
as  they  are  rather  extensive.  However,  the  equations  which 
define  the  math  model  are  given  to  provide  the  reader  with 
insight  into  the  details  of  the  analysis.  The  FORTRAN  listings 
for  each  of  the  computer  programs  are  incluaed  in  Volume  II 
and  serve  as  the  documentation  of  the  equations  of  motion. 

A  brief  discussion  of  the  computer  programs  is  also  given  in 
this  section. 

A.  FEATHER-STOP  ANALYSIS 

The  requirement  to  predict  proprotor  stability  and  whirl 
flutter  characteristics,  and  the  blade  and  pylon  response 
during  feathering  or  unfeathering  dictate  a  fully  coupled 
wing-pylon-proprotor  math  model  in  which  certain  dynamic  and 
aerodynamic  characteristics  are  rigorously  accounted  for. 

The  feather/unfeather  simulation  poses  unusual  requirements 
for  rotor  aeroelastic  analysis;  (1)  that  the  rotor  angular 
velocity  be  an  independent  variable,  and  (2)  that  the  blades 
be  treated  as  individual  blades  rather  than  assuming  identical 
response  for  each  blade.  The  second  requirement  is  the  result 
of  wing-rotor  aerodynamic  interference  considerations. 
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1.  Description  of  the  Math  Model 


An  overview  of  the  math  model  for  the  feather-stop  analysis  is 
shown  in  Figure  2.  The  degrees  of  freedom  include:  (1)  rigid- 
body  blade  flapping,  feathering,  and  rotation,  (2)  blade  beam 
and  chord  bending  and  torsion,  (3)  pylon  pitching  and  yawing 
with  respect  to  tne  wing-tip,  and  (4)  wing  beam  and  chord 
bending  and  torsion. 

The  development  of  the  equations  of  motion  are  greatly  simpli¬ 
fied  by  treating  as  separate  systems  the  wing,  pylon,  and  hub, 
the  proprotor  rigid-body  flapping,  feathering  and  rotation,  and 
the  blade  elastic  motions  as  illustrated  in  Figure  3.  The  equa¬ 
tions  of  motion  for  the  proprotor  rigid-body  motion  are  obtained 
using  the  moment  of  momentum  theorem  and  those  of  the  wing-pylon 
hub  and  blade  elastic  deformation  using  the  normal  mode  method. 
The  equations  are  joined  through  generalized  force  coupling. 

Four  Cartesian  coordinate  systems  are  employed  in  the  feather- 
stop  analysis.  The  right  hand  rule  is  used  in  all  cases. 

These  are  illustrated  in  Figure  4.  The  wing  coordinate  system, 

^F^F^F 

serves  as  the  reference  inertial  system  and  has  the  following 
sign  convention: 

Xp  is  positive  forward 

Yp  is  positive  outboard 

Zp  is  positive  downward 

The  pylon  coordinate  system,  XYZ,  is  located  at  the  inter¬ 
section  of  the  shaft  axis  and  the  pylon  conversion  axis  and 
rotates  as  the  pylon  converts.  The  sign  convention  of  the 
XYZ  system  is  such  that  in  helicopter  mode  (shaft  axis 
vertical),  X  is  positive  aft,  Y  outboard  and  Z  upward,  i.e. 
the  Z  axis  lies  along  the  shaft.  In  airplane  mode  (shaft 
axis  horizontal)  X  is  positive  upward,  Y  outboard,  and  Z 
forward.  Rotating  the  XYZ  coordinate  system  with  the  shaft 
generates  the  proprotor  rotating  coordinate  system,  xyz.  In 
this  third  coordinate  system  the  x  axis  is  parallel  to  the 
feathering  axis  and  the  negative  y  axis  parallel  to  the 
flapping  axis  of  the  designated  master  blade.  The  z  axis  lies 
on  the  rotor  shaft.  Flapping  the  xyz  system  about  the  negative 
y  axis  yields  a  fourth  coordinate  system  x'y'z'  which  is  used 
primarily  for  blade  airload  calculations. 

2 .  Proprotor  Equations  of  Motion 
a.  Rigid-Body  Motion 

The  equations  of  motion  for  the  rigid-body  motions  are  derived 
in  a  general  form  by  first  writing  down  the  position  vector 
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WING/PYLON  FLEXIBILITY 
•IS  REPRESENTED  BY  FOUR 
NORMAL  MODES  (STEADY 
WING  AERODYNAMICS) 


PYLON  Y 
HINGE 


7 


PYLON  "C\ 
PITCH— ^  \ 
HINGE 

(CONVERSION 

HINGE) 


PYLON  MASS- 


GIMBALLED  HUB 
WITH  FLAPPING 
RESTRAINT  AND 
PITCH- FLAP 
COUPLING, 
COLLECTIVE  AND 
CYCLIC  PITCH 


BLADE  HAS  RIGID  BODY 
FLAPPING,  FEATHERING,  AND 
ROTATIONAL  DEGREES  OF 
FREEDOM.  BLADE  FLEXIBILITY 
IS  REPRESENTED  BY  FOUR 
NORMAL  MODES  (QUASI-STEADY 
PROPROTOR  AERODYNAMICS) 


Figure  2.  Overview  of  Math  Model  for  Feather-Stop 
Analysis . 
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(1) 


(2) 


Wing-Pylon 

Dynamics 


-  NORMAL  MODE  APPROACH 
USED  TO  REPRESENT  WING- 
PYLON  SYSTEM 


Proprotor 
Rigid  Body 
Flapping  and 
Rotational 
Freedom 


FORCES  AND  MOMENTS  FROM 
PROPROTOR  TREATED  AS 
GENERALIZED  FORCES 
(COUPLING  WITH  PROPROTOR) 

-COUPLED  EQUATIONS  FOR 
TPP  PITCH  AND  YAW  AND 
ROTATION  USED  TO 
REPRESENT  GIMBAL  FREEDOM 


WING-PYLON  SYSTEM 
> COUPLING  REPRE¬ 
SENTED  BY  HUB 
MOTION 


(3) 


NORMAL  MODES  USED  TO 
REPRESENT  BLADE 
ELASTIC  MODES 


COUPLING  WITH  WING- 
PYLON  SYSTEM  AND 
PROPROTOR  RIGID  BODY 
MOTIONS  THROUGH 
GENERALIZED  FORCE 
(AERODYNAMIC  AND 
INERTIAL) 


Figure  3.  Variable  Groupings  for  Feather-Stop  Analysis. 
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of  a  discrete  mass  point  on  a  designated  master  blade.  The  mass 
point  velocity  and  acceleration  are  determined  by  successive 
differentiation  of  the  position  vector.  The  equations  of 
motion  are  then  formulated  by  applying  the  moment  of  momentum 
theorem,  which  states  that  the  time  derivative  of  the  total 
angular  momentum  of  the  system  is  equal  to  the  moment  of  the 
external  forces  about  a  given  point.  The  equations  of  motion 
for  the  remaining  blades  are  obtained  from  those  for  the  master 
blade  by  coordinate  transformation. 

Figure  5  shows  the  relationship  between  a  discrete  mass  point 
on  a  designated  master  blade  and  the  fixed  coordinate  system, 
XpYpZp.  The  position  of  the  pylon  conversion  axis-shaft  inter¬ 
section  (point  CA)  is  denoted  by  X0YoZo,  which  are  linear  dis¬ 
placements  referred  to  the  XpYpZp  coordinate  system.  0p  and 
are  the  pylon  pitching  and  yawing  respectively;  point  H  is 
tne  hub  centerline  (zero  flapping  hinge  offset  is  assumed) . 

The  discrete  mass,  mj_,  is  located  relative  to  point  H  by  the 
following  rigid-body  coordinates:  r^,  the  radial  distance 
along  the  feathering  axis,  and  6,  the  rigid-body  flapping,  yj., 
the  chordwise  position  from  the  feathering  axis,  and  $£,  the 
blade  element  pitch.  The  elastic  displacements  further  define 
the  location  of  m^:  ni  and  5i  are  the  blade  out-of-plane  and 
inplane  deflections,  respectively,  and  0£f  is  the  pitching 
deflection.  6^  and  account  for  the  radial  motion  of  the 
mass  point  as  the  blade  bends  elastically.  Since  ni  and  Ci  are 
both  perpendicular  to  the  undeflected  blade,  fictious  spanwise 
motion  is  introduced,  i.e.,  terms  indicating  "stretching"  of  the 
blade  are  generated  when  the  blade  deflects.  In  nonrotating 
beam  analysis  these  fictious  terms  are  not  important.  But,  for 
a  rotating  beam,  fictious  Coriolis  forces  are  generated.  To 
correctly  represent  the  blade  bending  deflections  and  ej.  are 
used  to  cancel  the  fictious  spanwise  motion,  i.e.,  to  represent 
the  "arcing"  when  the  blade  bends. 

The  elastic  deflections  are  dropped  out  of  the  rigid-body 
equations  when  the  normal  mode  approach  for  the  elastic 
deflection  is  applied  (orthogonality  condition) .  They  are 
included  in  this  development  for  generality. 

To  define  the  position  vector  of  the  point  mass,  m^,  the  five 
coordinate  rotations  illustrated  in  Figure  6  are  required: 

(1)  Pylon  pitch  -  rotating  inertia  axis  system  by 

8p  about  E2  at  the  point  CA  gives  an^axis  system 
defined  by  unit  vectors  @i,  *2 t  and  ®2* 

(2)  Pylon  yaw  -  rotating  by  about  @2  at  the  point  CA 
gives  the  axis  system  defined  by  unit  vectors  @3,  @4, 
and  @2» 

(3)  Rotor  rotation  -  rotating  by  \p  about  @3  at  the  point  H 
(rotor  hub)  gives  the  axis  system  defined  by  unit 
vectors  65,  @5,  and  @3. 
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Figure  5.  Relationship  of  Discrete  Mass  Point, 

mj ,  With  Fixed  Coordinate  System  XPYFZ 


I*]'' 


Figure  6 .  Coordinate  Rotations 


(4)  Blade  flap  -  rotating  by  3  about  -fig  at  the  point 

H  gives  the  axis  system  defined  by  unit  vectors  -©g, 

@7,  and  ©g.  This  is  the  blade  local  coordinate 
system  which  blade  airloads  calculation  is  referred  to. 

(5)  Blade  pitch  -  rotating  by  0i  about  ©7  at  the  point  H 

gives  the  axis  system  defined  by  unit  vectors  @7,  ©9, 
and  @io •  is  the  blade  total  pitch. 

The  vector  defining  the  mass  point  position  is  given  by 


Ri  =  RH  ♦  rb 


in  which 


RH  X°E1  +  YoE2  +  ZoE3  +  ZH®3 


and 


r.  = 


A  /  A 


=  r.e?  +  yie9  -  6-e/  ♦  -  e.e?  - 


(1) 

(2) 

(3) 


where  e7 '  and 


are  referred  to  preconed  blade  position. 


To  work  consistently  in  the  rotating  system  the  coordinate 
rotations  shown  above  (and  in  Figure  6)  are  employed.  From 
inspection, 

e7  =  cosp  e5  +  sinp  e g 
eg  =  -  sinp  eg  +  cosp  eg 

©9  =  -  sinG^  sinp  eg  +  cos0^  eg  +  sinP^  cosp  eg 

E,  =  -  (sinG  cos\J  +cos0  sin\|  sin\A  ec 
1  \  P  P  P  7  5 

+  /sin©  sin\J'-cos0  sin\)'  cos\A  e£  +  cos0  cosifr  e0 

\  P  P  p  /  o  p  Yp  3 


n  *  A  A  A 

E2  =  cos^p  sin\J/  e^  +  cosij/p  cos\J/  eg  +  sin\pp  e^ 

Eg  =  |sin0p  sin\(/p  sint|r-cos9p  costyj  eg 

+  |sin0p  sin\J.p  cos\J/+cos0p  sinty)  e^  -  sin0p  cos^  eg 


(4) 

(5) 

(6) 


(7) 

(8) 


(9) 
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When  equations  4  through  9  are  substituted  into  1  through  3, 
the  mass  position  vector  can  be  written  as 


R.  =  R.  4C  +  R.  &,  +  R.  e, 
l  lx  5  ly  6  iz  3 


(10) 


in  which 


R^x  =  -  XQ  |sin9p  cosiJ/+cos9p  sintp  sintj  +  Y0  costp  sint 

+  ZQ  |sin0p  sintp  sint-cosGp  cos\f'j  +  r^cosp  -  +  f^jcosa 

-  'n^sinaQ  -  y^sinfK  sinp  (1 


(11) 


R.  =  Xr,  /sine  sin\}'-cos9  sin\j'  cost)  +  Y0  cos']  cost 
iy  ‘o  \  p  p  p  7  °  p 

+  Z0  ^sin9p  s^-n'Jp  cos\j+cos9p  sint)  +  |y^  cose^-f^  (12) 

R.„  ~  X„  cos9  cost  +  Yn  sint'  -  ZD  sine  cost'  +  zu  +  r.sinB 
izJ  °  p  p°  pu  p  pHir 


-  /6.  +€.\sina  t  n.  cosa  +  y.sinG.  cosB 
V  i  1/  o  *i  o-'i  l 


(13) 


At  this  point  an  approximation  for  the  shortening  displace¬ 
ments,  6^  and  e^,  must  be  made.  If  the  flapping  angle  due  to 
blade  flexibility,  Pif,  is  assumed  to  be  small  and  the  dif¬ 
ference  between  the  deflected  and  undeflected  radial  position 
to  be  negligible,  one  can  write  (see  Figure  7) 


p.f  *  sinpif  = 

1 


(14) 


Similarly 


Pif  ]_ 

6i  =  ^i  tanT~  =  pif  =  7r7 


e  .V 


(15) 


(16) 
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The  velocity  of  the  mass  point,  mi,  is  obtained  by  differentiat¬ 
ing  equation  10  with  respect  to  time.  Thus 


DR .  9R . 

=  -5T  =  -5T  +  “  x 


(17) 


where 


or 


— ♦  •  A  •  A  •  A 

co  =  0E_  +  \l'eo  +  \t'eo 

p  Z.  p  ^  *3 


(18) 


(iT  =  ^  0p  cosit'p  sin\}'-\i/p  cos\|/)  e^  +  ^  cos\J»p  cosiJ'+il'pSin^ j  e^ 

+  (®p  sln4-p+i)43 


(19) 


or 


“  =  “x45  *  “y46  *  “z43  <20> 

D  a 

in  which  ^  and  denote  the  time  rate  of  change  of  a  vector 

referred  to  an  inertial  and  a  body  coordinate  system 
respectively. 


Figure  7.  Geometrical  Relation  Between  Foreshortening 
and  Generalized  Coordinate. 
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5W1 


r . ti  ~i~  r  r-^--  IWiMMWWi  irti,m’T4Cv<M*'‘Ai-*i*»-,(J*tV*i'**z*4Aek**-'M 


Performing  the  mathematical  manipulation  indicated  in  equation 
17,  the  result  is 


—  A  A  A 

v.  =  v.  ec  +  v.  ec  +  v.  e„ 
1  lx  5  ly  6  iz  3 


(21) 


where 


v .  =  -  r 

IX 


[0  cos\l/p  COS\|'  +  \J/ 


\j/p  sin\|() 


L  sinp  |p-  {l 

-  jsinS^  cosp|(3-|©p  cos\|p  cosf  +  ^  sintyjj 

+  cos^  sinp  +  sin^p  +  4f^j| 

-  XQ  ^cos0p  sin\J<p  sinij'  +  sin0p  cosij'j  +  Y0  cosily  simji 
+  ZQ  |sin0p  sin^p  sin^  -  cos0p  costJ'J 

+  zh(®p  cos^p  cosij/  +  \{/p  sin\J/j 
+  t)^  cosaQ^0p  cos>J/p  cosij;  +  \jp  sin\J/j 
+  (©p  sin'J/p  +  ^ )  -  ri£  sinaQ 

-  '277  (^i2  +  ^i2)sinao(®p  cos^p  cost]/  +  \j'p  sin'J/j 

-  —  /ri.f).  +  £.£.\cosa 
r.  (  i  i  “i“i/  o 


(22) 


v.  =  r . 
iy  i 


cosp  |©p  sin\J/p  +  -  sinp  ^9p  cos^  sin>|/  -  ^  cosi|/^j 

-  y^  sin5^  +  sinp  ^Gp  sin’|/p  +  i|/|  +  cosp  ^9p  cosily  sin\(/  -  ^  cosj/jj 

-  XQ  |cosGp  sinij/p  cos\|/  -  sinBp  sin\J/)  +  Y0  costj/p  cos>J/ 

+  Z0  /sin0  sin')/  cos\|r  +  cosG  sin>]/) 

-  z1?(0  cosij /  siml/  -  vjf  cosij/^ 

H\p  p  p  I 

sina  ( 0  sinxp  +  +  cosa  (&  cosiJ/  sinil/  -  J'  costJH 

o\p  pl  o\p  p  p  / 


-  7 


-  I . 
i 
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+  ‘277  ^i2  +^i.2)  [sina0  (®p  cos'l,p  sin'1'  "  cos'l') 

-  cosa  ('J*  +  G  sim1/ 
o  \  p  l 

v.  =  X0  cosG  cos»j'  +  Y0  sin*  -  ZQ  sinG  cos\J/ 

12  P  p  w  p  u  p  p 

+  r^  cosp  |[?-  ^Gp  cos\|/p  cosil/  +  Jp  sin\J/| 

+  y.lcos£.  ?.  cosG  +  /g  cosijf  siri'l/  -  \ii  cos^/l 
•'ll  i[i  r\p  P  P  I 

-  sin^  sinp|(3-|§p  cosily  cosv!'  +  ^  siniljjj 

+  ri.  sina  (£  cosil/  cosil/  +  U<  sin>l»\ 
i  o  \  p  p  p  / 

-  £.  (g  cosij/  sirul/  -  *5  cost]/)  +  p.  cosa 

i  \  p  p  p  /  i  o 

i  /  2  2\  • 

+  |r)£  +  )  cosa0  (®p  cos'!,p  cos\J/  +  ip  sinij/j 

— —  fa.  p.  +  i-iA  sina 
ri  \  1  'i  ^1^1/  o 


cosa. 


(23) 


(24) 


The  moment  of  momentum  theorem  is  now  applied.  If  "pi  denotes 
the  linear  momentum  and  v£  the  velocity,  I^i,  the  angular 


momentum  of  a  mass  particle  m£  on  the  master  blade,  by 
definition  is 

"p .  =  m .  v\ 

Ki  li 

(25) 

=  R^  x  p^ 

(26) 

where  R^ 

has  been  defined  by  equation  10.  Then 

°ri  -  d  — 

"DF  =  Ri  x  Dt<mivi> 

(27) 

or 

dl\  ^  Dv\ 

_  =  rn.R.  X 

(28) 
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For  a  system  of  particles, 


i 


(29) 


BS  _ 

where  denotes  summation  over  all  blades  and  N  is  the  total 
1 

moment  about  the  fixed  point  of  the  external  forces  acting  on 
the  system.  Differentiating  equation  21  with  respect  to 
time  gives 


Dt  “  at 


co 


X  V. 
1 


(30) 


or 

D^i  / 

(av. 

IX 

Dt  "  \ 

at 

+  I 

f9Viz 

{  at 

+  CO 


\a  /9viy  \ 

Viz  “  Viy  “zK  +  pT  +  Wz  Vix  ‘  Viz  axj 


v.  -  v„- 
ly  ix 


0)  1 

lx  yj 


or 


Dv  ■ 

a\  =  =  a-  e_  +  a.  e ,  +  a.  e 0 

l  Dt  lx  5  ly  6  iz  3 


(31) 


(32) 


in  which  aix>  ^ay,  and  a^z  are  the  components  of  the 
acceleration  of  a  blade  particle. 


Substituting  equation  32  into  28  and  using  equation  29,  the 
results  are 


iz 


a . 


IX 


iz 


IX 


(33) 


(34) 


(35) 
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The  components  of  the  total  external  torque  about  the  fixed 
origin,  0,  consist  of  aerodynamic  forces  and  moments,  gravita¬ 
tional  torque  and  moments  generated  by  inertial  loads  at  the 
rotor  hub.  Thus 


N 


+  NM  *  NC 


+  N, 


(36) 


where 


(37) 


(38) 


(39) 


=  ~Rh  x  (-T) 


(40) 


in  which 


gy  and  gz 
? 


is  the  displacement  vector  of  the  it*1 
blade  segment  aerodynamic  center  from 
the  rotor  hub 

is  the  aerodynamic  force 

is  the  aerodynamic  pitching  moment 

are  gravitational  acceleration  components 

is  the  inertia  load  of  all  blades 


Two  assumptions  are  now  made.  First,  it  is  assumed  that  the 
rigid-body  flapping  degrees  of  freedom  are  used  to  simulate 
the  flapping  motion  of  a  gimballed  semirigid  proprotor.  (The 
blade  elastic  degrees  of  freedom  are  used  to  simulate  hingeless 
and  fully  articulated  rotor  types.)  When  this  assumption  is 
made,  the  individual  blade  flapping  can  be  related  to  tip- path 
plane  motion.  For  example,  the  flapping  of  the  blades  of  a 
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four-bladed  rotor  can  be  expressed  in  terms  of  the  flapping  of 
opposite  pairs  of  blades.  Thus  the  number  of  flapping  degrees 
of  freedom  can  be  reduced  to  two. 


For  a  three-bladed  gimballed  rotor  the  individual  blade  flapping 
relationships  with  the  tip- path- plane  motion  are  more  complex. 

If  aQ  represents  the  blade  precone  and  a*  and  b'  orthonormal 
tip- path- plane  flapping  in  the  rotating  system,  the  individual 
blade  flapping  is 

PL  =  a0  +  a'  (41) 


and 


(42) 


(43) 


where  a*  is  the  flapping  of  the  master  blade  and  b’  is  the 
tip-path-plane  flapping  90  degrees  to  the  master  blade. 

This  assumption  results  in  the  number  of  rigid-body  equations 
being  reduced  to  three:  tip-path-plane  flapping  90  degrees  to 
the  master  blade,  master  blade  flapping  a',  and  proprotor 
rotation  iJj.  These  correspond  to  equations  33,  34,  and  35. 

The  second  assumption  is  that  the  blade  elastic  motion  will  be 
represented  using  normal  modes  obtained  from  a  vibration  analysis 
which  includes  the  rigid- body  flapping  and  rotational  freedoms. 
This  assumption  will  be  discussed  in  more  detail  below.  Here, 
the  implication  is  that  the  coupling  terms  involving  blade 
elastic  motion,  rj.  ,  £4,  and  0if  may  be  neglected  in  the  rigid- 
body  equations  of  motion. 

The  form  of  the  resulting  rigid-body  equations  of  motion  is 
given  at  the  end  of  this  section. 

b .  Blade  Elastic  Equations  of  Motion 

(1)  General  Formulation  -  The  use  of  the  normal  mode  approach 
to  represent  the  blade  elastic  motion  consists  of  employing 
Lagrange's  equation  plus  the  normal  mode  orthogonality  con¬ 
dition.  Lagrange's  equation  for  normal  coordinates,  £j,  is 
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where 


T»V  ,D 


where 


are  kinetic  energy,  potential  energy, 
and  dissipation  function,  respectively 

is  the  generalized  force,  defined  as 


is  the  displacement  vector  of  the  itl1 
element 


F.  is  the  external  force  acting  at  the 

element  (V  in  this  case  consists  of 
aerodynamic  forces  and  moments  as  well 
as  gravitational  and  inertial  forces) 


The  kinetic  energy  for  a  blade  has  the  following  form, 


T 


+  v 


iy 


+  V 


1Z 


in  which 


v .  v  • 
ix,  ly, 

and  v. 

1Z 


are  defined  by  equations  22,  23,  and  24, 
respectively 


NS 


is  the  number  of  discrete  mass  stations 
on  the  blade 


(45) 


(46) 


The  potential  energy  may  be  written  as 


N 


n 


(47) 


in  which 


men  and  are  effective  mass  and  natural  frequency 
w  of  the  nth  mode,  respectively 

n 


N  is  the  total  number  of  natural  modes 

under  consideration 
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The  energy  dissipation  arising  from  the  structural  damping  of 
the  blade  is 


D 


men  "n 


(48) 


where  t  is  the  equivalent  viscous  damping  ratio. 

The  kinetic  energy  expression  given  by  equation  46  contains 
the  generalized  coordinates  r)j_ ,  and  9-^f  representing  the 

element  out-of-plane ,  in-plane,  and  pitching  deflection, 
respectively,  and  their  time  derivatives.  In  terms  of  normal 
coordinates,  the  following  relationships  apply. 


N 


N 


N 


^  -E 


it]n  n 


4.  4> * r  £  e.~  =7  <t> . n  K 

i?n  n  if  /  -■  lOn  n 


n  n  n 

Successive  differentiating  with  respect  to  time  gives 


N 

r).  =y  <p.  t 

'l  /  -1  lqn  ^n 
n 


and 


N 


N 

4-  i 

^1  /  ^  i£n  n 

n 


N 


N 


if 


E  ^ien  «n 


n 


N 


T).  =  >  <J>.  £  K-  =7  4>.c  t 

'i  f  -j  ir]n  ^n  ^1  i£n  n 


gif  *ien 


n 


n 


n 


(49) 


(50) 


(51) 


Substituting  equations  45,  46,  47,  and  48  into  equation  44, 
using  expressions  49,  50,  and  51,  and  the  orthogonality  con¬ 
ditions  among  normal  modes,  the  following  equation  of  motion 
for  a  normal  mode  is  obtained 


m 


ej 


NS 


m 


r  j  t 

of  the  jth 

mode  has  the  form 

l  2 

2\ 

It  -  2  \  . 

(*ird 

+  *Uj  )  + 

Pd  +  miyi  f i@j 

|cosG^ 

1  It; j  *iej 

+  sinf5i  *iej)J 
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(52) 


(53) 


in  which  y^  is  the  distance  of  a  blade  segment  mass  to  the 
feathering  axis,  §i  is  measured  with  respect  to  the  feathering 
axis.  Only  blade  root  collective  pitch  and  blade  built-in 
twist  are  considered  in  9j_. 

The  Qjj  term  in  equation  52  represents  the  inertial 
forces  acting  on  the  blade  due  to  hub  and  rigid-body  blade 
motion.  The  expression  for  Qj.  is  extremely  long,  but 
to  illustrate  the  coupling  terms  included,  equation  54  is 
provided. 


o_  •  (x)  =  /  <j>  •  -sina  +  (ri .  <t> .  .  +  £  .  <J>  •  r  \ 

vIj  1T}J  O  r.  \  l  1T]J  i£ j/ 

•  1. 
l 

|  |r^sinp  +  r]^cosaoj cos'^  cosi  + 


•  ■ 

O  \!'  sin1  cos 
P  P  P 


-  0  i  cosi  sini  +  i  sin'1'  + 

P  P  P  P 


cos 


+  r^cosp  (5-  |p  cos^/p  cos']  +  i  sinij 

-  t]^cosao  |0p  cos ii'  cosi  +  ip  sinij 

+  X  (cos  6  sinij  sinij'  +  sinP  cosil) 
o  \  p  p  p  / 

+  X  (0  sinG  sin',1  sini  -  vj'  cos0  cosil'  sin11, 

o  \  P  P  P  P  P  P 

-  i  cos0  sin!  cos'!'  -  9  cos 0  cosil/ 

P  P  P  P 

+  i  sin0  sinij) 

P  P  ' 

+  Y  cosil  sinij'  +  Y  / vl<  sinil'  sinij  +  i!'  cos'!  cos  ) 
op  o\  p  p  p  / 

+  Z  (cos©  cosi  -  sin0  sini  sin'll 
o  \  p  p  p  / 

-  Z  (i  sin0  cosil/  sini  +  6  cos0  sini  sin'1' 

o\  p  p  P  P  P  P 

+  \f/  sin©  sini  cosi  +  0  sinP  cosi 

P  P  P  P 

+  i  cos 0  sini) 

P  / 
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+  z„  fe  costy  cosi|/  -  6  ti i  sirul/  cos']'  -  €>  i  cos\l  sinij/ 
H\  p  p  p  p  p  p  P 


•  • 


+  \I'  sin\J/  +  vl/  \l'  cos 

P  P 


*) 


+  4i(ep  sin*p  +  i)  *  £.(ep  sin'ip  +  6p  ip  cos»p) 


-  r).  sina 
'i  o 


TFT^i2  +  ^i2)sinao(‘e'i 


•  • 


costl'  cos\i'  -  0  \|/  siruj'  cos\!» 

P  P  P  P 


-6  \j/  cos\J  sin\J'  +  'j/  sin\J 

P  P  P 


+  ij:  \|;  cos\l' 

P 


) 


■  FrK^i  +  Mi)  sinao  (®p  co^p  COS']'  +  'Ip  siruj) 

1  /•  2  t  2\ 

-  —  ( r; -  +  £  .  cosa 

r±  V  i  1  /  ° 

+  v^z  cos\i  p  cos.!  -  ii'p  sin\l) 

-  viy  (*  +  sin^  p) 


(54) 


Equation  54  is  the  x  component  of  Qjj*  with  blade  mass  un¬ 
balance  and  blade  torsional  mass  moment  of  inertia,  assumed  to 
be  zero.  The  hub  motion  coupling  terms  are  those  involving 
X0,  Y0,  ZG,  Gp»  and  v|/p  and  their  derivatives.  The  blade  Coriolis 
terms  are  those  containing  a0  and  (3.  (The  blade  flapping 
acceleration  terms,  j?,  have  been  dropped  because  of  the 
orthogonality  condition.  In  this  formulation  the  rigid-body 
flapping  is  considered  to  be  the  first  normal  mode.) 
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Q_,.f  the  external  force,  has  the  following  form: 
fcj 


sinP  cosij 
P  P 


QEJ  {[("’inj  *  h  cos5i  *iej)co 

-  FT^iru  +  «i)sinao]  [Fl  -  ",iSzsi 

-  n>igy  Sin*pj  +  |(PiTlj  +  Vi  cosB.  tiej)sina0  ♦  ^-(*iT,j  n- 
4^)cosaQ|  jf^  +  m^^cosPp  cos\J  -  sinP^  siml'^  siml) 

*Uj  *  Fi  sin«i  *iej][Hi  -  mi<hlsi 

*  cosep  sin+)J  ♦  Qi  4>iej  | 


+  $  .  r  . 
i4j 


[sinO  sin\J/  cos<j' 
P  P 


in  which 


Fi,  fj. 
and 


are  components  of  the  aerodynamic  force  in 
shaft  coordinate  system 


is  aerodynamic  pitching  moments  acting  at 
the  blade  elements 


§2,  6y 


are  components  of  gravitational  acceleration 
along  coordinates  of  Zp  and  -Yp  respectively 


(2)  Rigid-Body  Coupling  Considerations  -  Application  of  the 
normal  mode  approach  is  complicated  by the  fact  that  the  blade 
root  boundary  condition  is  variable,  depending  upon  the  loading 
of  the  blades.  For  example,  the  blade  reaction  to  forces  acting 
in  the  plane  of  the  rotor  depends  upon  whether  the  forces  pro¬ 
duce  a  net  torque  about  the  shaft.  If  a  torque  is  produced, 
the  shaft  will  rotate  and  the  blade  root  constraint  must  be 
considered  as  being  pinned.  On  the  other  hand,  if  the  net 
torque  is  zero,  then  the  blade  root  condition  must  be  considered 
as  being  cantilevered.  Similar  reasoning  can  be  applied  to  the 
out-of-plane  response. 

In  this  program  the  blade  root  boundary  consideration  is 
resolved  by  considering  the  blade  response  as  being  either 
polar  symmetric  or  polar  asymmetric  about  the  shaft.  This  is 
essentially  equivalent  to  the  usual  symmetric -asymmetric 
assumption  applied  to  beams  symmetrical  about  a  center  axis, 
e.g.  elastic  wing  response  about  the  longitudinal  axis. 
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The  polar  symmetric  motion  is  commonly  referred  to  as  collective 
motion  of  the  blades;  the  polar  asymmetric  as  cyclic  motion. 
Figure  8  illustrates  this  definition  of  collective  and  cyclic 
modes.  These  names  arise  from  the  blade  elastic  response  to 
collective  and  cyclic  pitch  respectively.  It  can  be  shown  that 
any  rotor  type,  rigid  (hingeless)  gimballed  or  teetering  semi¬ 
rigid,  or  fully  articulated  will  respond  to  external  excitation 
in  a  combination  of  collect i\je  and  cyclic  motion,  depending  upon 
the  nature  of  the  excitation. 


For  the  collective  modes  the  normal  mode  representation  is 
simple.  One  normal  mode  is  used  to  represent  each  collective 
mode.  For  a  three  bladed  rotor  the  collective  mode  normal 
mode  equation  is  written  as 


♦  2t“3'S  *  "A j  *  E  (?Ej  *  «ij)| 

eJ  |bl=l 


(56) 


The  normal  coordinate  for  each  individual  blade  is  then 
given  as 


*1  =  =  ^3  = 

1  J  Collective 


(57) 


A  more  complex  relationship  exists  for  the  cyclic  modes.  In 
that  case,  the  two  normal  modes,  having  the  same  natural 
frequency  and  mode  shape,  must  be  used  to  represent  the  cyclic 
response  of  the  three  blades.  This  follows  the  same  reasoning 
as  employed  for  the  rigid-body  flapping  and  is  illustrated  in 
Figure  9. 


If  the  modes  are  designated  mode  a  and  mode  b,  the  normal  mode 
equation  for  a  cyclic  mode  becomes 


I^ka  *  21  “k  «ka  *  V  «ka) 

_  1  1 
"  mekl 

QIkl  +  ^Ekl  “  "2  (*^Ik2  +  <?Ek2)  "  ?(^Ik3  +  ^Ek3) 

and 

4(<kb  +  2t  a 

k  ^kb  +  “k  ^kb) 

1 

"  mek 

■r(QIk2  *  «Ek2)  -  *  9Ek3)J 

(58) 


(59) 
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(b)  About  x  Axis 


Figure  9 


Components  of 
BLaded  Rotor. 


The  normal  coordinates  for  each  blade  are  then  given  as 


«kl  =  ^ka  (6 

^k2  =  “  ¥ka  +  IT^kb  (6 

^k3  =  "  I^ka  "  T^kb  (6 

Table  I  gives  the  boundary  conditions  which  must  be  used  to 
determine  the  collective  and  cyclic  modes  for  the  Feather-Step 
Analysis.  Note  that  if  the  analysis  is  being  used  to  simulate 
blade  folding  individual  blade  modes  are  employed. 


TABLE  I. 

BLADE  BOUNDARY  CONDITION  REQUIREMENTS 

Collective 

Cyclic 

Rotor  Type 

Vertical 

Inplane 

Vertical  Inplane 

★  + 

Fully  Articulated  ’  P 

P 

P 

P 

Offset  Flapping 

Hinge  P 

P 

P 

C 

Simballed 

C 

P 

P 

C 

Hinge less 

C 

P 

C 

C 

Blade  Folding 

C 

C 

C 

C 

P  -  Pinned 

C  -  Cantilevered 

*  -  The  first  normal  mode  is  rigid-body  flapping. 

+  -  The  blade  lead-lag  motion  must 

be  represented  using  a 

normal  mode. 

c .  Proprotor  Aerodynamics 

Two  dimensional,  quasi-steady  strip  theory  is  used  to  calculate 
the  airloads  acting  on  the  blades. 1®  Each  blade  is  divided  into 
strips  centered  on  a  mass  point,  as  illustrated  in  Figure  10. 

Up  to  ten  strips  may  be  used.  The  width,  chord,  chordwise 
location  of  the  aerodynamic  center,  and  the  built-in  twist  is 
specified  for  each  individual  strip. 

The  local  velocities  at  a  blade  element,  i,  are  first  determined 
by  adding  the  free- stream  and  induced  velocity  to  the  mass  point 
velocity  vector,  (equation  21).  If  represents  the  _ 

induced  velocity  perpendicular  to  the  tip  path  plane,  and  U  the 
free  stream  velocity,  the  total  velocity  at  a  mass  point  on  a 
blade  can  be  written  as 
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Figure  10.  Blade  Aerodynamic  Strip  Relationships. 
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v.  =  -  v.  e,  -  v.  ^  *  v,-,  -  U  1 

l  lx  5  ly  6  iz  3  1 

-  a. cosa  eQ  -  v.sina0  e7 
where  ac  is  the  rotor  precone. 

The  radial  tangential  and  vertical  components  of  the  total 
velocity  can  then  be  obtained  by  forming  dot  products  of 
with  unit  vectors  £7,  -&6»  ar*d  eg,  respectively.  However,  the 
radial  flow  will  be  neglected.  The  other  two  components  take 
the  following  forms: 
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in  which  ecf  is  the  distance  from  the  blade  quarter  chord  to 
its  center  of  gravity,  positive  if  center  of  gravity  is  on  the 
leading  edge  side.  The  terms  containing  ecf  are  the  components 
of  the  relative  velocity  between  the  center  of  gravity  and  the 
quarter  chord. 

The  strip  total  velocity,  inflow  angle,  and  Mach  number  are 
then  obtained  from  the  following  relationships: 


'i =  f 


V_.  =  'vTi2  ♦  vPi2 


) 


1/2 


(66) 


-1  pi 
*i  =  tan  7“ 

Ti 


(67) 


Mi  =  T 

where  S  is  the  speed  of  sound. 


(68) 
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The  strip  angle  of  attack  is  then  calculated  as 

a.  =  0.  +  4).  (69) 

ill 

where 

=  6c  -  cos\j/  -  sini}/  ♦  -  aQ|tan6^  (70) 

in  which 

0C  is  blade  root  collective  pitch  angle 

is  the  lateral  cyclic 
is  the  longitudinal  cyclic 
0Ti  is  the  blade  built-in  twist 

©if  is  the  blade  flexible  pitching  deflection 

63  is  the  negative  pitch-flap  coupling 


The  aerodynamic  lift,  drag,  and  pitching  moment  at 
chord  of  the  strip  are  then  given  by 

the  quarter 

dLi  ■  ici?civi2(dri) 

(71) 

dDi  =  -iV°ivi2(dri) 

(72) 

dMi  =  •JCMpCi2vi2  (drl) 

(73) 

where 

c^  is  the  blade  chord 

dr^  is  the  segment  length 

CL  is  the  lift  coefficient  for  the  strip 

CD  is  the  drag  coefficient  for  the  strip 

is  the  pitching  moment  coefficient  for  the  strip 

The  lift,  drag,  and  pitching  moment  coefficients  for  the  strip 
are  obtained  from  tabular  forms,  where  the  coefficients  are 
given  versus  angle  of  attack  and  Mach  number.  Figure  11  is  an 
example  of  the  tabular  form  employed  in  the  computer  program. 

Note  that  reverse  flow  and  stall  are  accounted  for  in  the  tables. 
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MACH  NUMBERS  ll»  ANGLES  OF  ATTACK  3) 


ANGLES 

OF 

MACH  NUMBERS 
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0.40000 
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Figure  11.  Tabular  Form  of  Aerodynamic  Coefficients. 
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The  pitch  damping  of  the  strips  is  approximated  using  an 
approximation  which  includes  the  effect  of  stall. H  The  pitch 
damping  moment  is  given  by 


(74) 


where  is  the  rate  of  change  of  the  strip  angle  of  attack  and 
Gft.  is  a  coefficient  obtained  from  the  table  shown  in  Figure  12. 
Note  the  C&.  is  a  function  of  the  strip  angle  of  attack  af,  and 
a  reduced  velocity  Kf,  defined  as 


(75) 


where  o>g  is  the  frequency  of  a  selected  torsional  mode  of  the 
blade;  usually  the  fundamental  torsional  mode. 

The  strip  thrust,  inplane,  and  radial  forces  and  pitching 
moment,  in  the  blade  rotating  coordinate  system  x,  y,  z  are 


then  given 

as 

Fi 

= 

^dL^  cos4>^  ♦  dD^  sin^^jcosp 

(76) 

Hi 

= 

(-dl^  sin<j>^  ♦  dD^  cos<{k  j 

(77) 

fi 

= 

(dL^  cos<}>£  +  dD^  sin^^jsinp 

(78) 

dQ. 

= 

dM.  +  dM.  ’  +  e  .  (dL.  cosa.  ♦  dD.  sina.) 

l  l  ei  \  l  l  i  J./ 

(79) 

where  sci  is  the  distance  from  the  blade  aerodynamic  center 
to  the  blade  elastic  axis.  Positive  if  aerodynamic  .M&ntcr  is 
ahead  of  the  elastic  axis. 

3.  Wing -Pylon- Hub  Equations  of  Motion 

a.  Dynamics 

The  normal  mode  method  is  again  used  to  simulate  the  wing - 
pylon-hub  dynamics.  As  with  the  blade  elastic  equations,  this 
greatly  simplifies  the  development  of  the  equations  of  motion. 

Figure  13  shows  the  math  model  for  the  wing-pylon-hub  system. 
The  wing  is  represented  by  lumped  masses  located  spanwise 
along  the  wing  and  having  inertia  and  unbalance  properties. 
Beamwise ,  chordwise,  and  torsional  deflections  of  the  wing  are 
accounted  for,  as  well  as  pylon  pitching  and  yawing  with 
respect  to  the  wing. 
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NOTES : 

The  distance  from  the  pylon  conversion  axis  to  the  pylon 
center  of  gravity  is  defined  as  Zp.  The  signs  of 
eccentricities,  Xp,  Yp,  as  well  aS  Zp,  are  determined  in 
the  XYZ  coordinate  system.  The  distance  from  the  pylon 
conversion  axis  to  the  hub  is  defined  as  z^. 


Figure  13.  Math  Model  for  Wing-Pylon-Hub  System. 
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The  pylon  is  considered  to  be  inelastic  with  the  pylon  center 
of  gravity  located  a  distance  Zp  along  the  shaft  axis  from  the 
conversion  axis,  and  with  eccentricities  Xp  and  Yp  from  the 
shaft.  The  hub  is  treated  as  a  point  mass  located  on  the 
shaft  centerline,  a  distance  z^  from  the  conversion  axis. 

By  using  the  normal  mode  approach  only  one  derivation  is 
required,  the  generalized  forces  for  the  normal  modes.  These 
generalized  forces  consist  of  the  wing  airloads,  gravitational 
forces,  and  the  proprotor  aerodynamic  gravitational  and 
inertial  forces  and  moments  acting  at  the  hub. 

By  definition  the  generalized  force  for  the  jth  normal  mode  is 


where  is  the  external  force  acting  at  the  ith  element. 

The  external  force  applied  at  the  ith  segment  on  the  wing  is 


F  .  =  - 

Wl 

HwiEl  +  ^wi  E2 

-  (Twi  -  mwigz)E3  -  mwigy  E2 

in  which 

Twi 

is  the  lift 

force 

Hwi 

is  the  drag 

force 

^wi 

is  the  pitching  moment 

8z^dd  gy 

are  components  of  gravitational  acceleration 

is  the  wing 

discrete  mass 

(45) 


(80) 


The  external  force  applied  on  the  pylon  is 


Fp  =  VzS  -  mpgY  E2 


where  m^  is  the  pylon  mass. 

The  external  force  applied  on  the  pylon  tip  is 

"Fh  =  -  TFFX  &2  +  TFFy  e4  +  TFFZ  §3  -  TFMX  &2 

+  TFMy  e^  +  TFMZ  e3  +  mHgzE3  -  mHgy  E2 


(81) 


(82) 
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TFFV  =  FF„  +  Wv  +  F  cos  -  F  siml» 
X  X  X  x  y 

TFF„  =  FF„  +  Wv  +  F  sin'!1  +  F  cos 
Y  Y  Y  x  y 

TFFZ  =  FFZ  *  Wz  ♦  Fz 

TFMX  =  (Nx  *  Mx)  cosi!'  -  |  +  Myjsiruj/ 
TFMy  =(nx  +  Mx)sin:  +  (n  +  My)cos*j/ 


TFM„  =  N  +  M 

Z  z  z 


in  which 
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(84) 

(85) 

(86) 

(87) 

(88) 


FFX’  FFY’  ffz 
WX’  WY’  WZ 


The  position  vectors  for 
as 


are  airloads 

are  gravitational  forces 

are  inertial  forces 

are  aerodynamic  and  gravitational 
moments 

are  inertial  moments 

the  wing,  pylon,  and  hub  are  given 
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(89) 


Rp  ■  (xo  +  zp  <*»*p  cos0p  -  Xp  sinep  -  Yp  sin0p  oos0p)E1 
*  (Yo  +  zp  +  Yp  =°=*p)E2 

+  (z  -  z  cos^  sin0  -  x  cos0  +  Y  sini^  sin0  \£L(90) 
\°p  p  PP  P  P  p  P/  3 
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irH  -  (x0  +  zH  coai|jp  oosepjij 

+  (Yo  +  ZH  sin4'p)S2 

+  (z0  -  zH  cos*p  sinep)i3 


(91) 


where 

XQ,  YQ,  and  Z0  give  the  position  of  the  conversion 

axis-shaft  axis  intersection 


0 


P 


equals  0pc  +  6pp  with  0pp  as  the  pitch 
due  to  wing  torsional  deflections 


The  generalized  coordinates  are  now  replaced  by  the  orthonormal 
mode  shapes 

X0  “  S^XAnqn 
Y0  “  S^YAnqn 

Z0  ^^ZAnqn 

%  "  EVnqn  (9 

6PP  "  2*6Anqn 
^wi  2!)^Zinqn 
^wi  *  ^Xinqn 
®wi  *  ^0inqn 

in  which 

^XAn  is  t^ie  xf  comPonent  of  nth  coupled  mode 

at  the  CA 
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4>x^n  is  the  Xp  component  of  the  nth  coupled 

mode  at  the  ith  segment  on  the  wing 

qn  is  the  nth  normal  coordinate 

By  the  same  token  as  to  derive  equation  52,  the  equations  of 
motion  for  the  wing-pylon-hub  system  can  be  written  as 


ij.  ♦  2™^  +  -J*  Oj 

The  generalized  force,  Qj ,  has  the  following  form 


(93) 


w  . 

Qj  "£p*i*eij  _Hwi*xij  Twi  ^zi j 
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+  TPFX  (z„  4>eAj  -  sinepc  -  *ZA .  cos8po) 
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(*XAj  cos0pc  '  *ZAj  sin6pc) 

TFFZ  (ZH  *<iAj  +  *VAj)  *  TFMZ  *0Aj]  £*>An  9n  *94> 


The  generalized  mass,  Mgj,  may  be  obtained  directly  from  the 
vibration  analysis  used  to  generate  the  normal  modes,  or  may 
be  calculated  (as  is  done  in  ARAP06)  using  the  relationship 
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+  ^ 2  ZAj  +  z2H^%Aj  +  ^ 2  0Aj) 
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in  which 
I 


is  the  wing  torsional  mass  moments-of-inertia 
about  its  elastic  axis 


is  the  pylon  pitching  mass  moment-of-inertia 
about  the  pylon  center  of  gravity 
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b.  Aerodynamics 

Two-dimensional  steady  strip  theory  is  used  for  the  wing  air¬ 
loads.  Figure  14a  shows  the  wing  strip  arrangement.  Pylon 
and  hub  airloads  are  neglected. 


The  normal  and  tangential  components  of  the  relative  wind 
acting  at  the  wing  quarter  chord  are  obtained  by  inspection 
(neglecting  higher  order  terms) 


‘wi  wi 


cosa 

w 


(96) 


vm.  =  u  +  £  . 

Ti  swi 


“  Vwi  6wi  sinaw 


(97) 


where  aw  is  the  wing  steady  angle-of-attack 
The  inflow  angle  is 


^wi 


*  tan 


!pi 
vTi 


with  the  local  angle-of-attack  given  by 
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wi 
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wi 


+  i  +  ot 

*wi  w 


and 


The  local  lift  and  drag  force  are  then,  respectively 


Lwi  -  Vcwi  V  4Li 

Dwi  -  h  CDpcwi  SLi 


in  which 


Cwi 

is 
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is 
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-H 

is 

the 

wing  chord 
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total  relative  velocity 


(98) 

(99) 

(100) 

(101) 
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CL  and  CD  are  the  lift  and  drag  coefficients, 
defined  as  follows 


a  a 


wi 


C  =  b,  +  b.a  .  +  b.ot  . 2 
D  1  2  wi  3  wi 


(103) 

(104) 


The  local  thrust,  in-plane  (flight  path  plane)  drag,  and 
pitching  moment  are  then  (see  Figure  14b) 


Twi  -  Lwi  cos*wi  +  Dwi  sin»wi 

(105) 

Hwi  =  -Lwi  s*n<^wi  +  Dwi  cos<^wi 

(106) 

%i  *  Ewi  (Lwi  cos“wi  +  Dwi  sin“wi) 

(107) 

where  ew^  is  the  distance  from  the  wing  elastic  axis  to  its 
aerodynamic  center  (AC) ,  considered  positive  if  the  AC  is 
toward  the  leading  edge. 

4 .  Wing-Blade  Aerodynamic  Interference 

The  aerodynamic  interference  between  the  wing  and  proprotor 
blades  is  simulated  by  determining  the  change  in  blade  airloads 
caused  by  passage  of  the  blade  through  the  wing  upwash  field. 
The  influence  of  the  blade  downwash  on  the  wing  is  neglected. 
The  problem  is  treated  in  two  parts ,  the  approximation  of  the 
flow  field  ahead  of  the  wing  and  the  location  of  the  blade 
elements  in  that  flow  field. 

a.  Approximation  of  Flow  Field  Ahead  of  the  Wing 

The  flow  field  of  the  wing  is  approximated  by  superimposing 
the  flow  field  of  an  elliptical  profile  at  zero  angle  of 
attack  with  that  of  a  flat  plate  having  the  wing  angle  of 
attack.  The  ellipse  simulates  the  wing  thickness  and  the 
flat  plate  the  circulation.  Figure  15  illustrates  this 
representation . 

(1)  Flow  Field  of  an  Elliptical  Profile  -  Defining  Z-plane  as 
the  math  plane,  £ -plane  as  the  physical  plane,  as  shown  in 
Figure  16,  the  streaming  motion  past  a  circular  cylinder  of 
radius  Ri  in  the  Z-plane  is  given  by  the  complex  velocity 
potential 
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(a)  Flow  Problem 


ELLIPSE 


(b)  Math  Model 


Figure  15.  Approximation  of  Flow  Field  Ahead  of  the 
Wing. 
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Figure  16.  Conformal  Mapping  of  an  Ellipse  to 
a  Circle. 

The  velocity  field  in  the  c-plane  is  given  by 


n  _  d$  _  d4>  dZ 

~  a?  ~  37  a? 


(109) 


Let  the  Joukowsky  transformation  be  the  mapping  function, 
which  has  the  following  form 

R  * 

C  -  Z  +  (110) 

where  R2  is  a  parameter. 

The  inverse  transformation  is  then 


(111) 
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To  map  an  ellipse  of  thickness  ratio  6  in  the  c “Plane  to  a 
circle  of  radius  R^  in  the  Z-plane  by  using  Joukowsky's 
transformation,  the  following  relationships  are  established 

Rx  =  h(l  +  6)  (112) 


R2 


(113) 


Differentiating  equation  108  with  respect  to  Z,  and  differ¬ 
entiating  equation  111  with  respect  to  c,  substituting  the 
results  in  equation  109  the  flow  field  of  an  elliptical  profile 
at  zero  angle  of  attack  is  then  obtained 


U 

Ce 


u 

00 


(2  +  Rl)  (2 
(2  +  R2)  (2 


(114) 


where  Z  is  a  function  of  z,  defined  by  the  inverse  transforma¬ 
tion  equation  111. 


The  complex  perturbation  velocity  is  calculated  by 


AU 


Ce 


u 


5e 


-  u 


(115) 


(2)  Flow  Field  of  a  Flat  Plate 

For  a  flat  plate  with  an  angle  of  attack,  the  complex  velocity 
potential  is  obtained  by  superimposing  the  following, 


(a)  uniform  flow  along  x-axis  with  a  doublet 

(b)  uniform  flow  along  y-axis  with  a  doublet 

(c)  a  line  vortex  with  clockwise  circulation  at  the 
origin 


Thus 


*  (Z) 


■  --(■  *  ¥)  -  -  ¥)  *  e 


(116) 


where 


and  v^  are  the  x  and  y  components  of  the 

freestream  velocity 

r  is  the  circulation,  positive  if 

counterclockwise 
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If  the  Joukowsky  transformation,  equation  110,  is  used  to 
map  a  flat  plate  of  chord  length  c  in  the  5-plane  to  a  circle 
of  radius  R^  in  the  Z-plane,  the  following  relationships  are 
obtained 


R 


1 


(117) 


The  mapping  function  for  the  case  of  a  flat  plate  can, 
therefore,  be  written  as 


c  -  z  +  sL 

16Z 


(118) 


Thus 


(119) 


Differentiating  equation  116  with  respect  to  Z  yields 


d$> 

dZ  "  u°° 


(‘ 


iv 


00 


iT  1 
2  7T  Z 


(120) 


With  equations  119,  120  derived,  the  flow  field  of  a  flat 
plate  having  angle  of  attack  can  be  obtained  as 


The  Kutta- Joukowsky  condition  requires  that  the  value  of  the 
circulation  is  such  as  to  make  the  term  in  the  parentheses 
equal  to  zero  at  the  trailing  edge. 12  This  is  essential 
to  avoid  having  infinite  velocity  component  at  c  =  c/2. 
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Therefore , 


r  =  2'rrvoo^  =  TTVooC 


(122) 


Substituting  this  expression  into  equation  121,  for  a  flat 
plate  having  unit  chord  length,  we  have 


U 

CP 


u 

to 


(123) 


The  perturbation  velocity  due  to  a  flat  plate  with  angle  of 
attack  a  is  given  in  the  following  form 


AU 


CP  “ 


U 

CP 


+  iv  ) 
00 


(124) 


in  which 


c 

8 

II 

G 

cosctw 

(125) 

< 

8 

II 

G 

sinaw 

(126) 

Where  U  is  the  free  stream  velocity. 

b.  Calculations  of  Blade  Element  Locations  Ahead  of  the  Wing 

The  wing  upwash  velocity  is  vectorily  added  to  the  proprotor 
blade  element  velocity  to  obtain  the  change  in  angle  of 
attack  and  dynamic  pressure  caused  by  the  presence  of  the 
wing.  This  requires  the  location  of  the  blade  element 
relative  to  the  wing  be  known.  Figure  17  shows  the  top  and 
side  views  of  a  blade  element  location  ahead  of  the  wing. 

The  position  vector  of  a  blade  mass  element  mi  from  the  wing 
leading-edge  can  be  written 

q  =  d*i  +  ZHe3  +TL  (127) 


in  which 

di  is  a  vector  from  the  wing  leading-edge  to  the 

pylon  conversion  axis 
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Zjj@3  is  a  vector  from  pylon  conversion  axis  to 

rotor  hub 

"r ^  is  the  position  vector  of  a  blade  mass  element 

relative  to  the  rotor  hub,  which  is  given 
by  equation  3. 

In  terms  of  the  fixed  coordinate  system,  XFYpZp,  the  following 
is  obtained 


ZH®3  +7i  ■  qx/l  +  ?iyp®2  +  =iZpE3 


(128) 


in  which 


zH  cosi|)p  cos6p  +  r^  (-cosB  sini|;  sin^p  cos6p 


-  cos$  costy  sin0  +  sinB  cos4>  cose  \ 

P  P  P  J 


+  y^  sin0p^cos0^  sinij)  +  sin§^  sin$  cosij;) 


+  cos0  sin0,  cos$  cosij; 

pi  p 


+  sin^p  cos0p|sin0^  sinB  sinij>  -  cos0^  cosi^ 
-  ^-cosaQ  simp  sini^p  cos0p  -  cosaQ  cosi|;  sin0p 

+  sinaD  cos^p  cos0pj 

+  ni  (cosaG  cosi^p  cos0p  +  sina0  sini^  sintjjp  cos0p 

+  sinaQ  cos\p  sin0pj 


-  (sinij;  sin0p  -  cos\p  sinij/p  cosep) 


(129) 
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■  zH  sinipp  +  ^cosB  sirup  cosipp  +  sinB  sinipp) 

+  y^  (cos0^  cosip  cosipp  +  sin0^  cosB  sinipp 

-  sin0^  sin8  sinip  cosip^j 
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i  \  o  p  o  T  rp/ 

“  2rj(ni2  +  5i*)  (cosao  8in^  cos^p  +  sinao  sin^p) 

-  5^  cosip  cosipp 

«  -zH  cosipp  sin0p  +  (cosB  sinip  sinipp  sin0p 
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P  P  P  / 

+  cos0p  |cos0^  sinip  +  sin^  sinB  cosip) 

-  sinS^  cos8  cosipp  sin0p 

+  sinipp  sin0p  (cos0^  cosip  -  sine^  sinB  sinip) 

-  2r~(ni2  +  ^i2)  (cosao  sin^p  sin0p  "  cosa0  cos^  cos0p 
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(130) 


(131) 


A  A 


If  e^,  er,  represent  the  unit  vectors  in  the  wing-chord  co¬ 
ordinate  system,  referring  to  Figure  17,  we  have 


di  - 


+  (?iYF  +  YCa)  tanA 


h  -  'iVp  E2  <132> 


substituting  equations  128  and  132  into  equation  127,  £i  is 
given  in  the  wing-chord  coordinate  system  as 


*i  = 


L  -  ? 


cosaw  +  z,in  sina 


’iZ. 
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(133) 


The  two  components  in  er  and  £  directions  are  used  to  locate  a 
blade  element  ahead  of  the  wing  in  the  ?-plane.  With  these 
components  calculated,  the  flow  field  is  then  obtained  by 
using  equations  114  and  123. 


5 •  Special  Provisions 

a.  Blade  Feather/Unfeather  Capability 

To  stop  or  start  the  proprotor  the  blades  are  feathered  or 
unfeathered,  respectively.  In  the  feather/stop  analysis  the 
blade  collective  pitch  may  be  specified  as  a  function  of  time; 
constant  at  time  tj_,  then  varying  linearly  to  a  new  value  at 
time  t2 . 

The  windmilling  blade  collective  pitch  for  a  given  rpm  may  be 
estimated  using  the  following  approximation: 


fi 


c 


cose 


pc 


(134) 


where 

6pC  is  the  steady  shaft  angle-of-attack 

U  the  airspeed  in  inches  per  second 

R  the  blade  radius  in  inches 

Q  the  proprotor  angular  velocity  in  radians  per 

second 

Gt  the  blade  twist  at  the  73  percent  radius 

(washout  negative) 

(This  formula  is  based  on  zero  angle  of  attack  at  73 
percent  radius.) 
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b .  Hi’b  Restraint  and  Flapping  Stops 

Hub  restraint  is  commonly  employed  with  gimballed  proprotors. 

The  analysis  includes  the  capability  to  simulate  non-isotropic 
hub  restraint.  This  is  accomplished  by  introducing  moments 
in  the  rotor  fixed  coordinate  system,  XYZ  about  the  X  and  Y 
axes . 

Mx  -  -Kb  Slat  (135) 


and 

mY  =  *a  Biong  (136) 

where 

Ka  corresponds  to  longitudinal  hub 

restraint 

corresponds  to  lateral  hub  restraint 

3 iong  corresponds  to  the  longitudinal 

tip-path-plane  flapping 

3iat  corresponds  to  the  lateral 

tip-path-plane  flapping 

The  components  of  these  moments  in  the  blade  rotating 
coordinates,  xyz,  are 

mx  =  Mx  cos^  +  My  sinij>  (137) 

my  =  -Mx  sini|j  +  My  costy  (138) 

These  are  applied  to  the  right  hand  side  of  the  b'  and  a' 
flapping  equations,  respectively. 

The  longitudinal  and  lateral  restraint  rates,  Ka  and  K^,  may 
be  specified  as  a  function  of  rpm. 

An  approximation  to  flapping  stops  is  also  provided  for  the 
gimballed  rotor.  If  the  flapping  of  any  blade  exceeds  a 
specified  magnitude,  a  flapping  moment  is  added  to  equations 
33  and  34.  For  instance,  if  the  second  blade  flapping  of 
a  three  bladed  gimballed  rotor  satisfies  the  following 
condition 
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8,  -  B  .  <  0 

p2  pstop 


The  flapping  moment  components  are  calculated 

(^2  ^stop) 


Amx  2  Kstop  (^2  ^stop/ 


Amy  ~  ~  7  Kstop  ^2  “  ^stopj 


where 


K 


stop 


is  very  large 


c.  Gust  Simulation 


(139) 


(140) 

(141) 


Gust  encounters  with  head-on  and  vertical  gusts  of  sinusoidal, 
sine-square  or  sharp  edged  shape  may  be  specified.  The  magni¬ 
tude  and  time  for  the  encounter  may  be  specified  and  for  the 
sine-square  gust,  the  gust  length.  The  magnitude  and  fre¬ 
quency  are  specified  for  a  sinusoidal  gust.  Symbolically, 
the  gust  may  be  written  in  the  following  form 


where 

Vc  is  a  head-on  gust 

yh 


V_  is  a  vertical  gust 

yv 


The  coordinate  transformations  discussed  earlier  in  this 
section  are  used  to  transfer  the  \f  into  the  velocity  vectors 
in  the  local  blade  coordinate  system  for  the  proprotor  and 
the  fixed  coordinate  system  for  the  wing. 

The  gust  is  assumed  to  be  encountered  simultaneously  by  the 
wing  and  proprotor. 

d.  Frequency  Response  Excitation 

The  frequency  response  of  the  system  analyzed  can  be  deter¬ 
mined  using  two  means  of  aerodynamic  excitation  (1)  cyclic 
pitch  and,  (2)  aerodynamic  vane  located  on  the  pylon.  (The 
vane  may  be  used  to  simulate  aileron  excitation.) 
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(1)  Cyclic  Pitch  Excitation  -  Tha  cyclic  pitch  excitation  is 
to  determine  the  position  of  the  swashplate.  If  As  (t)  and 
(t)  are  the  lateral  and  longitudinal  cyclic  pitch  components 

for  t  <  tT 

Ax(t)  =  Ax  (143) 

B^t)  =  B1  (144) 


for  t  ^  tT 


Ax(t) 

=  A1 

+  Alc  COS 

“ac  (4  ■ 

*t) 

+  Als 

sin 

ks  (*  *  *T) 

Bj  (t) 

=  B1 

+  Blc  COS 

“be  (‘  - 

*t). 

+  Bls 

sin 

“bs  (*  -  ‘t) 

where 


(145) 

(146) 


A1 

B1 

The  four  us 


is  the  time  excitation  initiated 
is  the  steady  lateral  cyclic 
is  the  steady  longitudinal  cyclic 
are  frequencies  of  excitation,  defined  as 

3u. 


,  _  ,  _ac /,  .  \ 

wac  ~  wacO  +  9t  \  ~  tT) 


9u 


u 


as  “  wasO 


as 


3 1 


(*  -  St) 


u 


be  “  wbcO 


^  -  *x) 


u 


bs  *  wbsO 


^  -  0 


follows 

(147) 

(148) 

(149) 

(150) 


(2)  Aerodynamic  Vane  -  With  the  location  and  dimensions  of 
the  exciter  on  the  pylon  having  been  specified,  the  vertical 
and  tangential  velocity  components  can  be  computed.  The 
angle  of  attack  of  the  vane  is  then  calculated  by 
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(151) 


°e  "  9p  +  “Oe  +  “Os  •in[“e(t  '  *e)]  +  tan‘‘  7^ 
in  which 

0p  is  the  pylon  conversion  angle 

vpv  is  the  vertical  velocity  component 

vTv  is  the  tangential  velocity  component 

oiQe  is  the  steady  angle  of  attack  of  the  vane 
relative  to  the  pylon 

oiqs  is  the  oscillatory  angle  of  attack  of  the  vane 
relative  to  the  pylon 

te  is  the  time  when  excitation  of  the  vane  begins 

we  is  the  exciting  frequency  given  by 

“e  “  “eO  +  "  fce)  (152) 

The  lift  curve  slope  of  the  vane  is  specified.  By  the  same 
token  as  wing  aerodynamic  loads  and  moments  are  treated,  the 
ones  for  the  vane  are  calculated  and  used  as  part  of  the 
wing  pylon  forcing  function. 

6 .  ARAP06  Computer  Program 

The  final  form  of  the  equations  of  motion  are  shown  in  Figure 
18 .  Up  to  nineteen  nonlinear  simultaneous  equations  are 
involved,  depending  on  the  degrees  of  freedom  selected  by 
the  user. 

The  equations  of  motion  are  solved  using  a  step-by-step 
integration  technique.  Details  of  the  method  of  solution 
are  provided  in  Volume  II.  The  output  data  consists  of  a 
time  history  of  selected  accelerations,  displacements,  and 
forces  and  moments.  If  desired,  the  time  history  may  be 
machine  plotted  (CALCOMP) . 

a.  Program  Capabilities 

The  ARAP06  program  has  been  structured  to  provide  the  user 
with  considerable  latitude  in  the  math  model  to  be  analyzed. 
Table  II  gives  examples  of  the  math  models  available.  (The 
math  model  is  specified  by  control  cards.) 
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Figure  18.  Form  of  Equations  of  Motion. 
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In  addition  to  the  calculation  of  proprotor  stability  and 
dynamic  response  during  the  FPR  feather/stop  cycle,  the 
ARAP06  program  can  be  used  to  calculate: 

(1)  Proprotor  aerodynamic  derivatives,  including 
frequency  response  characteristics. 

(2)  Oscillatory  forces  and  moments  acting  at  the  hub. 
These  data,  when  used  with  a  finite  element 
structural  analysis  such  as  NASTRAN,  provides 
useful  design  information. 

(3)  Loads  and  response  during  blade  folding.  In  this 
respect  the  ARAP06  program  has  a  more  general 
capability  than  the  response  analysis  specifically 
developed  as  part  of  the  Folding  Dynamics  Analysis. 

(4)  Dynamic  characteristics  of  tilt-rotor  7T0L 
including  proprotor/pylon  stability,  blade  flapping 
and  loads,  and  pylon  and  fuselage  vibration. 

b.  Program  Usage 

(1)  Data  Requirements  -  An  input  format  and  us'  rs  guide  is 
given  in  Volume  II.  The  direct  data  requirements  consist 
of  the  proprotor  and  wing/pylon/hub  system  inertial  properties, 
geometry,  and  normal  mode  frequencies  and  mode  shapes.  The 
wing  lift  curve  slope  and  proprotor  hub  restraint,  if 
applicable,  must  also  be  supplied. 


Indirect  input  data  requirements  consist  of  the  stiffness 
properties  of  the  wing  and  blades  which  are  used  in  vibration 
analyses  to  generate  the  normal  modes.  The  tabulated  2-D 
airfoil  data  currently  incorporated  in  the  program  is  satis¬ 
factory  for  most  dynamics  work,  but  the  user  can  substitute 
his  own  tables  directly  in  the  source  deck  if  desired.  These 
tables  currently  include  the  following  airfoils: 


Airfoil 


Applies  to  Radial  Station 


t/c  *  0.25 
t/c  -  0.18 
t/c  -  0.12 
t/c  «  0.08 


0.075  <  r/R  <  0.45 
0.45  <  r/R  <  0.70 
0.70  <  r/R  <  0.90 
0.90  <  r/R  <  1.0 


(2)  Usage  Techniques 

(a)  Proprotor  Modeling  -  As  indicated  in  Table  II,  the  ARAP06 
program  can  be  used  to  model  a  number  of  rotor  types.  The 
gimballed  rotor  type  is  automatically  modeled  if  the  gimballed 
rotor  option  is  selected.  To  model  a  hingeless  rotor  (or  a 
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TABLE  II 

ARAP06  MATH  MODEL  CAPABILITY 


Math  Model 


Degrees  of 
Freedom  Used 


Ml)  M2)  *  (3) 

L.  Isolated  Proprotor  (Rigid  Blades) 


a . 

Rigid,  constant  rpm 

0 

0 

0 

b. 

Rigid,  variable  rpm 

0 

l 

0 

c . 

Gimballed,  constant  rpm 

0 

2 

0 

d. 

Gimballed,  variable  rpm 

0 

3 

0 

e . 

Offset  flapping  hinge, 

constant  rpm 

0 

0 

2 

f . 

Offset  flapping  hinge, 

variable  rpm 

0 

1 

2 

g- 

Fully  articulated,  constant 

rpm 

0 

0 

4 

h. 

Fully  articulated,  variable 

rpm 

0 

1 

4 

Wing/Pylon  System  Coupled  with 

Proprotor  with  Rigid  Blades 

up 

a . 

Rigid,  constant  rpm 

to  4 

0 

0 

b. 

Rigid,  variable  rpm 

II 

1 

0 

c . 

Gimballed,  constant  rpm 

II 

2 

0 

d. 

Gimballed,  variable  rpm 

:i 

3 

0 

e  . 

Offset  flapping  hinge, 

constant  rpm 

1 1 

0 

2 

f . 

Offset  flapping  hinge, 

variable  rpm 

n 

1 

2 

g. 

Fully  articulated,  constant 

rpm 

n 

0 

4 

h. 

Fully  articulated,  variable 

rpm 

1 1 

1 

4 

3.  Wing/Pylon  System  Coupled  with 
Proprotor  with  Elastic  Blades 

up 

a.  Rigid,  constant  rpm  to  4 

b.  Rigid,  variable  rpm  " 

c.  Gimballed,  constant  rpm  " 

d.  Gimballed,  variable  rpm  11 

e.  Offset  flapping  hinge, 

constant  rpm  " 

f.  Offset  flapping  hinge, 

variable  rpm  " 

g.  Fully  articulated,  constant  rpm  " 

h.  Fully  articulated,  variable  rpm  " 

Ml)  Wing-Pylon  Degrees  of  Freedom 
*(2)  Proprotor  Rigid-Body  Degrees  of  Freedom 
M3)  Rigid  Blade  Elastic  Degrees  of  Freedom 
+  Rigid 

6  3 


0  4 

1  4 

2  4 

3  4 

0  2  +  2 

1  2  +  2 

I.  • 


l 


4 

5 

6 
7 

6 

7 

8 
9 


8 

9 

10 

11 

8 

9 

8 

9 


l_n  -p-  CO  n>  Co  fo  i — 1  O 


propeller)  the  rigid  rotor  option  is  used  and  blade  normal 
modes  employed  to  represent  blade  flexibility. 

To  represent  the  articulated  and  offset  flapping  hinge  rotor 
types  the  rigid  rotor  option  is  used  in  conjunction  with 
normal  modes.  The  first  normal  mode  is  used  to  represent 
the  blade  rigid-body  flapping.  Accordingly  the  natural 
frequency  of  the  first  mode  is 


wf  =  n  \/l  +  ^  —  (153) 

r  2  R 

where  ef  is  the  flapping  hinge  offset. 

The  out-of-plane  mode  shape  is  zero  at  the  hinge  and  inboard 
of  the  hinge,  and  is  unity  at  the  tip,  A  straight  line  is 
assumed  between  the  hinge  and  the  tip.  The  inplane  shape  is 
zero,  root  to  tip. 

For  the  fully  articulated  rotor  the  lead-lag  motion  is 
simulated  using  the  second  normal  mode.  The  frequency  is 
given  by 


“i  -  "Vfx  <154) 

where  eL  is  the  lead-lag  hinge  offset. 

The  out-of-plane  mode  shape  is  zero  and  the  inplane  mode 
shape  is  zero  at  the  lag  hinge,  zero  inboard  of  the  lag 
hinge  and  unity  at  the  tip.  A  straight  line  approximation 
is  used  between  the  lag  hinge  and  the  tip. 

Since  four  normal  modes  can  be  specified,  the  first  and 
second  flapwise  and  first  edgewise  modes  can  also  be  input, 
for  offset  flapping  hinge  rotors  and  the  first  flapwise  and 
edgewise  modes  in  the  case  an  articulated  rotor  is  being 
simulated. 

(b)  Wing-Pylon  Modeling  -  Up  to  four  wing-pylon  modes  can  be 
specified.  These  can  be  used  to  represent  a  complete  wing- 
pylon-hub  system  or  if  desired  pylon  pitching  and/or  yawing 
freedom  only. 

When  the  wing-pylon-hub  system  is  modeled  either  cantilever 
wing  or  free-free  (symmetric  or  asymmetric)  modes  may  be  used. 
Usually  the  fundamental  beam,  chord,  and  torsion  modes,  and 
either  the  fuselage  vertical  bending  or  a  second  wing  bending 
mode  should  be  used.  These  modes  must  be  calculated  at  the 
pylon  conversion  angle  input  to  the  program, 
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The  pylon  pitching  and/or  yawing  modes  may  be  modeled  without 
having  to  first  perform  a  normal  mode  vibration  analysis. 

First  the  wing  mode  shape  is  set  to  zero.  Then  the  conversion 
axis  reference  point  is  assumed  to  be  located  at  the  pylon 
center  of  gravity.  (The  mast  length  is  then  input  as  the 
distance  between  the  pylon  center  of  gravity  and  the  hub.) 

The  pylon  pitch  and  yaw  angular  deflection  is  set  equal  to 
the  distance  from  the  pylon  center  of  gravity  to  the  hub, 
divided  by  the  actual  mast  length.  The  desired  pylon  natural 
frequencies  are  input  directly. 

The  inherent  capability  to  build  up  the  math  model  from  a 
limited  number  of  degrees  of  freedom  to  a  fully-coupled  wing- 
pylon-hub-proprotor  system  can  be  used  when  first  using  the 
analysis.  At  each  stage  of  build-up  the  math  model  should 
be  checked  to  verify  the  accuracy  of  the  representation. 

For  example,  the  in-vacuo  coupled  wing-pylon-hub  frequencies 
and  mode  shapes  can  be  determined  by  running  the  program 
with  the  air  density  and  structural  damping  set  to  zero. 

B.  BLADE  FOLDING  DYNAMICS  ANALYSIS 

The  development  of  a  dynamic  theory  for  blade  folding/unfolding 
is  divided  into  two  parts:  a  flutter  (and  divergence)  analysis 
and  a  dynamic  response  (loads)  analysis.  The  flutter  analysis 
is  the  basis  for  the  response  analysis. 

1.  Flutter  Analysis 

The  flutter  analysis  is  based  on  a  conventional  collocation 
method,  flutter  and  vibration  analysis.®  The  formulation 
utilizes  structural,  aerodynamic,  and  inertial  characteristics 
in  the  form  of  structural  and  aerodynamic  influence  coefficients 
and  mass  matrices. 

This  analysis  differs  from  the  original  analysis, 8  in  two 
ways:  (1)  a  specific  structural  influence  coefficient  matrix 

which  represents  the  folding-proprotor  wing-pylon-blade 
system  has  been  developed,  and  (2)  a  faster,  more  accurate, 
eigenvalue  routine  is  employed. 

a.  Description  of  the  Math  Model 

An  overview  of  the  math  model  is  shown  in  Figure  19.  This 
model  includes  the  following  degrees  of  freedom: 

1.  Wing  beam,  chord,  and  torsion  (5  elements) . 

2 .  Pylon  pitch  and  yaw  relative  to  the  wing-tip 

3.  Motions  about  the  fold  hinge  and  feathering  axis 

4.  Blade  beam,  chord,  and  torsion  (5  elements) 
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The  standard  sign  convention  is  used,  positive  starboard  down, 
and  aft.  Nose  up,  yaw  right,  and  roll  left  are  positive. 

The  wing  representation  consists  of  five  beam  elements 
located  along  the  elastic  axis  of  the  wing.  Lumped  masses 
represented  by  dumbells  are  located  at  each  node  point 
(element  junction).  These  reflect  the  mass,  inertia,  and 
static  unbalance  properties  of  the  wing.  The  sweep  of  the 
x  axis  orientation  follows  from  the  right  hand  rule.  The 
principal  benefit  of  this  coordinate  system  is  that  the 
aerodynamic  and  mass  matrices  can  be  formed  directly  without 
coordinate  transformation  (except  to  account  for  the  dumbbell 
representation) . 

b.  Formulation  of  Equations  of  Motion 

The  equation  of  motion  for  the  system  shown  in  Figure  19  may 
be  written  as 


[  m ] {q}  +  [1  +  ig]  [K]  {q}  =  [  F  ]  {q}  (155) 


where 


m 

K 

F 

g 


{q> 


is  a  mass  matrix 

is  a  stiffness  matrix 

is  a  complex  forcing  function 

is  the  complex  structural  damping  parameter 

is  the  generalized  coordinates 


The  generalized  coordinates  used  in  this  analysis  are  those 
compatible  with  the  aerodynamic  representation.  The  vertical 
(z)  axis  is  taken  normal  to  a  strip  and  the  y  axis  in  the 
streamwise  direction;  the  wing  elastic  axis  is  included  in  the 
structural  formulation.  (The  segment  of  the  elastic  axis 
continuing  from  the  fifth  node  point  is  considered  to  be 
identical  to  the  fifth  segment  and  of  one-half  the  fifth 
segment ' s  length . ) 

The  wing  aerodynamic  representation  consists  of  strips 
centered  on  the  node  points.  The  chord  and  width  of  each 
strip  are  input  variables.  The  strips  lie  parallel  to  the 
freestream.  Two  dimensional  incompressible,  oscillating,  thin 
airfoil  theory  is  used  to  approximate  the  airloads. 


The  pylon  (shaft)  is  attached  to  the  wing  by  an  inelastic 
chordwise  element  between  the  wing  elastic  axis  and  the 
pylon  conversion  axis  and  by  a  spanwise  element  to  the  shaft 
centerline.  The  torsional  stiffness  of  the  spanwise  element 
is  used  to  simulate  the  pylon  pitch  attachment  stiffness;  and 
its  chordwise  stiffness  the  pylon  yawing  stiffness. 
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TWO  AND  THREE  HAVE 


Figure  19.  Overview  of  Math  Model  for  Blade  Folding  Flutter  Analysis 


Two  strips,  located  on  the  pylon,  may  be  used  to  approximate 
the  vertical  and  lateral  aerodynamic  forces  acting  on  the 
pylon. 

The  hub  details  include  blade  feathering  (control  system) 
and  fold  hinge  flexibilities.  Feathering  flexibilities  are 
immediately  located  inboard  and  outboard  of  the  fold  hinge. 

The  fold  hinge  radial  location  is  specified  by  the  hub  offset. 

The  blades'  structural  and  aerodynamic  representation  is 
identical  to  that  of  the  wing  except  that  the  element  out¬ 
board  of  the  fifth  node  is  inelastic.  The  blade  is  assumed 
to  start  at  the  fold  hinge.  The  fold  angle  is  a  variable. 

Only  three  blades  are  considered  and  they  are  assumed  to 
fold  knife-like  into  the  nacelle. 

c .  Mass  and  Aerodynamic  Matrices 

The  mass  and  aerodynamic  formulation  of  the  original  flutter 
analysis  are  essentially  unchanged.  For  completeness, 
however,  a  brief  review  of  these  formulations  is  given  here. 

(1)  Mass  Matrix  -  Figure  20  shows  the  details  of  the  dumbbell 
arrangement.  The  dumbbells  are  parallel  to  the  y-axis 
(free  stream).  Working  backwards,  the  deflections  at  the 


element  elastic  axis  may  be  written  as 

6  =  hx  +  X£ 

(156) 

and 

h2  -  hl 

5  d 

(157) 

where 


h^  and  are  the  generalized  coordinates. 

Solving  in  terms  of  h^  and  h2  only, 


and 


hjY  +  Xh2 
=  d 


Y  =  d  -  X 

or  in  a  matrix  notation 


(158) 

(159) 
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(160) 


*  *2 


From  Lagrange's  equation  the  kinetic  energy  of  a  system  is 
given  by 


T  =  h 


m  S  |6 


S  Ie\  K 


(161) 


where 


is  the  mass 

is  the  unbalance  (positive  forward) 
is  the  inertia  about  the  elastic  axis 


Performing  the  differential  ,  gives  the  following  form 

for  the  mass  matrix  in  dumbbell  coordinates. 


4ir-  siri  £i  >  s 

=  [u]  [u  ]x 

is  T  _1  1  1  S  I 

dJL  dl”d  ® 


(162) 


For  the  chordwise  deflections  no  transformation  is  required. 

(2)  Aerodynamic  Matrix  -  The  matrix  of  aerodynamic  influence 
coefficients  for  a  strip  is  defined  in  terms  of  the 
generalized  coordinates  as 


in  which 


?  ■  Pw2br2slch]!  1 


(163) 


is  the  air  density 
is  the  oscillatory  frequency 
is  the  reference  semichord 
is  the  reference  semispan 


i 


The  matrix  of  aerodynamic  influence  coefficients  [C^]  are  in 
partitioned  form,  i,e. 


[chl 


(164) 


[Cw]  is  related  to  the  lift  and  moment  acting  on  the  strip  as 
follows.  Using  the  geometry  indicated  in  Figure  20b,  the 
relationship  between  the  lift  and  moment,  L  and  M,  respectively, 
may  be  written  as 


1  11  F 


0  dJ(F2 


(165) 


The  oscillatory  lift  and  moment  referred  to  the  blade  quarter 
chord  are  defined  as 


Trco*A(pwab2AY) 


1  0 1  [Lh  LJp 
0  bJlMu  mJIi 


(166) 


when  the  coefficients  Lh,  Lfl,  M^,  and  Ma  are  those  given  in 
Reference  10,  and  h  and  5  a re  related  to  h]^  and  h2  as  follows: 


1  Olfh 


i *  n-i  yu 

Solving  for  {F }  in  terms  of  {h}  an  expression  for  [Ch]is 
obtained 


(167) 


[Ch]  -  7rcosA(^- 


-sip*  L«  1 

_0  Mh  Ma  j 


(168) 


The  steady-state  aerodynamic  influence  coefficients  are  a 
limiting  case  of  the  oscillatory  coefficients,  i.e.  w+0. 


;;j  - 


(169) 


d.  Flexibility  Matrix 


The  flexibility  matrix  for  the  wing-pylon -blades  system  is 
formed  using  finite  element  theory.1-*  As  noted  earlier,  the 
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a _ 


structural  representation  consists  of  30  beam  elements  and  31 
nodes.  Figure  21  shows  the  structure  arrangement.  Each  beam 
element  is  considered  to  have  a  uniform  cross  section  and 
resist  moments  and  shears  about  the  principal  axes,  torsion 
about  the  normal  axis,  and  axial  loads.  Each  element  stiff¬ 
ness  matrix  is  12  x  12  and  is  of  the  form  shown  in  Figure  22. 
Each  node  is  assumed  to  have  six  degrees  of  freedom. 


To  construct  a  stiffness  matrix  for  the  complete  structure 
each  element  is  first  oriented  along  the  reference  x-axis 
(see  Figure  21)  and  then  rotated  into  its  structural  orienta¬ 
tion  through  three  successive  coordinate  system  rotations. 
These  rotations  are  in  order,  first  about  the  z-axis,  then 
about  the  y'-axis,  and  finally  about  the  z''-axis.  These 
rotations  are  designated  as  a,  \p,  and  B>  respectively,  and 
the  rotated  coordinate  systems  are  disgnated  as  x^yiz1, 
xIIyIIzII,  an(j  xIIIyIIIzIII , 


These  coordinate  transformations  are  given  by  the  following 
relationships  (for  brevity  only  the  xyz  transformations  are 
shown.  The  0X,  0y,  and  8Z  transformations  are  similar.) 


a  rotation 

I 


4>  rotation 


II 


B  rotation 

III 


(170) 


(171) 


(172) 


The  final  transformation  from  the  reference  coordinate  system 
to  the  final  structural  coordinate  system  is  then 
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Figure  22.  Element  Stiffness  Matrix 
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or 


For  reference  the  rotations  of  the  various  structural  elements 


shown  in  the  math 

model  (Figure  21) 

are 

Element (s) 

Segment  No. 

a 

J L 

B 

Wing 
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0 

0 

Wing  sweep 
angle 

EA  -  CA  Offset 
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Mast 
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60 
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0 

Blade  fold 

Blade  No .  2 

18  -  23 

180 

60 

-  Blade  fold 

Blade  No .  3 

25  -  30 

180 

-60 

-  Blade  fold 

The  stiffness  matrix  for  each  element,  in  the  beam  coordinate 
system  is  then 


fcj  =  [T]T[ke][T] 
2  2 


(175) 


where  [ke]  is  the  element  stiffness  matrix  shown  in  Figure  22. 


The  stiffness  matrix  for  the  complete  structural  system  is 
then  formed  by  adding  the  individual  stiffness  matrices  at 
each  node.  This  gives  a  186  x  186  matrix 


where  the  overlay  of  the  individual  element  stiffness  matrices 
indicate  addition  of  adjoining  element  stiffness  matrices. 

This  stiffness  matrix  (176)  must  now  be  transformed  into  a 
coordinate  system  compatible  with  the  inertial  and  aerodynamic 
math  models.  As  noted  earlier  these  are  normal  to  the  aero¬ 
dynamic  surface.  For  the  wing  and  blade  l(see  Figure  21)  the 
reference  xyz  coordinate  system  is  compatible  with  the 
inertial  and  aerodynamic  coordinate  systems.  However,  for 
blades  2  and  3  a  coordinate  transformation  to  the  xivyivziv 
system  is  required.  The  transformation  is  about  the  y1  axis 
and  is  (only  the  linear  transformation  is  shown  for  brevity) : 


or 


IV 


‘cos  8 

0 

-sing 

(*) 

I* 

- 

0 

1 

0 

M 

\  *  * 

sing 

0 

cosg  _ 

(z) 

(177) 


IV 
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For  a  three-bladed  rotor  6  =  -60°  for  the  second  blade  and 
+60°  for  the  third  blade . 

The  final  stiffness  matrix  is  obtained  by  pre  and  post  multi¬ 
plication  of  [K]by  [u  ] 3 .  This  gives 

IK)final  =  [u]  [K][U]3  (179) 

3 


and  is  a  186  x  186  symmetric  matrix. 


The  stiffness  matrix  is  now  condensed  to  the  number  of  degrees 
of  freedom  required  for  this  study.  This  is  accomplished  by 
transferring  the  contribution  of  the  undesired  coordinates  to 
the  retained  coordinates.  Single  coordinate  elimination  is 
employed  such  that 


reduced 


'  l*ul 


final 


l 


(180) 


Where  i  and  j  refer  to  columns  and  rows  and  n  is  the  coordinate 
of  the  row  and  column  being  eliminated.  Bv  successive 
application  of  180  the  stiffness  matrix  is  reduced  to  74  x  74. 
Figure  23  shows  the  coordinates  retained.  The  reader  can  refer 
to  Figure  23  for  identification  of  the  degrees  of  freedom 
which  are  retained  in  the  reduced  stiffness  matrix. 


The  reduced  stiffness  matrix  is  now  inverted  to  form  the 
flexibility  matrix.  A  final  transformation  into  the  colloca¬ 
tion  (dumbbell)  coordinates  is  now  made. 

By  definition  the  system  deflections  are  given  as 

{6}  -  [SIC]  {F}  (181) 

where 

{6}  is  the  deflection  matrix,  consisting  of 
linear  and  angular  deflections 

[SIC]  is  the  flexibility  matrix 

{F}  is  the  applied  force  matrix 

The  relationship  between  the  two  linear  deflections  at  the 
collocations  points  and  the  linear  and  angular  deflections  at 
the  element  elastic  axis,  derived  for  the  mass  matrix  also 
apply  for  the  stiffness  matrix.  Thus 


77 


1 

X 

49 

X 

97 

X 

* 

31 

145 

X 

2 

Y 

50 

Y 

* 

11  98 

Y 

* 

32 

146 

Y 

©  4 

Z 

eX 

© 

51 

52 

Z 

®x 

* 

* 

12(f *  99 
13^100 

Z 

eX 

* 

* 

33 

34 

© 

147 

148 

Z 

0x 

5 

9  Y 

53 

0Y 

101 

0y 

149 

6  Y 

6 

ez 

54 

6  z 

* 

14  102 

ez 

150 

9  z 

7 

X 

55 

xz 

103 

xz 

151 

xz 

8 

Y 

56 

Y 

104 

Y 

152 

Y 

©  9 
10 

Z 

ex 

* 

« 

1 

2 

© 

57 

58 

z 

ex 

©105 

^106 

Z 

6X 

@ 

153 

154 

Z 

6X 

11 

0y 

59 

0  Y 

107 

0  Y 

155 

®  Y 

12 

ez 

60 

eZ 

108 

®Z 

156 

®Z 

13 

X 

61 

x 

109 

X 

157 

X 

* 

55 

14 

Y 

62 

Y 

110 

Y 

158 

Y 

* 

56 

®  15 
W  16 

Z 

ex 

* 

« 

3 

4 

© 

63 

64 

Z 

ex 

©HI 

w112 

Z 

ex 

@ 

159 

160 

Z 

ex 

* 

* 

57 

58 

17 

6y 

65 

0  Y 

113 

e  y 

161 

®  Y 

18 

0Z 

66 

6  Z 

114 

ez 

162 

ez 

19 

X 

67 

X 

115 

X 

* 

35 

163 

x 

* 

59 

20 

Y 

68 

Y 

116 

Y 

* 

36 

164 

Y 

* 

60 

_  21 

Z 

* 

5 

69 

Z 

117 

Z 

* 

37 

165 

Z 

* 

61 

@22 

eX 

* 

6 

© 

70 

eX 

©118 

9X 

* 

38 

© 

166 

0X 

* 

62 

23 

eY 

71 

0  Y 

119 

®  Y 

167 

6  Y 

24 

eZ 

72 

ez 

120 

6Z 

168 

0Z 

25 

X 

73 

X 

★ 

15  121 

X 

* 

39 

169 

X 

* 

63 

26 

Y 

74 

Y 

* 

16  122 

Y 

* 

40 

170 

Y 

* 

64 

@  27 
V2'  28 

Z 

* 

7 

© 

75 

Z 

* 

!W23 

18^124 

Z 

* 

41 

© 

171 

Z 

* 

65 

ex 

* 

8 

76 

6X 

* 

®x 

* 

42 

172 

ex 

* 

66 

29 

0Y 

77 

6  Y 

125 

6  Y 

173 

0Y 

30 

ez 

78 

ez 

126 

ez 

174 

9 

31 

xz 

79 

xz 

* 

19  127 

xz 

* 

43 

175 

xz 

* 

67 

32 

Y 

80 

Y 

* 

20  128 

Y 

* 

44 

176 

Y 

* 

68 

33 

Z 

« 

9 

81 

Z 

* 

21  129 

Z 

* 

45 

177 

Z 

* 

69 

®34 
^  35 

ex 

0Y 

* 

10 

© 

82 

83 

ex 

0  Y 

* 

2W30 

^131 

ex 

0  Y 

* 

46 

© 

178 

179 

®X 
0  Y 

* 

70 

36 

®z 

84 

0  z 

132 

6  z 

180 

®Z 

37 

X 

85 

* 

23  133 

xz 

* 

47 

181 

X 

* 

71 

38 

Y 

86 

Y 

* 

24  134 

Y 

* 

48 

182 

Y 

* 

72 

©  40 

Z 

0X 

© 

87 

88 

Z 

ex 

* 

* 

25— .135 
26©136 

Z 

ex 

* 

* 

49 

50 

© 

183 

184 

Z 

ex 

* 

* 

73 

74 

41 

0Y 

89 

0  Y 

137 

ey 

185 

0  y 

42 

0Z 

90 

0Z 

138 

6Z 

186 

0Z 

43 

X 

91 

X 

* 

27  139 

X 

* 

51 

44 

Y 

92 

Y 

* 

28  140 

Y 

* 

52 

®tl 

Z 

ex 

© 

93 

94 

Z 

ex 

* 

* 

2W41 

30^142 

Z 

0X 

* 

* 

53 

54 

47 

eY 

95 

0  Y 

143 

0  y 

48 

ez 

96 

ez 

144 

eZ 

♦Indicates  coordinates  retained 


Figure  23.  Coordinates  Retained  in  Stiffness  Matrix. 
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Hence,  the  final  form  of  the  flexibility  matrix  is 
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(e)  Method  of  Solution  -  Solution  of  the  matrix  of  differential 
equations  follows  the  usual  flutter  practice.  Equation  155 
is  rearranged  such  that 


or 


r  1  j 

SIC 


|  m  +  f]  J  -  (1.  *2i2~"  =  ° 


{q}  =  0 


(187) 


(188) 
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where 


X  =  iA  *  ^9)  =  XR  +  Xj  (Real  and  imaginary  parts) 

to^ 

[m]  «[SIC]([m]  +  [F]) 

The  complex  eigenvalue  X  is  solved  for  using  a  qr  algorithm 
called  ALLMAT  which  uses  Given's  rotations  for  the  elementary 
unitary  transformations.  The  Weilandt  inverse  power  method 
is  used  to  compute  the  complex  vectors . 

The  frequency  and  damping  is  then  obtained  from  the  following 
relationships . 

As  assumed, 


therefore 


and 


(189) 


(190) 

(191) 


(192) 


As  is  the  usual  flutter  practice,  a  wide  range  of  reduced 
frequencies  must  be  swept  to  construct  a  V-g  diagram. 

2 .  Response  Analysis 

The  Folding  Dynamics  Response  Analysis  solves  for  the  internal 
forces,  moments,  and  torques  acting  at  a  specified  node  on 
the  wing-pylon-blade  system  due  to  externally  applied  airloads. 
These  externally  applied  airloads  may  be  due  to  an  atmospheric 
gust  or  a  maneuver. 

It  must  be  pointed  out  that  the  response  analysis  determines 
only  the  incremental  loads.  The  steady  structural  loading 
must  be  added  to  these  incremental  loads  to  arrive  at  the  net 
load.  The  Feather-Stop  Analysis  (ARAP06)  may  be  used  to 
determine  the  steady  loads  for  a  trimmed  flight  condition. 

The  Feather-Stop  Analysis  also  has  the  capability  to  compute 
the  response  during  blade  folding  and  is  recommended  for  pre¬ 
liminary  studies  where  unsteady  aerodynamic  effects  can  be 
neglected.  80 


a.  Formulation  of  the  Response  Analysis 


The  response  analysis  employs  the  Laplace  transform-transfer 
function  method  to  determine  the  response  of  the  generalized 
coordinates  to  external  forces.  These  responses  are  then 
used  to  calculate  the  transient  internal  loads.  The  general 
method  may  be  explained  as  follows : 


The  dynamics  of  the  wing-pylon-proprotor  system  with  respect 
to  a  trim  or  static  position  can  be  represented  by  differential 
equations : 


2  dq 

[m]£_2n+  |d]  Jjli  +  [K]qn  =  [  F  (t)  j  +  [Fe(w2)qn 
dt 


(193) 


where 


[  m] 
I  D] 
IK] 

Fe 


is  the  mass  matrix 
is  the  damping  matrix 
is  the  stiffness  matrix 

is  the  external  aerodynamic  resulting  from 
wing  motion 


F  (t)  is  the  externally  applied  forces 


q  is  the  generalized  coordinates 

Assuming  the  classical  flutter  form  and  harmonic  motion,  then 
equation  193  is  rewritten  as 


(B  -  AI ) q  =  F' 


(194) 


where 


B 

I 

X 

F  ' 


_  l  _  l  _  l 

K  m  +  iD  IL  -  K  Fe 

w 


identity  matrix 
eigenvalue 


1  +  ig 
w 2 


K~  1 

“5T  F  (t) 


The  damping  matrix  [d]  will  be  assumed  equal  to  zero,  hence 

B  =  K-1m  -  K_1Fe 
81 


Note  that  if  F'  is  set  to  zero  equation  194-will  ha.  recognized 
as  the  classical  flutter  equation.  The  generalized  coordinate 
q^  can  be  solved  for  by  applying  Cramer's  rule 


(195) 


where  the  numerator  |(B  -  XI)  Ip'  is  the  determinant  of  the 
matrix  (B  -  XI)  with  the  i  column  replaced  by  F'. 

Since  only  a  finite  number  of  modes  will  define  the  system 
response,  equation  195  can  be  simplified  by  evaluating  it  at 
its  natural  frequencies  and  determining  the  residue.  The 
technique  can  be  explained  by  writing  equation  195  in  matrix 
form  and  evaluating  the  residue  at  the  first  pole  or  eigenvalue, 
i .  e . 


where 


Si  = 


|(B  -  bil )  |  F  t 
(b!  -  X)  j(bi  -  bx)| 


(196) 


is  the  first  eigenvalue 

b^  is  the  remainder  of  eigenvalues 

X  is  an  integration  operator 

The  expression  is  thus  evaluated  at  one  frequency  minus  the 
first.  Equation  196  can  now  be  written  as 


where 


(197) 


is  the  contribution  of  the  first  mode  to  the 
total  response  of  coordinate  q^ 

(198) 


|(B  -  bjl )  1  f 

PI  -  bx) 


The  total  response  of  q^  is  then  given  by  equation  199 

fP 

Dn  =  D1  +  °2  +  d3  ••• 

n=l 

where  NP  is  the  total  number  of  eigenvalues. 


(199) 


82 


However  A,  the  integration  operator,  is  difficult  to  evaluate 
as  a  function  of  time.  It  is,  therefore,  necessary  to  express 
199  as  a  transfer  function  to  determine  the  coordinate  transient 
response 


(200) 


This  can  be  expressed  in  terms  of  an  inverse  Laplace  transform 


q±(t)  =  F  (t)  L-VEill)  (201) 

q(Z) 

where  Z  is  the  differential  operator  . 

dt 

It  is  now  necessary  that  a  relationship  be  established  between 
the  operators  A  and  Z.  By  definition, 


A  =  1  i-2-?  and  Z  =  ia) 

to 

where 

w  is  the  frequency 

2a  is  the  damping  factor,  g 

hence 

1  a)2  J2  +  Z2ow 

A  *  1  +  i2a  *  "  1  +  4 a* 

Then  rewriting  equation  197 


then 


Since 


h. 

bxA 

m 

A  b^ 


1  w*2  u>i2  (  1  -  i2ai) 

hi  =  1  +  i2ai  =  TTTop 


f  Z2  +  Z2aw 
b^  '  1  +  4a2 
"Z2  +  "Z2  ow\  1 


(202) 


(203) 


(197) 


(204) 


(205) 
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We  have  the  form 


K1  / 

Ef  (z2  +  2aw1z) 

Z2  +  2au>]V  +  wz  (1  -  j2a) 

We  neglect  the  term  j2au2,  since  it  represents  an  imaginary 
response  and  steady-state  component  of  no  interest.  Thus, 

K1 

r  z  (2  +  2cio)i) 

Rl  =  bl _  (; 

(Z  +  aw) 2  +  w2  (1  -  a2 ) 

Equation  207  is  of  the  form  P ( z)  . 

q(z) 

To  inverse  transform  equation  207,  s  in  the  numerator  is 
neglected  and  differentiated  with  respect  to  time.  Only  the 
real  part  of  the  inverse  transform  is  used: 

K, 

Re  (H  w2 


(206) 


(207) 


D  =  —  *L 

1  (S  +  oiw)  2  +  w2  (1  -  a2) 


(208) 


This  manipulation  is  made  for  all  modes  of  importance 
Equation  189  can  be  written  as 


q(t)  N0 

Fit!  “  Dq  +  DjS  +  7? 


(209) 


where  N0,  DQ,  and  are  linear  coefficients  from  equation  197. 

The  expression  for  the  coordinate  motion  in  finite  differencing 
notation  is  then 

At2(NQF(t)  +  E.i)  +  4 (E_x  +  E_1At)  +  RiAt(2E_1  +  E.1At) 


q  (t)  = 


where 


D0At 2  +  2DiAt  +  4 


(210) 


is  the  time  increment 
is  A.  F  (t) 


For  finite  differencing 


MF<t>  - 

indicates  the  previous  time  step 
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The  advantages  of  this  method  are: 


1.  The  steady-state  gain  value  is  immediately  known  as 


a!  «  e4°) 

F  qToT 

2.  The  response  can  be  obtained  very  rapidly  by 
evaluating  only  a  few  roots  rather  than  a  large 
system. 


(211) 


3.  Large  multiroot  systems  can  safely  be  reduced  by 
making  simple  approximations. 

4.  The  magnitude  of  the  highest  retained  frequency  can 
be  used  to  define  an  optimum  value  of  At  used  for 
the  time  response . 

5.  The  solution  is  valid  for  any  time  form  of  the 
forcing  function  without  need  to  resolve. 

6.  The  solution  form  is  flexible  and  can  be  changed 
with  a  minimum  of  difficulty. 


Approximations  made  during  solution  to  reduce  the  root  system 
are: 

1.  If  roots  of  p(S)  are  equal  to  roots  of  q(S)  within 
one  percent,  their  roots  effectively  cancel  from  the 
transfer  function  and  are  thus  removed. 


2.  If  natural  frequencies  of  the  system  (roots  of  q(S)) 
are  greater  than  some  value,  usually  200  HZ, 
they  are  considered  sufficiently  far  removed  from  the 
forcing  frequencies  to  have  little  effect  on  response. 
These  roots  are  neglected. 

Special  Notes: 

1.  The  evaluation  of  the  determinant  of  a  complex  74th 
order  matrix  is  a  time  consuming  operation.  Thus 
the  matrix  is  reduced  to  a  6th  order  for  each  input 
made  by  matrix  pivotal  reduction  techniques.  The 
reduced  6x6  matrix  then  provides  an  equivalent 
system  as  far  as  evaluation  by  Cramer's  rule  is 
concerned.  The  forcing  vector  is  also  condensed  to 
retain  system  eigenvalence . 
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2.  The  roots  of  the  unforced  system  are  available  directly 
from  the  flutter  analysis.  These  are,  however,  based 
on  the  assumption  of  simple  harmonic  motion,  i.e. 

the  damping  parameter  g  is  the  value  of  structural 
damping  required  to  maintain  neutral  stability. 

These  may  be  modified  for  use  in  the  response 
analysis  using  the  method  developed  by  Frueh  and 
Millerl4 

“response  °  “flutter 
sresponse  =  **  [l1  (^)(|f)]  Aflutter 

where 

V  =  velocity 

to  =  model  natural  frequency 

^  »  rate  of  change  of  model  natural 

frequency  with  velocity  evaluated  at 
the  airspeed  V 

This  transformation  must  be  made  manually.  However, 
it  is  seldom  required  for  FPR  studies  since  the  rate 
of  change  of  frequency  is  small  in  the  transition 
mode . 

3.  The  external  forcing  function  F(t)  is  written  as  a 
unit  vector  and  the  superposition  theorem  employed 
to  facilitate  the  representation  of  complicated 
forcing  functions.  The  vector  is  calculated 
assuming  steady-state  aerodynamic  theory,  i.e. 

Cu  =  1.  The  gust  induced  lift  for  aerodynamic 

i 

section  is  then: 


(212) 

(213) 


L.  =  W  PCT  AY.  C. 
1  Q  Jj  1  1 

y  a 


(214) 


where 

Vg  *  gust  velocity 

AYj.  =  width  of  strip  i.  For  blades  2  and  3 
this  input  is  \  AY^  for  blade  1  to 
account  for  blade  orientation 

p  =  air  density 
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lift  curve  slope 


Cl  " 

a 

■  strip  semichord 
V  ■  forward  velocity 

4.  The  R.  T.  Jones  approximation  to  Wagner's  function10 

-0.0455T  -0 . 3t 

4(t)  -  1  -  0.165e  -  0 . 335e  (215) 

is  used  to  account  for  aerodynamic  lag  (t  is  the 
distance  traveled  in  semichords) . 

b.  Calculation  of  the  Internal  Loads 


In  order  to  determine  the  internal  loads  at  a  given  node  the 
linear  displacements,  qn,  at  the  dumbbell  coordinates  the 
node  and  at  the  node  to  its  left  are  computed.  These  are 
then  transformed  into  the  displacements  and  rotations  at  the 
node  and  node  to  its  left. 

The  torque  is  calculated  directly  from  the  differential  twist 
of  the  nodes. 


,216> 

The  moment  and  shear  are  approximated  by  solving  simultaneous 
equations  for  the  deflection  and  slope  at  the  tip  of  a 
cantilevered  beam. 


6  ■  +  ^2liZ  (217) 

and 

5  -  3e1Jli2  +  2e2Jli  (218) 

Note  that  the  5  expression  is  equal  to  36/3 i. 

The  displacement  at  the  node  and  node  lift  are  then  substituted 
into  equation  217  to  determine  the  constant  and  e2.  Once 
e^  and  e2  are  known  they  are  used  to  solve  for  slope,  9.  In 
general  form  the  slopes  at  the  nodes  are 
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(219) 


(•‘1 

- 

£12/2  -  i2a1  -a ! 

Vl’  -  V/«  Y/2" 
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k 

“^2  2/2  -12 

l23/3  H2!/2 

m  “ 

kl 

The  moment  and  shear  are  then  given  as 


M 


2EIg 


6l)+  26R  +  0L 


and 


(220) 


6Eiflr 


SII  = 


2 

F7 

i 


(6r  ■  6l)+  eR  +  eL 


(221) 


(This  approximation  is  exact  for  a  uniform  beam.  For  a 
nonuniform  beam,  with  nonuniform  loading  it  is  only  approximate. 
However,  numerous  checks  have  established  the  error  as  small.) 
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SECTION  III 


CORRELATION  OF  THEORY  WITH  EXPERIMENTAL  DATA 


Prior  to  this  program  a  considerable  amount  of  experimental 
data  on  proprotor  stability  and  dynamic  response  were  available. 
These  data  proved  valuable  during  the  development  and  checkout 
of  the  computer  programs.  A  wind  tunnel  test  of  a  dynamically 
scaled  model  was  made  during  this  program  to  obtain  additional 
data  on  proprotor  stability  at  low  rpm,  and  on  the  effect  of 
feathering  rate  on  dynamic  response.  A  joint  NASA-Bell 
Helicopter  Company  test,  using  the  same  model,  provided 
experimental  data  on  the  stability  and  load  characteristics 
during  blade  folding. 

Appendix  I  contains  descriptive  data  for  the  models  used  to 
obtain  the  experimental  data  shown  in  this  section.  These 
descriptive  data  are  sufficient  to  allow  the  reader  to  conduct 
independent  dynamic  analysis. 

A.  WING/BLADE  AERODYNAMIC  INTERFERENCE 

Data  on  wing/blade  aerodynamic  interference  were  available  from 
a  test  of  the  model  shown  in  Figure  24.  This  model  was  con¬ 
structed  specifically  to  investigate  wing/blade  aerodynamic 
interference.  Provisions  incorporated  in  the  model  to  vary 
the  proprotor  to  wing  spacing,  wing  chord  length,  wing  leading 
edge  shape ,  and  pylon  and  wing  angle  of  attack .  A  three- 
component  strain  gage  balance  was  employed  to  measure  the  hub 
normal  and  side  force,  pitching  and  yawing  moments,  and  axial 
force.  The  wind-tunnel  balance  was  used  to  measure  the  combined 
wing  and  proprotor  forces  and  moments.  A  strain  gage  was 
located  at  the  blade  root  fitting  to  measure  the  blade  root 
bending  moment  in  the  plane  of  the  rotor. 

Figure  25  shows  the  correlation  between  theory  (Program  ARAP06 ) 
and  the  measured  variation  in  the  blade  inplane  bending  moment 
as  the  blade  passes  in  front  of  the  wing.  The  experimental 
data  were  obtained  by  indexing  the  blade  in  five  degree  incre¬ 
ments  of  azimuth  and  recording  the  bending  moment.  (Zero 
azimuth  is  defined  to  be  the  designated  master  blade  normal  to 
the  upper  surface  of  the  wing.)  The  data  have  been  normalized 
on  dynamic  pressure  since  the  moments  were  found  to  vary 
linearly  with  dynamic  pressure. 

The  oscillatory  normal  force  correlation  is  shown  in  Figures 
26  and  27.  The  relative  influences  of  wing-rotor  spacing, 
shaft  angle  of  attack,  and  wing  angle  of  attack  on  the 
oscillatory  normal  force  are  indicated. 
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The  normal  force  correlation  with  regard  to  waveform  is  shown 
in  Figure  28.  Although  there  appears  to  be  a  sizeable  varia¬ 
tion  between  the  calculated  and  measured  waveform,  the  harmonic 
content  at  the  blade  passage  frequency  is  nearly  equal.  Table 
III  compares  the  frequency  content  in  terms  of  the  blade 
passage  frequency.  The  phase  angle  is  also  shown.  It  is 
evident  that  the  principal  difference  lies  with  the  magnitudes 
and  phase  angles  of  the  higher  harmonics  of  the  blade 
passage  frequency. 


TABLE  III 


Calculated 
in2 ,  60  degrees 
in2 ,  40  degrees 
in2 ,  5  degrees 


lift  and 

pitching  moment  due  to  wing-blade  aerodynamic  interference. 

The  decision  to  neglect  the  effect  of  the  blades  on  the  wing 
in  the  Feather-Stop  Analysis  was  based  on  this  experimental 
evidence . 

B.  DYNAMIC  STABILITY  CORRELATION 

The  model  shown  in  Figure  29  was  used  to  obtain  dynamic 
stability  data  for  the  proprotor  over  a  range  of  rpm  from  zero 
to  a  scale  tip  speed  of  600  feet  per  second,  and  with  the 
proprotor  stopped,  over  the  full  range  of  blade  fold  angles. 

It  should  be  noted  that  the  model  configuration  employed  for 
the  fold  test  differed  somewhat  from  that  of  the  rotating 
test.  This  was  the  result  of  having  to  insert  a  fold-hinge 
fitting  between  the  blade  and  the  hub  for  the  blade  folding 
test.  This  fitting  could  not  be  used  for  the  rotating  tests 
since  it  was  not  designed  to  carry  centrifugal  force. 

The  model  is  not  scaled  to  represent  a  specific  FPR  design 
but  is  representative  of  the  dynamic  and  aeroelastic  character¬ 
istics  of  the  FPR  concept.  The  basis  for  the  scaling  is 
discussed  in  Appendix  I.  The  wing  and  blades  are  dynamically 
scaled  in  beam,  chord,  and  torsion.  The  pylon  is  not  elastically 
scaled  but  does  have  representative  inertial  properties.  The 
proprotor  is  of  the  gimballed,  semirigid  type,  i.e.  the  blades 
are  rigidly  attached  to  a  central  yoke  which  is  in  turn 


NORMAL 

Harmonic 

1  BP  (3  per  rev) 

2  BP  (6  per  rev) 

3  BP  (9  per  rev) 

Measured  data  show  very 


FORCE  HARMONIC  CONTENT 
Measured 

26  in2,  60  degrees  25 

9.5  in2,  60  degrees  14 
15  in2 ,  10  degrees  10 

little  variation  in  wing 
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attached  to  the  hub  with  a  gimbal.  The  hub  may  be  elastically 
restrained  and  the  blade  pitch-flap  coupling  varied.  The 
collective  pitch  of  the  blades  can  be  remotely  adjusted  to 
control  rpm. 

For  the  folding  tests  a  fold  hinge  fitting  is  added  between 
the  blades  and  the  yoke  and  the  gimbal  is  locked  out.  Con¬ 
trolled  flexibility  is  included  in  the  fitting  to  simulate 
fold  hinge  and  blade  pitch  control  flexibility. 

1.  Proprotor/Pylon  Stability 

Proprotor  stability  was  determined  for  two  model  configura¬ 
tions.  In  the  first  configuration,  the  proprotor  gimbal  was 
restrained  by  a  hub  restraint  spring  with  a  rate  typical  of 
that  employed  in  Bell's  proprotor  designs.  This  spring  rate 
was  found  to  be  a  significant  factor  in  the  stability  at  low 
rpm  and  was  varied  over  a  range  to  determine  its  influence. 

The  blade  pitch-flap  coupling  was  also  varied  to  establish  its 
influence . 

The  model  was  also  tested  with  the  flapping  locked  out.  How¬ 
ever,  this  configuration  was  not  typical  of  current  helicopter 
hingeless  rotors  as  the  blade  first  elastic  flapping  natural 
frequency  was  approximately  1.8  per  rev,  and  these  rotors 
generally  have  natural  frequency  ratios  of  1.15  to  1.2  per  rev. 
The  configuration  was,  of  course,  typical  of  a  gimbal led  rotor 
configuration,  where  the  gimbal  freedom  is  locked  out  prior 
to  blade  feathering  and  folding. 

The  experimental  stability  data  were  obtained  by  maintaining 
constant  rpm  and  incrementally  increasing  the  wind-tunnel 
velocity  until  the  model  motion  was  lightly  damped.  At  each 
velocity  the  model  was  manually  excited  using  thin  cables 
attached  to  the  wing  tip.  The  decay  of  the  wing  and  proprotor 
motion  was  recorded  and  used  to  obtain  frequency  and  damping 
information.  Figure  30  shows  a  typical  decay  history.  Figure 
31  a  typical  variation  in  wing  frequency  and  damping  with 
velocity.  Similar  plots  were  obtained  for  the  wing  chord  mode. 

The  calculated  stability  data  were  obtained  in  a  similar  manner. 
The  computer  program  was  started  with  the  wing  undeflected 
so  that  the  simulated  gravitational  force  and  the  airloads 
would  excite  the  system  natural  modes.  Frequency  and  damping 
were  obtained  from  the  time  history  of  the  model  motions. 

a.  Proprotor  Stability  with  Gimbal  Free 

Figure  32  shows  the  stability  correlation  with  the  model 
gimbal  free.  For  reference  the  proprotor  rpm  at  which  wing- 
pylon  modes  are  in  resonance  with  proprotor  excitation  are 
indicated. 
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Where  possible,  the  tunnel  velocity  was  increased  until  the 
model  was  neutrally  stable  or  until  the  neutral  stability 
velocity  was  evident  from  extrapolation  of  a  damping  versus 
velocity  plot.  However,  in  many  cases,  the  instability 
velocity  could  not  be  reached  because  of  excessive  vibration 
and/or  loads.  In  these  cases  the  maximum  airspeed  tested  is 
indicated  by  an  arrow  at  the  data  point. 

In  the  range  from  250  to  350  rpm,  instability  occurs  in  the 
wing  beamwise  bending  mode.  In  this  instability  simultaneous 
wing  beamwise  bending  and  torsion  are  the  dominant  model 
motions.  There  is  very  little  blade  flapping  with  respect  to 
space  since  the  rotor  cannot  follow  pylon  pitching  or  yawing 
motions  at  the  wing  beam  natural  frequency. 

The  wing  beamwise  instability  is  well  understood  and  readily 
predicted  as  evidenced  by  the  correlation.  The  mechanism  of 
this  instability  is  discussed  in  References  4,  5,  and  6. 

From  50  to  250  rpm  wing  torsion  mode  instability  is  predicted. 
This  instability  is  similar  to  the  wing  beam  mode  instability 
in  that  it  involves  blade  rigid-body  flapping  coupled  with 
elastic  wing  motion.  However,  the  frequency  of  instability  is 
much  higher;  occurring  at  the  wing  torsion  mode  frequency, 
about  2*5  times  that  of  the  wing  beam  mode. 

At  very  low  rpm  and  near  zero  rpm  the  blade  rigid-body  flapping 
motion  becomes  unstable.  The  frequency  of  this  instability  is 
close  to  the  uncoupled  frequency  of  the  proprotor  flapping  on 
the  hub  flapping  restraint.  Wing  vertical  bending  and  orsion 
are  in-phase  with  the  flapping.  At  zero  rpm  the  instability 
appears  as  a  pitching  of  the  proprotor  tip-path-plane;  at  low 
rpm  it  appears  as  a  forward  precession  of  the  tip-path-plane. 

A  time  history  of  this  flapping  instability  is  shown  in  Figure  33. 

The  mechanism  of  the  low  rpm  instability  was  not  understood  at 
the  time  of  the  test.  To  gain  insight  into  the  instability 
mechanism  and  to  provide  additional  correlation  data  the  model 
pitch-flap  coupling  and  hub  restraint  were  varied. 

Figures  34  and  35  show  the  variation  in  stability  with  these 
parameters.  The  principal  influence  appears  to  be  the  flapping 
restraint,  with  pitch-flap  coupling  also  being  destabilizing. 

At  the  point  of  minimum  stability,  KHUB  ■  2500  foot-pound/degree, 
the  uncoupled  flapping  natural  frequency  of  the  rotor  on  the 
hub  restraint  spring  is  equal  to  the  uncoupled  wing  beamwise 
natural  frequency.  Thus,  the  instability  mechanism  appears 
to  be  the  frequency ..  coincidence  of  the  flapping  and  wing  beam- 
wise  bending  modes.  Correlation  with  the  low  rpm  data  is  good 
from  a  trend  standpoint  and  slightly  conservative. 

b.  Proprotor  Stability  with  Flapping  Locked  Out 

With  the  blade  flapping  locked  out  blade  loads  limited  the 
maximum  speed  at  which  the  model  could  be  tested.  Consequently, 
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the  instability  boundary  could  not  be  defined  experimentally. 
Figure  36  compares  the  maximum  velocity  reached,  with  the 
calculated  proprotor/pylon  stability  boundary.  Note  that  the 
zero  rpm  instability  is  effectively  suppressed  by  locking  out 
flapping. 

Correlation  between  the  calculated  and  measured  frequency  and 
damping  serves  as  a  measure  of  the  accuracy  of  the  program  with 
flapping  locked  out.  Figures  37  and  38  compare  the  calculated 
and  measured  frequency  and  damping  versus  airspeed  for 
£2R  =  600  feet  per  second  and  zero  feet  per  second,  respectively. 
Note  that  with  flapping  locked  out  the  model  is  more  stable 
than  it  was  with  the  blade  free  to  flap.  However,  care  must 
be  taken  in  concluding  that  a  hingeless  rotor  type  would 
necessarily  exhibit  the  same  stability.  As  noted  earlier, 
the  flapping  natural  frequency  of  the  blades  with  flapping 
locked  out  was  considerably  higher  than  that  typical  of  hingeless 
rotors.  Since  the  flapping  natural  frequency  is  known  to  be  an 
important  consideration  in  proprotor/pylon  instability  this 
model  characteristic  must  be  taken  into  account. 15 

c .  Instability  Frequency  and  Mode  Shape  Correlation 

Another  measure  of  program  accuracy  is  instability  frequency 
and  mode  shape  correlation.  Figure  39  compares  the  measured 
and  calculated  frequencies,  the  1:1  line  being  perfect  correla¬ 
tion.  Some  of  the  differences  between  the  calculated  and 
measured  frequency  is  due  to  the  use  of  calculated  model  mass 
and  stiffness  values  in  the  calculations.  These  do  not  give 
exactly  the  same  natural  frequencies  as  measured  in  a  ground 
vibration  test  of  the  model  (see  Appendix  I) . 

The  calculated  instability  mode  shapes  were  also  compared 
to  the  measured  instability  mode  shapes,  and  found  to  be  very 
close.  In  no  case  was  the  calculated  mode  of  instability  other 
than  that  measured. 

2 .  Flutter  and  Divergence  During  Blade  Folding 

A  model  configuration  which  was  felt  to  be  prone  to  flutter 
was  chosen  for  the  blade  folding  tests.  This  involved  using 
a  blade  feathering  restraint  with  a  stiffness  such  that  the 
blade  bending  natural  frequencies  traverse  the  wing  funda¬ 
mental  beam  and  chord  modes  as  the  blades  fold.  A  representa¬ 
tive  fold  hinge  stiffness  was  used. 

Prior  to  the  wind  tunnel  test  a  model  vibration  survey  was 
made.  The  correlation  was  found  to  be  poor  and  as  a  result, 
the  math  model  for  the  Folding  Dynamics  Analysis  was  expanded 
from  44  to  74  degrees  of  freedom.  This  modification  improved 
the  correlation  considerably. 
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During  the  wind  tunnel  test  it  was  not  possible  to  approach  a 
flutter  or  divergence  condition.  This  was  the  result  of  a 
limitation  on  maximum  velocity  imposed  by  fold  hinge  and  wing 
load  considerations.  Consequently,  frequency  and  damping 
variation  with  velocity  and  fold  angle  must  serve  to  indicate 
the  accuracy  of  the  Folding  Dynamics  program. 

a.  Natural  Frequency  Correlation 

Figures  40  and  41  show  the  calculated  and  measured  wing-pylon- 
blade  natural  frequency  variation  with  fold  angle,  respectively. 
Identification  of  the  modes  was  extremely  difficult  because  of 
the  high  coupling  between  the  blade  and  wing-pylon  motions. 

As  a  consequence,  the  modes  are  somewhat  arbitrarily  defined 
on  the  basis  of  the  frequency  of  the  uncoupled  modes.  The  modes 
are  identified  in  Table  IV. 


Mode 

1 


3 

4 


5 

6,7 

8 

9 

10 


Description 


TABLE  IV 

MODE  IDENTIFICATION 
Remarks 


Wing  Fundamental  Beam 
Bending 

Mast  Torsion 


Wing  Fundamental  Chord 
Bending 

Blade  Beam  Bending 


This  mode  changes  to  blade 
flapping  on  the  feathering 
hinge  as  the  blades  are 
folded. 


This  mode  couples  with  mast 
torsion  as  the  blades  are 
folded.  The  blade  bending 
is  in-phase. 


Blade  Chord  Bending 
Blade  Beam  Bending 
Blade  Chord  Bending 

Wing  Fundamental  Torsion  - 

Blade  Second  Beam  Bending  - 


Blades  out-of-phase 
Blades  out-of-phase 
Blades  in-phase 
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The  general  trends  of  the  calculated  frequency  as  the  blades 
are  folded  is  in  agreement  with  the  measured  trends. 
(Unfortunately  many  of  the  modes  could  not  be  identified 
adequately  during  the  vibration  survey  and  are  thus  not  shown 
for  all  fold  angles  in  Figure  40.) 

Frequency  and  damping  was  measured  versus  fold  angle  and 
velocity.  Figure  42  shows  the  measured  frequency  and  damping 
variation  with  velocity.  The  calculated  characteristics  are 
given  for  both  the  Feather-Stop  Analysis  (ARAP06)  and  the 
Folding-Dynamics  Analysis  (DFAL17) .  The  higher  damping 
calculated  with  the  latter  is  due  to  neglecting  stall  in  the 
blade  root  region.  ARAP06  included  stall  effects.  Adjustment 
of  the  strip  width  in  the  root  region  to  compensate  for  stall 
would  reduce  the  damping  calculated  with  DFAL17. 

The  correlation  with  regard  to  fold  angle,  shown  in  Figure  43, 
shows  the  same  damping  trend  for  both  analyses. 

C.  DYNAMIC  RESPONSE 

During  the  wind  tunnel  test  the  steady-state  response  of  the 
model  was  determined  as  a  function  of  proprotor  rpm.  At  each 
rpm  from  zero  to  350,  in  increments  of  10  rpm,  the  pylon 
accelerations,  wing  and  blade  loads,  and  flapping  were 
recorded.  Then  the  feathering  rate  was  progressively  increased 
to  the  power  limit  of  the  collective  pitch  control  motor. 
Testing  time  limitations  permitted  only  a  selected  number  of 
speed  and  angle  of  attack  combination  to  be  evaluated. 
Consequently  only  limited  experimental  data  are  available  to 
indicate  the  influence  of  airspeed  and  angle  of  attack. 

During  the  correlation  study  the  system  steady-state  response 
was  first  calculated  by  making  runs  at  selected  proprotor  rpm 
and  allowing  the  system  to  stabilize.  The  system  steady-state 
response  was  then  determined  from  the  calculated  time  history. 
These  calculations  were  made  early  in  the  correlation  study 
and  showed  the  computer  representation  of  wing-blade  inter¬ 
ference  to  be  in  error.  After  the  program  was  corrected,  the 
steady-state  response  was  not  recalculated  because  of  the  large 
amount  of  computer  time  required.  However,  extrapolation  of 
the  response  at  various  rates  of  feathering  to  zero  rate , 
indicates  reasonable  correlation  with  the  steady-state  response 


95 


A  flight  condition  felt  to  be  typical  of  the  FPR  transition 
mode  was  chosen  for  the  model  test;  an  airspeed  of  180  knots 
and  a  shaft  angle  of  attack  of  zero  degrees.  The  wing  angle 
of  attack  was  3.6  degrees,  less  than  that  required  for  lg 
flight  at  80  knots.  Flaps,  which  were  not  available  on 
the  model,  would  be  used  to  achieve  the  required  lift 
coefficient  while  maintaining  zero  shaft  angle  of  attack. 

Figures  44  and  45  show  the  steady-state  pylon  vibration  at 
the  blade  passage  frequency  and  harmonics  of  the  blade 
passage  frequency,  respectively,  versus  rpm.  Note  that  the 
response  at  the  blade  passage  frequency  is  dominant.  The 
one-third  blade  passage  frequency  (one  per  rev)  response 
shown  for  the  wing  beam  and  chord  modes  in  Figure  45  is 
due  to  model  rotor  out-of-balance  and  out-of-track  and  is 
felt  to  be  higher  than  representative  of  full-scale  since 
small  scale  model  tolerances  are  generally  proportionally 
larger  than  full-scale. 

Figure  46  shows  the  measured  steady-state  blade  loads  versus 
rpm.  Noteworthy  is  the  absence  of  large  buildups  in  load 
when  the  rotor  rpm  is  such  that  a  rotor  blade  natural 
frequency  is  in  resonance  with  aerodynamic  excitation. 

As  noted  earlier,  the  basic  configuration  involved  a 
gimballed  rotor  with  moderate  hub  restraint.  For  the 
steady-state  response  run  the  mast  was  aligned  close  to 
parallel  to  the  test  section  centerline.  (The  measured 
mast  pitch  was  0.6  degrees;  the  mast  yaw  was  not  measured.) 
Nevertheless,  large  amplitude  flapping  was  encountered. 
Figure  47  shows  the  total  flapping  of  the  red  blade  as  a 
function  of  proprotor  rpm. 

1.  Response  Correlation 

The  influence  of  blade  feathering  rate  was  evaluated  by 
progressively  increasing  the  feathering  rate  up  to  the 
collective  motor  limit.  The  slowest  rate  attained  was 
approximately  2  degrees  per  second,  which  caused  the  rotor 
to  stop  in  30  seconds  (0.75  degrees  per  second  and  82 
seconds  respectively,  full-scale).  The  highest  rate  was 
approximately  150  degrees  per  second,  which  caused  the  rotor 
to  stop  in  0.4  seconds  (55  degrees  per  second  and  1.1 
seconds  respectively,  full-scale). 
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Several  observations  during  the  tests  should  be  noted  at 
this  times  (1)  once  the  zero  rpm  collective  pitch  setting 
was  determined,  and  the  collective  limit  set,  there  was  no 
requirement  to  change  the  setting  for  different  tunnel 
speeds  or  shaft  angles.  (2)  The  rotor  blades  were  virtually 
self-indexing.  The  rotor  always  stopped  with  one  blade 
pointed  outward  from  the  wing.  This  is  a  result  of  the 
wing-blade  aerodynamic  interference. 

Figure  48  shows  a  time  history  of  model  accelerations  and 
loads  during  a  slow  stop  (20  seconds  full-scale) .  The 
buildups  caused  by  traversing  resonant  rpm  are  noted.  For 
comparison  time  histories  of  a  fast  stop  (1.5  seconds  full- 
acale)  and  of  a  fast  start  (1.1  seconds)  are  shown  in  Figures 
49  and  50,  respectively.  The  response  during  unfeathering 
of  the  proprotor  was  observed  to  be,  in  general,  si .ghtly 
less  than  during  feathering. 

A  computed  time  history  of  a  fast  stop  (1  second,  full- 
scale)  is  shown  in  Figure  51.  Note  that  the  general  character 
of  the  response  is  similar  to  that  of  Figure  49.  The 
primary  difference  between  the  calculated  and  measured 
response  is  the  absence  of  one-per-rev  response  in  the 
calculated  history. 

Sensitivity  to  rate  was  evaluated  by  plotting  the  response 
occurring  at  the  aystem  natural  frequencies,  versus  feather¬ 
ing  rate.  Figures  52  through  56  compare  the  calculated  and 
measured  rate  sensitivity  to  the  inverse  of  the  time  to 
feather  or  unfeather  the  rotor.  The  feathering  rate  is 
related  to  the  inverse  of  the  time  to  feather,  AT,  as  follows: 


Feather  Rate  Degrees  Per  Second  =  (108.8°) (1/AT)  (222) 


In  general,  the  correlation  between  theory  and  measured 
response  is  excellent. 

Several  noteworthy  observations  can  be  made  from  the  feathering 
rate  data:  (1)  Even  at  the  fastest  feathering  rates,  blade 
stall  was  not  evident.  Figure  54  shows  that  the  wing  chord- 
wise  peak  transient  bending  moment  varied  almost  linearly  with 
rate.  Had  the  blades  stalled  at  the  higher  rates  the  transient 
moment  would  fall  off.  (2)  The  valley  in  the  wing  chordwise 
moment  was  the  only  indicator  of  an  "optimum  feathering  rate" 
observed  during  the  test.  Based  on  the  chordwise  moment  valley, 
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the  optimum  feather  time  would  be  about  2.75  to  3  seconds. 
Figures  52  and  53  indicate  that  faster  feathering  rates  do 
not  provide  a  significant  reduction  in  other  system  responses. 
(3)  Increasing  the  feathering  rate  reduces  the  resonant  buildup 
in  blade  loads.  But  very  fast  rates  produce  an  impulsive 
response  in  the  wing  chordwise  bending  moment.  This  results 
from  the  impulsive  nature  of  the  proprotor  aerodynamic  loading 
at  the  higher  rates.  However,  over  the  range  of  rates  tested, 
blade  loads  were  very  low  and  were  not  a  limitation  on  the 
testing.  (4)  Flapping  freedom  must  be  restrained  prior  to 
feathering  the  blades.  Figure  56  shows  that  faster  rates 
increase  the  transient  flapping  during  proprotor  stopping. 

The  characteristics  may  be  explained  by  considering  the  way  in 
which  the  blade  flap-damping  varies  during  feathering.  The 
flap-damping  coefficient  is  given  as? 

C  -  &  cos*3/4  (223) 


where 


<J>3/4  =  TAN 


(224) 


the  inflow  angle  at  the  blade  3/4  radius,  y  is  blade  lock  number. 


As  the  blade  is  feathered  the  inflow  angle  varies  approximately 
in  proportion  to  the  blade  collective  pitch.  According  to 
equation  223  the  flap  damping  then  decreases  as  the  blade  is 
feathered  and  increases  as  the  blade  is  unfeathered.  And 
according  to  Reference  16,  this  time  variation  in  the  damping 
parameter  of  a  single  degree-of-freedom  system  will  cause  the 
apparent  damping  to  be  less  when  the  damping  decreases  with 
time,  and  more  when  the  damping  increases  with  time.  Since 
flapping  amplitude  is  controlled  by  the  flap  damping,  the 
characteristics  shown  in  Figure  56  are  to  be  expected. 

2 .  Other  Factors  Influencing  Dynamic  Response  Characteristics 

(1)  Flapping  Lockout  -  Prior  to  this  study,  it  was  thought 
that  flapping  could  be  locked  out  during  the  feather/stop 
sequence,  at  about  25-30  percent  of  the  hover  rpm.  The  effect 
of  feathering  rate  on  the  dynamic  flapping  response  (Figure  56) 
and  practical  considerations,  however,  indicate  flapping  should 
be  locked  out  prior  to  feathering. 

Locking  out  flapping  did  not  significantly  influence  the  wing 
and  pylon  dynamic  response.  Both  test  data  and  calculated 
pylon  responses  show  the  acceleration  to  be  close  to  that  with 
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the  flapping  free.  The  blade  tip  deflections  were  not  measured 
during  the  model  test,  but  were  observed  to  be  small. 

The  principal  influence  of  locking  out  flapping  freedom  is  a 
mast  angle  of  attack  limit  resulting  from  blade  loads.  (Locking 
out  the  flapping  freedom  increases  the  one  per  rev  out-of-plane 
loads  since  flapping  acts  to  relieve  the  one-per-rev  airloads.) 
Figure  57  shows  the  measured  angle  of  attack  limit  established 
by  the  model  test.  Note  that  this  limit  would  only  apply 
during  the  feather/stop  sequence.  The  angle  of  attack  limit 
with  the  gimbal  free  is  much  higher. 

(2)  Angle  of  Attack  -  Shaft  angle  of  attack  has  a  significant 
influence  on  the  dynamic  response.  Figure  58  shows  a  time 
history  of  a  blade  feathering  at  a  shaft  angle  of  attack  of 
5.6  degrees.  The  most  apparent  difference  with  the  zero  angle- 
of-attack  history,  Figure  49,  is  the  impulsive  response  at  the 
wing  beam  natural  frequency,  seen  in  the  wing  beam  bending 
moment  and  torsion  traces.  This  is  the  result  of  the  sudden 
buildup  in  normal  force  as  the  proprotor  stops.  (The  normal 
force  is  essentially  proportional  to  the  inflow  angle.)  Figure 
59  shows  how  the  beam  response  varies  with  feathering  rate. 

Note  that  while  a  slow  feather  rate  is  indicated  to  be  desirable, 
a  shaft  angle  of  attack  of  5.6  degrees  during  the  feather/stop 
sequence  is  considered  of  transient  nature  only.  Operationally, 
the  feather/stop  sequence  should  be  made  with  a  shaft  angle  of 
attack  near  zero  degrees. 

Another  noteworthy  observation  can  be  made  from  Figure  58.  The 
blade  loads  buildup  remain  moderate  even  at  large  shaft  angles 
of  attack. 

(3)  Airspeed  -  Three  airspeeds  were  investigated  during  the 
feather/stop  test.  From  the  test  data  it  appears  that  the 
response  amplitude  increased  essentially  linearly  with  airspeed 
for  a  constant  wing  lift  coefficient.  This  is  logical,  since 
the  system  damping  is  proportional  to  velocity,  whereas  the 
airload  forcing  function  is  proportional  to  velocity  squared. 

At  resonance  the  response  amplitude  is  proportional  to  the 
ratio  of  the  forcing  function  divided  by  the  damping.  Hence 
the  response  would  be  expected  to  be  proportional  to  velocity. 

In  Section  IV,  the  influence  of  airspeed  on  response  is  studied 
in  more  detail.  When  the  lower  wing  lift  coefficient  with 
increasing  airspeed  and  consequent  lower  wing  upwash  is 
accounted  for,  the  net  effect  is  a  negligible  influence  of 
airspeed  on  response. 
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3 .  Response  During  Blade  Folding 

In  the  blade  folding  configuration  the  model  was  tested  over 
the  range  of  airspeed,  fold  angle,  and  shaft  angle-of-attack 
indicated  in  Figure  60.  Wing  and  blade  load  considerations 
limited  the  envelope  investigated. 

As  noted  earlier  in  this  section,  there  was  no  evidence  of 
flutter,  divergence,  or  buffet  during  the  testing.  The  data 
presented  in  this  report  therefore  deals  only  with  the  steady 
load  correlation. 

The  airstream  oscillation  capability  of  the  16-foot  TDT  wind 
tunnel  was  employed  to  investigate  the  frequency  response 
characteristics  versus  fold  angle.  Unfortunately  the  measured 
airstream  flow  angularity  data  appears  questionable.  The 
flow  angle  data  do  not  agree  with  previous  flow  angle  measure¬ 
ments  and  the  zero  frequency  response  per  degree  of  indicated 
flow  angle  is  too  low.  Consequently,  no  correlation  with 
measured  frequency  response  is  presented  here.  To  indicate 
the  correlation  between  theory  and  measured  loads  data,  the 
wing  beam  bending  at  15  percent  semispan,  blade  beam  and  chord 
bending  at  35  percent  radius,  and  the  blade  spindle  chord 
moments  and  torques  were  used. 

Figures  61,  62,  and  63  compare  the  calculated  and  measured 
loads  for  the  wing  beam  versus  fold  angle,  velocity,  and  shaft 
angle  of  attack.  The  correlation  is  felt  to  be  very  good. 

(The  wing  lift  curve  slope  used  in  the  analysis  is  based  on 
measured  data  for  the  model.) 

The  calculated  and  measured  spindle  torque  is  shown  as  repre¬ 
sentative  of  the  blade  load  correlation.  Figures  64,  65,  and 
66  show  how  the  spindle  torque  varies  with  fold  angle,  angle 
of  attack,  and  airspeed.  The  peak  at  a  fold  angle  of  35  to 
45  degrees  is  evident.  A  significant  load  characteristic 
during  blade  folding  is  evident  in  Figure  66 .  The  blade  loads 
variation  with  airspeed,  at  partial  fold  angles,  does  not 
follow  the  usual  velocity-square  ratio.  This  is  the  result  of 
deflections  of  the  blade  and  control  system  reducing  the  blade 
angle  of  attack.  This  is  similar  to  the  gust  load  alleviation 
characteristic  of  swept  wings. 

Another  characteristic  noted  in  the  tests  is  higher  than  pre¬ 
dicted  blade  loads  at  fold  angles  near  90  degrees  as  shown  in 
Figure  64.  This  is  theorized  as  being  caused  by  the  wing-tip 
vortex  tending  to  lift  the  outboard  blade  as  the  blade  folds 
close  to  the  vortex.  The  analysis  does  not  include  a  tip 
vortex  representation,  hence  the  lower  calculated  loads. 

Figure  67  further  illustrates  this  effect. 
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NOT  REPROnUCIBLE 


Figure  24 


Model  Used  to  Obtain  Wing- Blade 
Aerodynamic  Interference  Data 


V 


YELLOW  BLADE  ROOT  BENDING  MOMENT/q 


Figure  25.  Correlation  Between  Theory  and  Measured 
Blade  Root  Inplane  Bending  Moment 
(Model  Parameters) 
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FROPROTOR  THREE- PER-REV  H-FORCE/q 


MAST  ANGLE -OF-ATTACK  -  DEGREES 


Figure  26,  Correlation  Between  Theory  and  Measured 

Effect  of  Shaft  Angle  of  Attack  on  Oscil¬ 
latory  Normal  Force  (Model  Parameters) 
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PROPROTOR  THREE-PER-REV  H  FORCE/q 
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RED  BLADE  AZIMUTH  -  DEGREES 


Figure  28.  Correlation  Between  Theory  and  Measured 
Shapes  of  Oscillatory  Normal  Force 
(Model  Parameters)  . 
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b.  Folding  Configuration 


Figure  29.  Model  Used  to  Obtain  Dynamic 
Stability  Data 
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D.  Time  History  of  Proprotor  Pylon 
Instability  at  Wing  Beam  Natural 
Frequency 
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Figure  31.  Example  of  Measured  Wing-Beam  Frequency 
and  Damping  Variation  with  Airspeed, 

300  rpm. 
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Figure  32.  Correlation  Between  Theory  and  Measured 
Proprotor  Stability  Boundary,  Gimbal 
Free . 
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igure  33.  Time  Hlatory  of  Zero  RPM  Flapping 
Instability 
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Figure  34 .  Correlation  Between  Theory  and  Measured 
Effect  of  Hub  Restraint  on  Zero  RPM 
Flapping  Stability. 
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Figure  36.  Airspeed  -  RPM  Envelope  Explored  With 
Girabal  Flapping  Freedom  Locked  Out. 
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Figure  37.  Correlation  Between  Theory  and  Measured 
Wing  Beam  and  Chord  Damping,  Gimbal 
Flapping  Freedom  Locked  Out,  300  RPM. 
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DAMPING  RATIO  -  C/C 


SIMULATED  TRUE  AIRSPEED  -  KNOTS 


Figure  38.  Correlation  Between  Theory  and  Measured 

Wing  Beam  and  Chord  Model  Damping,  Gimbal 
Flapping  Locked  Out,  Zero  RPM. 
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CALCULATED  FREQUENCY 


Figure  39.  Correlation  Between  Theory  and  Measured 
Instability  Frequency. 
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NATURAL  FREQUENCY 


THE  MODES  ARE  IDENTIFIED  IN 
_  TABLE  IV.  PAGE  94  IN  THE  TEXT 


FOLD  ANGLE  -  DEGREES 


Figure  41.  Calculated  Blade  Natural  Frequencies  Versus 
Blade  Fold  Angle . 
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FREQUENCY  -  CPS 


SIMULATED  TRUE  AIRSPEED  -  KT 


Figure  42.  Correlation  of  Theory  With  Measured 

Frequency  and  Damping  Versus  Velocity 
Zero  Fold  Angle. 
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DAMPING  RATIO  -  C/C 


Figure  42.  Concluded. 
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Figure  43.  Correlation  Between  Theory  and  Measured 
Frequency  and  Damping  Versus  Fold  Angle, 
211  Knots  Simulated  Airspeed. 
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Figure  43.  Concluded. 
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PYLON  VIBRATION 


Figure  44.  Measured  Steady  State  Pylon 
Vibration  at  Blade  Passage 
Frequency. 
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PYLON  VIBRATION 


BENDING  MOMENT  -  IN-LB  X  10" 


Figure  46.  Measured  Steady  State  Blade  Loads  at  Blade 
35-Percent  Radius. 
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BLADE  FLAPPING  -  DEG 


Figure  47.  Measured  Steady-State  Flapping  Response. 
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Time  History  of  Model  Vibration  During  Slow 
Unfeathering  of  the  Blades* 
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WIND  CHORD  (28%) 
582  IN-LB/IN 


Figure  50.  Time  History  of  Model  Response  to  Rapid  Unfeathering 
(AT  =  0.4  Seconds),  175  Knots  Simulated  Airspeed. 


MAXIMUM  VERTICAL  ACCELERATION  AT  CONVERSION  AXIS 


1/TIME  TO  FEATHER/UNFEATHER  -  1/SEC 


Figure  52.  Correlation  Between  Theory  and  Measured 
Response  at  Wing  Torsional  Natural 
Frequency. 
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Figure  53.  Correlation  Between  Theory  and  Measured 
Response  at  Wing  Beam  Natural  Frequency. 
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BENDING  MOMENT  -  IN- LB  X  10 
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I/TIME*  TO  FEATHER  -  1/SEC 


Figure  55.  Correlation  Between  Theory  and  Measured 
Blade  Beam  Bending  Moments  During 
Feathering . 
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TOTAL  FLAPPING  -  DEG 


Figure  56.  Correlation  Between  Theory  and  Measured 
Effect  of  Feathering  Rate  on  Flapping. 
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MAST  ANGLE  OF  ATTACK 


Figure  57.  Measured  Mast  Angle  of  Attack  Limit  with 
Flapping  Freedom  Locked  Out,  300  RPM. 
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Figure  58.  Time  History  of  Model  Vibration  During  Feathering 
with  Shaft  Angle  of  Attack  of  5.6°. 
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Figure  59.  Correlation  Between  Theory  and  Measured 
Response  at  Wing  Beam  Natural  Frequency 
With  Shaft  Angle-Of-Attack  Of  5.6  Degrees. 
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Figure  60.  Angle-of-Attack  and  Airspeed  Envelope 
Explored  in  Blade  Folding  Test. 
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BENDING  MOMENT  -  IN-LB 


BENDING  MOMENT  -  IN-LB 
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Figure  62.  Correlation  Between  Theory  and  Measured 
Wing  Root  Beam  Bending  Moment  Variation 
With  Airspeed. 
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Correlation  Between  Theory  and  Measured 
Wing  Root  Beam  Bending  Moment  Variation 
With  Shaft  Angle  of  Attack. 
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Figure  63. 


Figure  64.  Correlation  Between  Theory  and  Measured 
Pitch  Change  Torque,  Variation  with 
Fold  Angle . 
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Figure  65.  Correlation  Between  Theory  and  Measured 
Pitch  Change  Torque,  Variation  with 
Shaft  Angle  of  Attack 
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SECTION  IV 


DYNAMIC  ANALYSIS  AND  PARAMETRIC  STUDY 

A  dynamic  analysis  of  a  representative,  66,000-pound  gross 
weight,  FPR  VTOL  design  was  made  to  indicate  dynamic  stability 
and  response  characteristics  during  the  feathering  and  fold¬ 
ing  of  the  blades  of  that  class  aircraft.  It  also  serves  as 
the  basis  for  a  parametric  study  to  identify  design  factors 
which  control  the  dynamic  stability  and  response  during  blade 
feathering/folding . 

The  design  used  for  this  analysis  is  one  developed  by  Bell 
Helicopter  Company  under  contract  to  the  Air  Force  Flight 
Dyanmics  Laboratory . I?  It  carries  the  Bell  Helicopter  design 
designation  D270A.  Descriptive  data  and  major  component 
drawings  are  given  in  Appendix  II.  The  data  provided  are 
sufficient  to  allow  the  reader  to  conduct  an  independent 
analysis . 

A.  D270A  DYNAMIC  ANALYSIS 

The  following  D270A  feather/fold  dynamic  characteristics  were 
identified  in  the  course  of  this  analysis. 

-  proprotor/pylon  stability  boundaries 

-  wing/blade  flutter  and  divergence  boundaries 

-  dynamic  response  during  blade  feathering 
gust  and  maneuver  limits 

Although  the  analysis  was  focused  on  the  feather/fold  sequence, 
the  proprotor  stability  boundary  in  proprotor  cruise  mode  and 
system  natural  frequencies  in  all  modes  were  determined. 

A  number  of  D270A  design  features  are  incorporated  on  the 
basis  of  dynamic  considerations.  These  will  be  discussed 
first. 

1 .  Design  Features  Influenced  by  Dynamic  Considerations 
a .  Wing 

The  wing  thickness  and  planform  are  selected  primarily  on  the 
basis  of  proprotor  dynamic  stability  considerations.  A  wing 
thickness  ratio  of  0.18,  constant  root  to  tip,  combined  with  a 
tapered  chord  section  and  spars  located  at  5-  and  50-percent 
chord,  provide  a  large  torque  box.  All  of  the  required 
bending  material  is  concentrated  in  the  skins  and  spar  webs 
to  maximize  the  effective  torque  box  skin  thickness.  The 
skin  gage  in  the  outboard  50-percent  span  is  greater  than 
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required  for  strength,  further  increasing  the  wing  torsional 
stiffness.  The  high  torsional  stiffness  is  reflected  in  a 
tip  torsional  spring  rate  of  6.25  x  10<5  in-lb/radian.  The 
wing  bending  stiffnesses,  beamwise  and  chordwise,  were 
established  by  strength  requirements. 

A  wing  forward  sweep  of  six  degrees  was  chosen  to  minimize  the 
distance  from  the  tip-path-plane  of  the  proprotor  to  the  wing 
elastic  axis,  thus  minimizing  the  destabilizing  influence  of 
the  proprotor  shear  forces,  while  maintaining  adequate  flapping 
clearance . 

b .  Pylon 

The  pylon  is  rigidly  attached  to  the  wing  tip  to  take  advantage 
of  the  wing  torsional  stiffness.  The  shaft  is  as  short  as 
possible,  consistent  with  flapping  clearance  requirements. 

c .  Proprotor 

A  gimbaled  hub,  three-bladed ,  semirigid  type  proprotor  was 
chosen  on  the  basis  of  blade  loads,  blade  motion  stability, 
and  mechanical  stability  considerations.  When  the  gimbal 
freedom  is  locked  out,  as  it  is  before  feathering  the  blades, 
this  proprotor  configuration  approaches  a  propeller  in  stiff¬ 
ness  characteristics. 

The  proprotor  blades  are  extremely  strong  and  stiff  to  resist 
oscillatory  loads  in  helicopter  mode  and  gust  and  maneuver 
loads  in  proprotor  mode.  The  blades  are  preconed  2%  degrees 
and  employ  positive  pitch-flap  coupling  (<S3  =  -25  degrees)  to 
preclude  blade  motion  instability. t 

2 .  Feather-Fold  Flight  Procedure 

For  this  analysis  it  was  assumed  the  following  flight  pro¬ 
cedures  will  be  followed  when  feathering  and  folding  the 
blades : 


-  The  airspeed  will  be  typically  between  150  and 
200  knots.  In  this  analysis,  175  knots  were 
used  as  a  representative  airspeed. 

The  aircraft  attitude  will  be  trimmed,  using  the  wing 
flaps,  to  zero  degrees  fuselage  angle  of  attack. 

This  places  the  shaft  at  zero  angle  of  attack  to 
minimize  the  variation  in  proprotor  steady  aerodynamic 
forces  and  moments  during  feathering  and  folding. 
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Flapping  will  be  automatically  locked  out  prior  to 
feathering  the  blades,  and  restored  only  after  the 
blades  are  unfeathered. 

This  operational  procedure  was  established  in  light  of  the 
results  of  exploratory  tests  and  studies,  including  the 
wind  tunnel  test  conducted  under  this  program. 

3 .  Airframe  and  Proprotor  Natural  Frequencies 

a.  Airframe  Natural  Frequencies 

The  natural  frequencies  of  the  wing-pylon-fuselage  system 
were  calculated  using  a  finite  element  (beam  representation 
only)  vibration  analysis.  The  empennage  and  proprotor  were 
assumed  to  be  lumped  rigid  masses. 

Figures  68  and  69  show  the  symmetric  and  asymmetric  natural 
frequencies  versus  mast  conversion  angles.  The  changes  in 
natural  frequency  when  the  proprotor  gimbal  is  locked  out 
prior  to  and  during  the  feather/fold  sequence  are  also 
indicated.  Table  V,  below,  identifies  the  modes. 

TABLE  V 


AIRFRAME  NATURAL  FREQUENCY  IDENTIFICATION 


Mode  No. 

Symmetric 

Asymmetric 

1 

Wing  1st 

Beam 

Wing  1st 

Torsion 

2 

Wing  1st 

Chord 

Wing  1st 

i 

Chord 

3 

Wing  1st 

Torsion 

Fuselage 

Lateral 

4 

Fuselage 

Vertical 

Wing  1st 

Beam 

5 

Wing  2nd 

Beam 

Wing  2nd 

Torsion 

6 

Wing  2nd 

Torsion 

Wing  2nd 

Beam 

7 

Wing  2nd 

Chord 

Wing  2nd 

Chord 

Several  resonance  problems  are  indicated:  1)  The  proximity 
to  one-per-rev  of  the  first  symmetric  chord  and  first 
asymmetric  torsion  modes  and  resonance  during  partial  con¬ 
version  of  the  first  asymmetric  chord  mode,  indicate  that 
rpm  may  have  to  be  scheduled  during  conversion  to  avoid 
prolonged  operation  at  resonance.  However,  this  may  not  be 
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as  serious  a  problem  as  it  appears  since  recent  full-scale 
tests  indicate  operation  in  and  near  one-per-rev  resonance  is 
possible  with  proprotors .19  2)  The  wing  second  beam  and 

torsion  modes  are  indicated  to  be  in  three-per-rev  resonance. 
These  modes  have  nodal  points  near  the  hub,  thus  they  are 
weakly  excited  by  the  proprotors.  These  modes  can  probably 
be  moved  out  of  resonance  by  spanwise  relocation  of  the 
engines . 

b.  Proprotor  Natural  Frequencies 

The  proprotor  blade  natural  frequencies  were  calculated  using 
a  Myklestad-type  analysis  for  a  rotating-twisted  beam.  The 
analysis  includes  the  coupling  between  beamwise  and  chord- 
wise  deflections  that  results  from  the  built-in  twist  and 
collective  pitch.  Good  correlation  has  been  achieved  with 
the  measured  frequencies  of  many  rotor  designs,  including 
three-bladed,  semirigid  rotors  such  as  those  of  the  Model 
D270A. 

Figures  70  and  71  show  the  calculated  natural  frequencies  as 
a  function  of  rpm  and  collective  pitch.  Figure  72  shows  the 
change  in  the  frequencies  when  gimbal  freedom  is  locked  out. 

The  frequencies  are  presented  here  in  terms  of  collective  and 
cyclic  modes.  The  collective  modes  are  those  excited  by  air¬ 
loads  whose  frequency  per  revolution  is  an  integer  multiple  of 
the  number  of  blades.  Noninteger  multiple  harmonic  airloads 
excite  the  cyclic  modes.  The  coupled  natural  frequencies  are 
noted  by  crosses  (x)  and  diamonds  (O) .  The  crosses  denote 
modes  whose  largest  deflection  is  normal  to  the  plane  of 
rotation.  The  solid  lines  denote  the  frequency  of  an  untwisted 
blade  at  zero  collective  pitch  (i.e.,  uncoupled  beamwise  and 
chordwise  frequencies)  and  are  provided  for  reference.  By 
comparing  the  uncoupled  and  coupled  frequencies,  the  effect 
of  built-in  twist  on  the  blade  frequencies  is  apparent. 

Past  experience  of  the  design  and  testing  of  the  three-bladed 
semirigid  rotors  has  shown  that  frequency  placement  of  the 
first  inplane  (cyclic)  and  second  beamwise  (collective)  modes 
poses  the  most  severe  requirements.  The  first  inplane  mode 
must  be  sufficiently  removed  from  one-per-rev  to  avoid  high 
one-per-rev  loads,  but  cannot  be  too  high  or  two-per-rev  loads 
will  be  a  problem.  A  frequency  of  1.5-per-rev  results  in  the 
lowest  oscillatory  loads.  However,  frequencies  as  low  as 
1.25-per-rev  are  acceptable.  Meeting  this  requirement  can  be 
difficult  in  proprotors  because  of  the  wide  rpm  and  collective 
range  required  for  efficient  operation.  The  second  beam  mode 
must  be  located  above  three-per-rev  to  avoid  high  loads  and 
airframe  vibration.  In  helicopters,  keeping  this  mode  out 
of  resonance  has  been  a  problem;  however,  it  is  less  of  a 
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problem  with  proprotors  because  of  the  thicker  and  stiffer 
blade  root  sections  required  for  static  strength. 

The  frequency  location  of  the  major  blade  modes  of  the  D270A 
proprotor  is  satisfactory:  the  first  inplane  mode  varies  from 
1.8  to  1.38-per-rev  in  helicopter  mode  (134  rpm)  and  from  1.5 
to  1.85-per-rev  in  airplane  mode  (225  rpm).  The  second  beam 
mode  is  about  4.0-per-rev  in  both  helicopter  and  airplane  mode. 
Thus,  these  two  modes  are  well  located.  Close  proximity  to 
4-per-rev  resonance  is  indicated  for  the  second  cyclic  mode 
at  high  pitch  in  airplane  mode  (225  rpm)  and  6-per-rev  resonance 
for  the  third  collective  mode  at  high  collective  pitch.  How¬ 
ever,  low  airload  excitation  at  these  resonances  would  be 
expected . 

When  gimbal  freedom  is  locked  out,  the  second  cyclic  mode  is 
in  resonance  with  2-per-rev  at  tip  speeds  above  750  feet  per 
second.  Consequently,  when  flapping  is  locked  out,  the 
proprotor  may  have  to  be  restricted  to  250  rpm,  maximum.  The 
collective  modes  are  not  affected  by  locking  out  flapping. 

The  frequency  of  the  blade  first  torsional  mode  was  not  calcu¬ 
lated  but  is  estimated  to  be  located  at  4.5-per-rev.  This 
mode  is  rigid  body  blade  pitching  based  on  the  control  system 
flexibility.  The  second  torsional  frequency  which  involves 
blade  torsional  deflection  is  much  higher  in  semirigid  rotor 
types  due  to  the  inherent  torsional  stiffness  of  the  blades. 

The  design  blade  loads  for  the  D270A  were  scaled  from  loads 
obtained  in  tilt-proprotor  design  studies.  These  past  studies 
have  shown  that  for  tilt-proprotors  two  flight  conditions 
impose  the  designing  loads.  For  oscillatory  loads,  the 
maximum  level  flight  airspeed  in  helicopter  mode  is  the  most 
severe.  These  must  be  below  the  endurance  limit  to  achieve  a 
reasonable  fatigue  life  for  the  proprotor.  Design  limit  loads 
are  generally  established  by  a  vertical  gust  encounter  in  air¬ 
plane  mode. 

4 .  Proprotor  Stability,  Flutter,  and  Divergence  Character¬ 
istics 


a.  Proprotor/Pylon  Stability 

The  proprotor/pylon  stability  (whirl  flutter)  boundaries  were 
calculated  for  proprotor  mode  as  well  as  during  blade 
feathering.  In  proprotor  mode  the  gimbal  is  restrained  with 
a  moderate,  isotropic,  hub  restraint  of  3850  ft-lb/deg  of  hub 
flapping.  Before  feathering  the  blades,  the  hub  gimbal 
freedom  is  locked  out  with  a  hydraulic  flapping  lockout  (see 
Drawing  D270A-410-002 ) . 
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The  mathematical  model  for  the  proprotor  stability  analysis 
included  the  wing  fundamental  beam,  chord,  and  torsion  modes 
and  the  pylon  yaw  (second  chord  mode) .  Both  symmetric  and 
asymmetric  modes  were  considered.  The  symmetric  modes  were 
found  to  be  the  less  stable. 

Figure  73  shows  the  proprotor/pylon  stability  boundary  in 
proprotor  mode  versus  rpm.  The  mode  of  instability  that 
first  becomes  unstable  is  the  symmetric  wing  chord  mode.  The 
symmetric  wing  beam  mode  becomes  unstable  at  a  slightly  higher 
airspeed.  The  root  motion  with  airspeed,  calculated  using  a 
Bell  Helicopter  Company  linearized  proprotor  stability  analysis, 
is  shown  in  Figure  74. 

The  stability  margin  requirement  of  1.15VL  as  specified  in  MIL-8870 
(ASG)  is  met  by  the  D270.  If  it  is  desired  to  increase  the  sta¬ 
bility  margin  in  proprotor  mode,  reducing  the  pitch-flap  coupling 
is  the  most  effective  means.  The  effect  of  reducing  the  pitch- 
flap  coupling  to  63  =  -15°  or  -20°,  is  shown  in  Figure  73.  This 
would  require  some  redesign  of  the  proprotor  rotating  controls 
but  could  be  accomplished  with  a  small  weight  increase. 

When  the  gimbal  flapping  freedom  is  locked  out  during  blade 
feathering,  the  stability  margin  is  increased.  Figure  75 
shows  the  proprotor  stability  boundary  with  the  flapping 
locked  out.  The  mode  of  instability  in  the  100  to  300  rpm 
range  is  that  of  the  wing  beam  mode  but  also  involves  wing 
chord  motion.  Figure  76  shows  the  root  motion  with  airspeed. 

The  beam  mode  in  this  case  corresponds  to  the  retrograde 
whirl  mode  of  classical  whirl  flutter  theory.20 

At  very  low  rpm,  divergence  of  the  wing  beam  mode  is  the  mode 
of  instability.  The  indicated  stability  margin  is  more  than 
ample. 


Flutter  and  Divergence 


Flutter  and  divergence  characteristics  during  blade  folding 
were  calculated  using  the  Folding-Dynamics  Analysis,  DFAL17 . 
Since  the  math  model  consists  of  a  cantilever  wing,  the 
divergence  boundary  in  this  case  is  somewhat  conservative. 


Figure  77  shows  the  D270  flutter  and  divergence  boundaries  as 
a  function  of  the  blade  fold  angle.  For  zero  fold  the  mode 
of  instability  is  wing  torsional  divergence,  a  result  of  the 
proprotor  aerodynamic  center  being  well  forward  of  the  wing 
elastic  axis.  Figure  78  shows  the  natural  frequency  variation 
with  airspeed  with  the  blades  folded.  While  the  mode  denoted 
wing  beam  bending  is  indicated  to  diverge,  inspection  of 
the  mode  shape  shows  increasing  wing  torsional  deflection  as 
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airspeed  increases.  As  the  blades  are  folded,  the  divergence 
speed  increases  rapidly  since  the  proprotor  aerodynamic 
center  moves  aft. 

The  flutter  boundary  indicated  in  Figure  77  is  a  weak,  wing 
chord  mode  flutter.  It  results  from  coupling  between  the 
blade  beamwise  bending  and  the  second  wing  chord  mode.  The 
second  chord  mode  contains  a  large  amount  of  nacelle  yawing 
coupling  the  mode  with  blade  beam  bending  at  large  fold  angles. 
Figure  79  shows  the  natural  frequency  variation  with  airspeed 
for  a  fold  angle  of  60  degrees.  Figure  80  shows  the  damping 
versus  airspeed.  Note  that  the  blade  natural  frequencies 
increase  rapidly  with  airspeed  and  coalesce  with  the  second 
chord  mode  at  a  speed  of  about  890  knots.  Since  aerodynamic 
damping  of  wing  chord  motions  is  small,  the  blade  drives  the 
chord  mode  unstable.  It  should  be  noted  that  only  a  small 
amount  of  damping  is  indicated  to  be  required  to  stabilize 
this  flutter.  Coupling  of  the  first  chord  mode  with  blade 
beam  bending  was  investigated  experimentally  as  discussed  in 
Section  III. 

The  blade  over  mass  balancing  required  for  helicopter  mode 
effectively  prevents  blade  flutter  during  folding. 

The  relatively  low  divergence  airspeed  suggested  that  aileron 
reversal  might  be  a  problem.  The  reversal  speed  was  estimated 
using  the  method  outlined  in  Reference  22  but  using  the  wing 
tip  torsional  spring  rate  rather  than  the  0.7  span  value  used 
for  conventional  aircraft.  This  gave  a  reversal  speed  of  320 
knots  for  the  zero-degree  fold  angle.  This  low  reversal  speed 
emphasizes  the  requirement  for  including  aeroelastic  effects 
in  stability  and  control  analysis  of  the  feather/fold  sequence. 

5 .  Dynamic  Response  During  Blade  Feathering 

a.  Math  Model  and  Method 

The  first  four  coupled  wing-pylon-fuselage  modes  are  included 
in  the  math  model  as  well  as  the  proprotor  first  collective 
and  first  and  second  cyclic  modes.  The  coupled  wing-pylon- 
fuselage  modes  consist  of  the  wing  fundamental  beam  and  chord 
bending  and  torsion  modes,  and  the  first  fuselage  vertical 
bending  mode.  As  indicated  in  Figures  68  and  69,  the  wing 
second  beam  bending  and  torsion  modes  are  also  transited  by 
three-per-rev  during  blade  feathering.  However,  these  modes 
are  relatively  weak  compared  to  the  fundamental  modes  and 
were,  therefore,  neglected.  Since  the  proprotors  are  inter¬ 
connected  and  indexed  such  that  only  symmetric  excitation  is 
possible,  the  asymmetric  mode  response  was  not  analyzed. 
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A  time  history  of  the  pylon  response  during  blade  feathering 
is  shown  in  Figure  81.  The  flight  condition  follows  the 
procedure  discussed  earlier;  airspeed  175  knots,  flapping 
freedom  locked  out. 

b.  Response  Amplitude 

Figures  82  and  83  show  the  calculated  vibratory  response 
during  blade  feathering  at  the  pylon  conversion  axis-shaft 
intersection  and  at  the  crew  station,  respectively,  as  a 
function  of  feathering  rate.  Also,  shown  in  Figure  83  is 
the  longitudinal  acceleration  caused  by  the  transient  thrust 
generated  by  blade  feathering.  This  in  combination  with  the 
vibratory  response  characteristics  suggests  that  a  feathering 
time  of  three  to  five  seconds  is  desirable. 

The  acceleration  levels  on  the  pylon,  at  the  engine  locations, 
and  on  the  fuselage,  are  all  well  below  normal  design  limits. 
The  stresses  in  the  wing  and  at  the  pylon  to  wing  attachment 
are  also  below  design  limits.  For  example,  the  peak  wing  tip 
deflection  (referenced  to  the  shaft  centerline)  is  ±1.6  inches. 
For  reference,  the  limit  beamwise  deflection  is  32  inches. 
Oscillatory  blade  loads  are  also  well  below  endurance  limits 
for  this  condition.  It  was  concluded  that  wing  and  pylon 
loads  and  vibration  are  not  a  limitation  on  blade  feathering. 

However,  the  crew  station  response  at  the  wing  beam  frequency, 
does  indicate  a  potential  human  factors  problem.  At  the  beam 
frequency  the  acceleration  level  is  close  to  the  level 
established  in  Reference  23  as  alarming.  The  transient  nature 
may  reduce  the  significance  of  this  observation,  but  a  require¬ 
ment  for  a  study  of  the  crew  acceptance  of  this  type  of 
response  is  indicated. 

The  response  sensitivity  to  airspeed  and  load  factor  was  also 
investigated  and  is  shown  in  Figures  84  and  85.  The  relative 
insensitivity  to  airspeed  is  explained  by  the  reduction  in 
wing  angle  of  attack  as  airspeed  is  increased.  Thus,  the 
wing  upwash  is  reduced  as  the  dynamic  pressure  increases. 

(In  making  the  airspeed  sensitivity  calculation,  it  was 
assumed  the  flaps  were  adjusted  to  maintain  a  shaft  angle  of 
attack  of  zero.) 

In  making  the  load  factor  calculation,  the  shaft  angle  of 
attack  was  increased  along  with  the  wing  angle  of  attack. 
Apparently  the  increased  shaft  angle  offsets  the  increase  in 
wing  upwash  at  increased  load  factor  to  keep  the  response 
constant  and  combines  with  reduced  upwash  at  reduced  load 
factors  to  decrease  the  response. 
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Load  factor  does  have  a  marked  effect  on  blade  loads.  This 
is  discussed  next. 

6.  Gust  and  Maneuver  Limits 


a .  Critical  Components 

The  blade  fold  hinge  loads  determine  the  gust  and  maneuver 
limits  during  feathering  and  folding  the  blades.  Drawing 
D270-410-003 ,  in  Appendix  II,  shows  the  fold  hinge  arrange¬ 
ment.  Tearout  of  the  spindle  folding  knuckle  lugs  is  the 
estimated  failure  mode.  The  spindle  is  titanium;  an  ultimate 
stress  of  130,000  was  assumed.  The  limit  bending  moment  is 
calculated  to  be  2.2  x  10 6  inch-pounds  (nonrotating).  Since 
the  folding  hinge  is  inboard  of  the  blade  spindle,  blade 
bending  moments  do  not  enter  the  control  system. 

The  wing  root  becomes  critical  when  the  blades  are  folded 
past  20-degree  fold  angle. 

b.  Methods 

A  discrete,  25-wing  chord  length,  1-cosine  shape,  vertical 
gust  was  used  in  the  analysis.  A  sharp  edged  gust  was  also 
investigated;  the  loads  imposed  were  about  15  percent  higher 
than  those  using  the  1-cosine  shaped  gust.  The  rigid  body 
vertical  motion  of  fuselage  was  included  in  the  analysis  to 
introduce  inertia  and  aerodynamic  relief.  However,  rigid 
body  pitching  was  neglected. 

Program  ARAP06  was  used  to  evaluate  the  limits  during  blade 
feathering.  For  gusts,  the  aircraft  was  assumed  to  be  in 
unaccelerated  flight.  The  gust  was  imposed  with  the  prop- 
rotor  at  constant  rpm,  with  the  rpm  varied  from  0  to  80 
percent  of  hover  rpm.  The  allowable  limit  load  was  then  used 
to  define  the  limit  gust.  For  maneuvers  the  wing  and  shaft 
angles  of  attack  were  adjusted  to  provide  the  required  lift. 

A  fatigue  oscillatory  load  limit  was  used  to  define  the 
maneuver  limit. 

The  limits  during  folding  were  determined  neglecting  stall 
and  are  thus  conservative.  The  blade  lift  curve  slope  was 
adjusted  with  fold  angle  in  accordance  with  the  lift  curve 
slope  variation  measured  in  the  dynamic  model  test.  The 
allowable  limit  loads  were  used  to  define  both  gust  and 
maneuver  limits. 

c .  Limits 

Figure  86  shows  the  gust  and  maneuver  limits  versus  rpm  and 
fold  angle.  The  influence  of  airspeed  on  the  limits  is  shown 
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in  Figure  87.  As  shown,  the  gust  limit  and  the  maneuver 
limit  during  blade  folding  are  acceptable.  A  small  increase 
in  fold  hinge  strength  would  result  in  meeting  the  required 
50  foot  per  second  gust  at  Vlr,  230  knots. 

Determination  of  the  maneuver  load  factor  limits  during 
feathering  raised  several  questions:  (1)  What  oscillatory 
load  limit  should  be  used?  And  (2),  at  what  rpm  should  the 
static  limit  be  used  rather  than  an  oscillatory  limit?  For 
example,  if  the  proprotor  stopping  time  is  five  seconds,  only 
ten  load  cycles  will  occur  during  feathering,  and  thus  a  low 
cyclic  limit  can  be  applied.  However,  when  flapping  is 
locked  out  prior  to  feathering,  a  large  number  of  cycles  can 
accumulate  so  the  endurance  limit  must  be  used  for  that 
condition.  And  at  some  low  rpm,  the  number  of  cycles  will  be 
very  low,  so  the  static  limit  may  be  used. 

For  this  study  the  following  procedure  was  followed:  (1) 
Between  80-  and  40-percent  hover  rpm,  the  endurance  limit 
(16,000  psi  for  titanium)  was  used.  The  limit  was  increased 
to  that  for  10,000  cycles  (30,000  psi)  at  20-percent  rpm,  and 
further  increased  to  the  static  limit  at  zero  rpm. 

The  resulting  maneuver  limits,  based  on  the  fold  hinge  lug, 
are  considered  too  low.  When  the  strength  of  the  fold  hinge 
lug  increased  to  that  of  the  blade  and  hub,  a  more  acceptable 
maneuver  limit  results.  Further  increases  (or  possibly  a 
decrease)  in  the  maneuver  limit  could  result  from  a  revised 
allowable  limit  schedule. 

7 .  Blade  Deflections 

Because  the  D270A  blades  and  hub  assembly  are  extremely  stiff 
relatively  small  blade  deflection  occurs  during  feathering  and 
folding.  This  is  significant  from  the  standpoint  of  wing-blade 
clearance.  For  reference,  the  calculated  peak  deflections 
under  extreme  conditions  during  folding  and  feathering  are 
tabulated  in  Table  VI.  The  inplane  deflections  are  greater 
than  the  out-of-plane  deflections  because  of  the  high  blade 
collective  pitch  during  folding  and  feathering. 


TABLE  VI 

BLADE  TIP  PEAK  DEFLECTIONS  -  AIRSPEED  175  KNOTS 


Condition 

Inplane 

(inches) 

Out-of-Plane 

(inches) 

Feathering, 

n  =  1 

5.5 

2.4 

Feathering, 

n  =  2 

20.1 

4.8 

Feathering, 

gust 

50-fps 

24.4 

4.8 
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TABLE  VI  -  Continued 


Condition  Inplane 

(inches) 

Folding,  n  =  1  2.6 

Folding,  n  =  2  9.6 

Folding,  50-fps  19.4 

gust 


B.  PARAMETRIC  STUDY 

The  sensitivity  of  the  dynamic  characteristics  during  blade 
feathering  and  folding  was  assessed  for  several  parameters. 

These  include: 

-  design  disc  loading 

-  proprotor  type 

-  wing  stiffness 

-  wing  to  proprotor  blade  spacing 

The  effect  of  feathering  rate  was  investigated  in  the 
correlation  study  and  the  dynamic  analysis  of  the  D270A. 

The  influence  of  detail  design  parameters  such  as  pylon 
center  of  gravity  and  pylon  incidence  was  also  investigated 
and  is  reported  in  this  section. 

1 .  Sensitivity  to  Design  Disc  Loading 

The  design  disc  loading  was  found  to  have  a  strong  influence 
on  proprotor  stability,  a  lesser  influence  on  dynamic  response 
during  feathering,  and  no  significant  influence  on  flutter  and 
divergence  during  folding,  or  on  the  gust  and  maneuver  limits. 

The  D270A  design  disc  loading  of  16.8  was  varied  by  ±25  per¬ 
cent  using  equal  mission  capability  as  a  design  constraint.  The 
design  value  of  16.8  was  selected  to  minimize  design  gross  weight, 
using  the  technique  described  in  Reference  24.  Table  VII  shows 
the  mass  and  stiffness  properties  for  the  disc  loadings  investi¬ 
gated.  These  were  estimated  using  the  methods  outlined  in  Ref¬ 
erence  1G.  It  was  assumed  that  design  factors  such  as  placement 
of  the  blade  natural  frequencies  were  the  same  regardless  of 
disc  loading. 

Figure  88  shows  the  influence  of  disc  loading  on  proprotor/ 
pylon  stability.  Note  that  the  influence  is  strongest  in 
proprotor  mode,  when  the  hub  gimbal  is  free  to  flap.  However, 
a  significant  influence  is  also  evident  when  the  gimbal  is 

157 


Out-of-Plane 

(inches) 

0.2 

0.52 

7.6 


locked  out.  The  high  sensitivity  to  disc  loading  is  explained 
by  the  variation  in  blade  Lock  number.  The  D270A  blade  Lock 
number  is  6.2;  for  a  disc  loading  of  12  the  Lock  number 
increases  to  8.7  and  for  20  it  decreases  to  5.0.  Since  the 
Lock  number  represents  the  ratio  of  aerodynamic  and  inertia 
forces  acting  on  the  blades  and  consequently  the  wing/pylon, 
a  higher  Lock  number  generally  implies  reduced  stability. 

The  variation  in  response  during  feathering  is  shown  in 
Figure  89.  Note  that  lower  disc  loading  leads  to  a  higher 
response.  The  reason  for  this  is  partially  explained  by  the 
higher  Lock  number,  and  partially  by  the  lower  wing  loading 
at  the  same  lower  disc  loading. 

The  reason  for  the  negligible  influence  of  disc  loading  on 
flutter  and  divergence  and  gust  and  maneuver  limits  is  the 
wing  and  blade  strength  to  airload  ratio  remains  a  design 
constant  regardless  of  disc  loading. 

2 .  Proprotor  Type 

There  are  several  configurations  possible  for  a  folding  prop- 
rotor.  The  gimbaled,  semirigid  configuration  such  as  used 
on  the  D270A  is  one  configuration.  Other  possible  configua- 
tions  are  the  hingeless  or  "rigid"  rotor  which  may  be  either 
soft  or  stiff  in  the  plane  of  the  rotor,  and  the  fully 
articulated  and  free-to-flap  configurations. 

Analysis  of  all  of  these  configurations  was  not  practical. 

The  hingeless,  soft-inplane  type  was  chosen  for  analysis  since 
it  represents  an  extreme  in  flexibility  compared  to  the 
gimbaled  semirigid  type.  It  should  be  noted  that  locking  out 
the  flapping  freedom  before  feathering  on  the  fully  articulated 
(and  lag  freedom)  and  the  free-to-flap  configurations  appears 
necessary  to  prevent  excessive  blade  motions.  When  this  is 
done  these  rotor  configurations  approach  either  the  gimbaled- 
semirigid  or  the  hingeless-soft  inplane  configurations  in 
dynamic  characteristics,  depending  on  the  blade  flexibility. 

a.  Hingeless-Soft  Inplane  Proprotor  Parameters 

The  para  ters  of  a  49.2-foot  diameter,  four-bladed,  hinge- 
less  folding  proprotor,  developed  under  contract  to  the  Air 
Force  we^*  used  in  this  study.  Figure  90  shows  the  blade 
stiffne**  mass  distributions  of  the  proprotor.  The 
blade  na  at  frequencies  are  shown  in  Figure  91. 

In  comparison  to  the  D270  proprotor  blades  (see  Appendix  II)  , 
this  proprotor  may  be  characterized  as  extremely  flexible. 

The  blade  weight  is  approximately  55  percent  of  the  D270 
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blade  weight,  and  the  stiffness  at  10  percent  radius  only 
about  7.8  percent  of  that  of  the  D270A  blade  in  the  feather¬ 
ing  spindle  region,  lh  percent  radius. 

For  this  study  a  three-bladed  hingeless  soft-inplane  prop- 
rotor  was  substituted  for  the  D270A  gimballed  proprotor. 

The  wing/pylon  dynamic  characteristics  were  not  modified. 


TABLE  VII 

WEIGHT  AND  STIFFNESS  VERSUS  DISC  LOADING* 

DL  =  12.0  DL  =  16.8  DL  ■  20.0 


Gross  Weight 

lb 

82500 

66000 

69300 

Weight  empty 

lb 

55000 

44900 

45200 

Rotor  radius 

ft 

33.0 

25.0 

23.5 

Blade  chord 

in . 

30.3 

40.0 

42.5 

Wing  span 

ft 

80.25 

64.25 

61.25 

Wing  area 

ft2 

1101 

706 

641 

Wing  loading 

lb/ft2 

74.9 

93.5 

108.1 

Wing  root 
stiffness 

lb-in2 

133  x  109 

68  x  109 

65  x 

Wing  tip 
stiffness 

lb-in2 

15.6  x  109 

10  x  109 

10  x 

Rotor  weight 

lb 

9005 

7521 

7819 

Pylon  weight 

lb 

4601 

3849 

4002 

Wing  weight 

lb 

8016 

5137 

4963 

*Basis : 

Equal  mission  capability 

b.  Dynamic  Stability  Characteristics 

The  proprotor/pylon  stability  boundary  with  the  hingeless- 
soft-inplane  proprotor  is  shown  in  Figure  92  as  a  function 
of  rpm.  In  the  proprotor  cruise  mode  rpm  range  the  mode  of 
instability  is  the  same  as  that  with  the  gimbaled  proprotor, 
the  wing  chordwise  bending  mode.  The  stability  boundary  is 
slightly  higher,  however.  Below  150  rpm  the  mode  of  in¬ 
stability  changes  to  one  involving  coupled  flapping,  blade 
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inplane,  and  wing  beam  motion.  The  boundary  decreases  with 
rpm,  and  is  a  minimum  near  zero  rpm. 

The  low  rpm  instability  appears  to  be  of  the  same  nature  as 
that  measured  with  the  dynamic  model  when  the  gimbal  was  free 
to  flap;  coalescence  of  the  blade  flapping  and  wing  beam 
modes.  A  plot  of  flapping  natural  frequency  versus  wing  beam 
natural  frequency  indicates  the  frequencies  would  coalesce  at 
about  80  rpm. 

With  the  blades  stopped,  blade  flutter  and  divergence  are 
more  of  a  problem  than  with  the  D270A  proprotor.  The  blade 
divergence  boundary  is  greatly  reduced.  Figure  93  shows  the 
variation  in  blade  beam  bending  natural  frequency  with  air¬ 
speed  for  sideslip  angles  of  5,  10,  and  20  degrees.  This 
characteristic  is  similar  to  that  measured  with  hingeless 
stiff  inplane  edgewise  stoppable  rotors. 24 

Since  the  divergence  dynamic  pressure  is  essentially  propor¬ 
tional  to  the  blade  stiffness,  it  appears  divergence  is  a 
fundamental  problem  associated  with  a  hingeless  type  folding 
proprotor. 

A  second  aeroelastic  characteristic  during  folding,  is  also 
a  potential  problem  for  the  gimballed  type  folding  proprotor, 
a  coupled  blade  beam  bending  wing-chord  bending  flutter. 
Figures  94  and  95  show  the  calculated  variation  of  frequency 
and  damping  with  airspeed  for  a  60-degree  fold  angle.  As  was 
observed  with  the  gimballed  rotor,  the  blade  bending  modes 
increase  in  frequency  and  transit  the  wing  beam  and  chord 
bending,  and  torsion  modes.  At  each  crossing  the  damping  of 
the  wing  modes  decreased.  For  the  beam  and  torsion  modes  the 
reduction  in  damping  does  not  overcome  the  available  aero¬ 
dynamic  damping.  But  for  the  wing  chord  mode  negative  damping 
is  predicted.  Whether  or  not  this  would  result  in  flutter 
depends  on  the  structural  damping. 

c.  Response  Characteristics 

The  pylon  vibration  levels  during  blade  feathering  are 
similar  to  those  with  the  gimballed  proprotor,  but  are  more 
sensitive  to  load  factor.  Figure  96  compares  the  vibration 
levels  calculated  for  both  types.  The  increased  sensitivity 
of  the  hingeless  proprotor  is  caused  by  the  larger  blade  tip 
deflections  associated  with  the  blade  flexibility.  Figure  97 
compares  the  blade  tip  peak  out-of-plane  deflections  during 
feathering  for  the  two  proprotor  types  as  a  function  of  load 
fcictors.  Blade  tip  deflections  toward  the  wing  leading  edge 
represent  an  effective  decrease  in  wing-rotor  spacing,  which 
is  shown  later  in  this  section  to  increase  response. 


160 


Blade  oscillatory  deflections  during  feathering  were  also 
investigated.  Since  the  fundamental  inplane  natural  frequency 
is  below  operating  rpm  in  the  normal  operating  rpm  range,  about 
0.7-per-rev,  1-per-rev  resonance  is  transmitted  during  feather¬ 
ing.  For  very  low  rates  of  feathering,  1-per-rev  resonance  of 
the  first  inplane  mode  could  be  a  load  problem.  However  for 
rapid  rates  the  resonant  build  up  is  relatively  small.  Figure 
98  shows  the  blade  tip  peak  inplane  deflections  versus  rpm 
during  a  four  second  blade  feathering.  The  oscillation  limit 
for  the  rotor  is  not  known  but  the  magnitude  shown  is  clearly 
within  the  capability  of  such  a  flexible  blade. 

d.  Gust  and  Maneuver  Limit 


Since  the  strength  characteristics  of  the  subject  proprotor 
are  not  known,  it  is  not  possible  to  accurately  assess  a  change 
in  the  gust  and  maneuver  limits.  The  calculated  blade  tip 
deflections  for  several  gust  and  maneuver  conditions  are  given 
in  Table  VIII  and  are  large  compared  to  those  for  the  gimbaled 
semirigid  type  but  appear  to  be  within  the  strength  capability 
of  the  hingeless  blade.  The  associated  fold  hinge  loads  are 
somewhat  lower,  due  to  the  relieving  effect  of  the  blade 
motion.  Thus  slightly  higher  gust  and  maneuver  limits  (during 
folding)  could  result  with  the  hingeless  proprotor. 

An  increase  in  the  limit  maneuver  load  factor  during  feather¬ 
ing  might  appear  to  be  a  benefit  of  using  the  hingeless  prop¬ 
rotor  type  since  flapping  freedom  is  retained  during  feather¬ 
ing.  However,  the  very  large  deflections  at  low  rpm  indicate 
that  high,  low  cycle,  blade  loads  may  result  in  maneuver  limit 
during  feathering,  in  the  same  range  as  those  of  the  stiffer, 
gimbaled  semirigid  proprotor.  Figure  99  shows  the  blade  tip 
peak  deflection  envelope  versus  rpm  for  a  load  factor  of  1.3 
to  illustrate  this  consideration. 

TABLE  VIII 


BLADE  TIP  PEAK  DEFLECTIONS 
HINGELESS,  SOFT-INPLANE  PROPROTOR 
AIRSPEED  175  KNOTS 


Condition 

Inplane 

Out-of-Plane 

Feathering, 

n  s 

=  1 

15.5 

9.8 

Feathering , 

n  * 

=  1.5 

65.0 

22.5 

Feathering , 

50- 

-fps 

83.6 

26.2 

gust 
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TABLE  VIII  -  Continued 


Condition 

Inplane 

Out-of-Plane 

Folding,  n  =  1 

2.9 

0.91 

Folding,  n  *  1.5 

57.4 

18.0 

Folding,  50-fps 
gust 

153.0 

47.8 

3 .  Wing  Stiffness 

Wing  stiffness  has  a  controlling  influence  on  proprotor 
stability  and  flutter  and  divergence.  Wing  stiffness  can  be 
used  to  a  lesser  extent  to  control  dynamic  response  during 
blade  feathering.  The  influence  on  gust  and  maneuver  limits 
is  small,  since  in  general  the  blade  loads  are  more  critical. 

a.  Influence  on  Stability 

Figure  100  shows  the  effect  of  varying  the  wing  beam,  chord, 
and  torsional  stiffness  (reference  to  the  wing  tip  spring 
rates)  separately  and  simultaneously,  on  stability  in  prop- 
rotor  mode  (flapping  free) .  Wing  torsion  has  the  strongest 
influence,  up  to  the  design  stiffness.  Above  that  stiffness, 
it  is  ineffective  since  the  mode  of  instability  becomes  that 
of  the  wing  chord  mode.  To  effect  a  further  increase  in 
stability,  the  wing  torsion  and  chord  stiffness  would  have  to 
be  increased  simultaneously.  Note  that  wing  beam  stiffness 
does  not  have  a  strong  influence. 

When  flapping  is  locked  out  the  proprotor  stability  becomes  a 
propeller  whirl  flutter  problem  and  relative  isotropicity  of 
the  wing/pylon  modes  becomes  a  controlling  factor. 20  Figure 
101  shows  the  effect  of  varying  the  wing  stiffness  when 
flapping  is  locked  out.  Since  the  fundamental  beam  and 
chordwise  modes  are  close  together  in  frequency,  varying  the  beam 
and  chord  stiffness  so  as  to  reduce  the  isotropicity  increas¬ 
ing  the  stability.  Increasing  the  wing  torsional  stiffness 
increases  stability  since  it  reduces  the  amount  of  shaft 
pitching  in  the  beamwise  mode.  When  the  torsion  stiffness  is 
increased  more  than  25  percent,  the  mode  of  instability  changes 
to  a  forward  whirl. 

The  influence  of  wing  stiffness  on  the  flutter  and  divergence 
characteristics  during  folding  was  found  to  be  straightforward. 
The  divergence  speed  varied  proportionally  to  the  square  root  of 
the  wing  torsional  stiffness,  and  was  not  significantly 
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influenced  by  the  wing  beamwise  or  chordwise  stiffness.  The 
flutter  velocity  was  not  sensitive  to  nominal  variation  in 
the  wing  stiffness  properties. 

b.  Influence  on  Response 

The  influence  of  wing  stiffness  on  the  dynamic  response  is 
somewhat  difficult  to  assess.  Figure  102  shows  the  influence 
of  wing  beamwise  and  torsional  stiffness  on  the  dynamic 
response.  Increasing  the  wing  torsional  stiffness  has  the 
opposite  effect  of  what  one  might  expect  on  the  basis  of 
superficial  considerations.  However,  the  reason  for  the 
increased  response  with  increased  torsional  stiffness  is 
logical  if  the  effect  of  torsional  stiffness  on  the  system 
damping  is  considered.  For  example,  for  the  wing  beam  mode 
an  increase  in  torsional  stiffness  decreases  the  amount  of 
coupled  pylon  pitching  motion.  This  decrease  in  pitching 
decreases  the  proprotor  aerodynamic  damping  of  the  wing  beam 
mode  (Onq  contribution) .  Since  the  forced  response  is 
essentially  proportional  to  the  system  damping  the  response 
is  increased. 

Varying  the  wing  beam  stiffness  has  a  smaller  influence  since 
the  system  damping  is  not  significantly  affected  by  beam 
stiffness.  However,  it  should  be  pointed  out  that  increasing 
the  beam  stiffness,  and  thus  the  beam  frequency  may  be  desir¬ 
able  from  a  human  factor  standpoint,  since  for  a  given  vibra¬ 
tional  level  the  displacement  is  reduced  in  proportion  to  the 
frequency  squared;  i.e.,  directly  proportional  to  the  stiff¬ 
ness.  Furthermore,  the  human  tolerance  to  vibration  is 
frequency  dependent  and  adjustment  of  the  beam  stiffness  could 
be  used  to  tune  the  wing  beam  mode  to  a  more  "comfortable" 
frequency. 

c .  Wing  Thickness  Influence 

Increasing  the  wing  thickness  is  the  most  efficient  means  of 
increasing  wing  stiffness  for  a  given  material/structural  con¬ 
cept.  The  influence  of  wing  thickness  on  the  response  during 
feathering  is  shown  in  Figure  103.  The  relatively  small 
increase  in  response  with  increased  thickness  suggests  that 
the  wing  circulation  is  the  strongest  source  of  excitation. 

4 .  Wing  to  Proprotor  Blade  Spacing 

Wing  to  rotor  spacing  has  a  strong  influence  on  proprotor 
stability,  divergence  and  dynamic  response.  The  influence 
on  gust  and  maneuver  limits  is  small  since  the  blade  is 
critical.  There  are  two  means  of  varying  the  spacing;  wing 
sweep  and  mast  length. 
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a.  Influence  on  Stability 

The  influence  of  mast  length  on  stability  is  indicated  in 
Figure  104.  The  influence  of  forward  sweep  on  stability,  for 
sweep  angles  of  up  to  12  degrees  (twice  that  of  the  D270A's 
forward  sweep  of  6  degrees)  was  found  to  be  negligible.  This 
is  because  the  wing  beam  stiffness  is  very  high  and  thus  does 
not  impact  on  the  effective  torsional  stiffness  as  strongly  as 
for  conventional  aircraft  designs. 

With  flapping  free,  for  proprotor  mode,  increasing  the  mast 
stiffness  decreases  the  proprotor  stability  and  would  require 
an  increase  in  wing  stiffness  to  maintain  the  stability  margin. 
With  flapping  locked  out,  however,  the  stability  increases  up 
to  a  point,  with  increasing  mast  length.  This  increase  in 
stability  with  mast  length  follows  the  travel  for  propellers 
noted  in  Reference  20.  Further  increases  in  mast  length  cause 
a  reduction  in  the  wing  divergence  speed. 

b.  Influence  on  Response 

Increasing  the  wing  to  rotor  spacing  reduces  the  dynamic 
response  during  feathering.  Figure  105  shows  the  effect  of 
these  parameters  on  the  response  at  the  wing  beam  frequency. 

It  should  be  noted  at  this  point  that  the  wing  to  rotor  spacing 
has  a  strong  influence  on  the  aircraft  stability  and  control 
characteristics  and  may  not  be  readily  used  to  control  dynamic 
characteristics . 

5 .  Other  Parameters  Investigated 

Several  secondary  design  parameters  were  also  varied  to 
investigate  their  influence. 

a .  Pylon  Incidence 

The  influence  of  load  factor  on  dynamic  response  suggested 
that  pylon  incidence  might  be  a  means  of  reducing  dynamic 
response.  Figure  106  shows  that  ±5  degrees  of  incidence 
reduces  the  response  by  about  25  percent.  However  this 
approach  to  reducing  response  is  limited  by  blade  load  con¬ 
siderations  . 

b.  Wing  Structural  Damping 

The  wing  structural  damping  was  assumed  to  be  one  percent  of 
critical  throughout  the  dynamic  analysis.  An  increase 
in  structural  damping  (or  the  application  of  external  damp¬ 
ing)  was  investigated  to  determine  if  that  could  be  used  to 
reduce  dynamic  response.  The  results,  shown  in  Figure  107 
indicate  that  a  resonable  increase  in  damping  is  ineffective. 
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c.  Pylon  Center  of  Gravity 

Since  the  location  of  the  pylon  center  of  gravity  is  con¬ 
trollable  to  some  extent,  an  investigation  was  made  to 
determine  its  influence.  This  investigation  showed  the  pylon 
center  of  gravity  location  has  a  negligible  influence  on 
stability  and  response. 

d .  Blade  Fold-Hinge  Arrangement 

As  noted  earlier,  the  D270A  fold  hinge  is  located  inboard  of 
the  pitch  change  spindle  to  prevent  the  blade  airloads  from 
acting  through  the  control  system  during  blade  folding. 

If  the  fold  hinge  were  located  outboard  of  the  pitch  change 
spindle,  the  gust  limit  (and  maneuver  load  factor  limit) 
would  be  greatly  reduced.  Figure  108  compares  the  limit 
gust  for  the  alternate  configuration  with  that  of  the  D270A. 

To  raise  the  allowable  gust  limit  to  the  level  of  the  D270A, 
a  feathering  lockout  would  be  required  to  prevent  the  blade 
airloads  from  acting  on  the  control  system. 
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Figure  68.  Airframe  Symmetric  Natural  Frequencies 
Versus  Pylon  Conversion  Angle. 
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NATURAL  FREQUENCY  -  CPS 
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Figure  69.  Airframe  Asymmetric  Natural  Frequencies  Versus 
Pylon  Conversion  Angle. 
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ure  70.  Proprotor  Collective  Natural 

Frequencies,  Girabal  Free  to  Flap 
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Figure  71.  Proprotor  Cyclic  Natural 

Frequencies,  Gimbal  Free-To 


n 

- 

r.;  i 

■t  •  •  *  y*  *  ^  ^  ** 

. 

.  ,  .  . 

+£f 

r—  ». 

h  ii 

PROPROTOR  RPM 


Figure  73.  Proprotor  Stability  Boundary,  Proprotor  Cruise 
Mode? . 
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Figure  74.  Root  Motion  With  Airspeed,  Proprotor 
Mode . 
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Figure  75.  Proprotor  Stability  Boundary  During  Blade 
Feathering. 
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Figure  77.  Flutter  and  Divergence  Boundary  During 
Blade  Folding. 
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Figure  79.  Frequency  Variation  with  Airspeed, 
60  Degrees  Blade  Fold. 
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Figure  78.  Frequency  Variation  with  Airspeed, 
Zero  Degree  Blade  Fold. 
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DAMPING 


Figure  80.  Damping  Variations  with  Airspeed,  60-Degrees 
Blade  Fold. 
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Figure  81. 
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Figure  82.  Peak  Vibration  at  Pylon  Conversion  Axis 
During  Blade  Feathering. 
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CREW  STATION  VIBRATION 
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Figure  83.  Peak  Vibration  at  Crew  Station  During  Blade 
Feathering . 
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Figure  84.  Influence  of  Airspeed  on  Crew  Station 
Vibration. 
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CREW  STATION  VIBRATION 
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LIMIT  VERTICAL  GUST  VELOCITY  -  FT/SEC 
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Figure  86.  Gust  and  Maneuver  Limits  During  Feathering 
and  Folding  at  175  Knots  Airspeed. 
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Figure  87.  Influence  of  Airspeed  on  Gust  and  Maneuver 
Limits  During  Feathering  and  Folding. 
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Figure  89.  Influence  of  Disc  Loading  on 

PyLon  Vibration  During  Feathering. 
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Figure  90.  Mass  and  Stiffness  Parameters 
for  HingeLess,  Soft-InpLane 
Folding  Proprotor. 
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Figure  91.  Hingeless  Proprotor  Natural 
Frequency  Spectrum  During 
Blade  Feathering. 
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Figure  93.  Blade  Divergence  Airspeed  Versus  Sideslip 
Angle,  Hingeless  Soft- Inplane  Proprotor, 
Zero  Degrees  Fold. 
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Figure  94.  Natural  Frequency  Variation  With  Airspeed 
During  Blade  Folding,  Hingeless  Soft- 
Inplane  Proprotor,  60  Degrees  Fold. 
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Figure  95.  Damping  Variation  with  Airspeed  During 
Blade  Folding,  Hingeless  Soft  Airplane 
Proprotor,  60  Degrees  Fold. 
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VERTICAL  VIBRATION  AT  CREW  STATION 
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Figure  96.  Comparison  of  Crew  Station  Vibration 
During  Blade  Feathering. 
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Figure  97.  Comparison  of  Out-of -Plane  Tip  Deflections, 
Hingeless  Soft-Inplane  and  Gimbaled,  Semi- 
Rigid,  Proprotors. 
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Figure  98.  Envelope  of  Blade  Inplane  Deflections  During 
Blade  Feathering,  Hingeless  Soft-In-Plane 
Proprotor. 
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Figure  99.  Influence  of  Load  Factor  on  Blade  Inplane 
Deflections,  Hingeless,  Soft-Inplane, 
Proprotor. 
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PERCENT  DESIGN  WING  STIFFNESS 


Figure  100.  Influence  of  Wing  Stiffness  on 
Proprotor  Stability,  Gimbal 
Free  to  Flap. 
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PERCENT  DESIGN  WING  STIFFNESS 
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Figure  I0Io  Influence  of  Wing  Stiffness  on  Proprotor 
Stability,  Gimbal  Flapping  Freedom 
Locked  Out. 
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PYLON  CONVERSION  AXIS  VIBRATION 


Figure  102,  Influence  of  Wing  Stiffness  on 
Vibration  During  Feathering. 
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Figure  103.  Influence  of  Wing  Thickness  on 
Vibration  During  Feathering. 
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PYLON  CONVERSION  AXIS  VIBRATION 


WING  FORWARD  SWEEP  ANGLE  -  DEG 

Figure  L05.  Influence  of  Wing  to  Proprotor  Spacing  and 
Wing  Forward  Sweep  on  Vibration  During 
Blade  Feathering. 
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Figure  106.  Influence  of  Pylon  Incidence  on 
Vibration  During  Feathering. 
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Figure  107.  Influence  of  Wing  Structural  Damping 
on  Vibration  During  Feathering. 
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LIMIT  VERTICAL  GUST  VELOCITY  -  FT/SEC 
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Figure  108.  Influence  of  Fold  Hinge  Location  on 
Gust  Limit  During  Blade  Folding. 
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SECTION  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

A.  THEORY  AND  COMPUTER  PROGRAMS 

A  theory  has  been  developed  to  predict  the  stability  and 
response  during  feathering  and  folding  the  blades  of  a  fold¬ 
ing  proprotor  VTOL.  The  theory,  manifested  as  digital  computer 
programs,  is  in  good  agreement  with  measured  stability  and 
response  characteristics. 

The  computer  programs  are  presently  restricted  to  some  extent 
in  the  configurations  they  are  capable  of  simulating.  For 
example,  the  proprotor  simulation  is  restricted  to  three 
blades.  These  restrictions  arose  from  computer  storage 
limitations  and  may  be  relaxed  if  more  storage  is  available. 
Information  is  provided  in  this  volume  and  in  Volume  II  to 
assist  a  user  in  modifying  and  extending  the  programs. 

It  should  be  noted  that  the  Feathering-Dynamics  Analysis 
Program  ARAP06 ,  is  also  applicable  to  tilt-proprotor  dynamic 
analysis .  Capability  in  this  regard  includes  prediction  of 
proprotor  stability,  blade  loads,  and  vibration  character¬ 
istics.  The  program  may  also  be  used  to  investigate  the 
dynamic  characteristics  in  helicopter  and  conversion  mode  for 
both  proprotor  VTOL  types. 

B.  CONCLUSIONS  FROM  THE  DYNAMIC  MODEL  TEST 

Several  conclusions  regarding  the  folding  proprotor  configura¬ 
tion  can  be  made  directly  from  the  dynamic  model  test. 

The  most  significant  conclusions  is  that  flapping  freedom  must 
be  locked  out  prior  to  feathering  the  blades.  During  tests  with 
the  gimbal  flapping  only  lightly  restrained,  large  flapping 
amplitudes  were  encountered  even  near  zero  shaft  angle  of  attack. 
As  the  blade  feathering  rate  was  increased,  the  peak  flapping 
excursion  increased.  At  a  shaft  angle  of  attack  of  5.6  degrees 
the  blades  contacted  the  mechanical  flapping  stopt,  of  ±12 
degrees.  Furthermore  a  flapping  instability  was  encountered 
at  low  rpm,  in  the  transition  mode  airspeed  range. 

When  the  gimbal  flapping  freedom  was  locked  out  the  blade  tip 
deflections  were  observed  to  be  small  and  the  flapping 
instability  eliminated. 

Locking  out  the  blade  flapping  freedom  appears  necessary 
regardless  of  how  the  freedom  is  attained,  i.e.,  with  a  gimbal, 
a  true  flapping  hinge,  or  using  an  elastic  virtual  hinge, 
as  with  a  hingeless  rotor. 
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C.  CONCLUSIONS  FROM  THE  DYNAMIC  ANALYSIS  AND  PARAMETRIC 
STUDY 

The  dynamic  analysis  indicates  that  acceptable  dynamic  stabil¬ 
ity  and  response  characteristics  during  blade  feathering  and 
folding  can  be  obtained  in  a  66,000-pound  gross  weight  FPR 
design.  The  most  serious  problem  area  identified  in  the 
course  of  the  analysis,  is  a  low  maneuver  load  factor  capa¬ 
bility  during  blade  feathering.  This  arises  from  the  require¬ 
ment  to  use  fatigue  allowable  load  limits  during  feathering 
rather  than  static  limit  loads.  For  the  D270A,  another  factor 
is  the  strength  of  the  critical  component,  the  fold  hinge  lug. 
When  the  lug  strength  is  increased  to  the  same  level  as  that 
of  the  hub  and  blades,  the  maneuver  load  limit  is  considerably 
higher. 

It  is  apparent  from  this  study  that  a  criterion  for  allowable 
loads  during  blade  feathering  must  be  established. 

A  second  problem  area  is  the  acceptability  of  the  crew  station 
vibration  during  blade  feathering.  The  largest  response  occurs 
at  a  frequency  of  approximately  two  cycles  per  second,  the 
wing  beam  bending  natural  frequency.  At  this  frequency  the 
peak  accelerations  reach  ±0.65  g,  which  exceeds  an  "alarming" 
subjective  discomfort  level. 2 3  However,  the  reference 
subject  discomfort  level  was  established  for  sinusoidal  vibra¬ 
tion  whereas  the  response  during  feathering  is  a  damped 
transient.  A  requirement  for  a  criteria  for  an  acceptable 
level  of  crew  station  response  is  evident. 

The  results  of  the  parametric  study  are  summarized  in  terms  of 
design  guidelines,  Table  IX.  A  brief  discussion  of  the  influ¬ 
ence  of  major  design  parameters  is  given  below. 

(1)  Disc  Loading  -  A  higher  disc  loading  has  a  favorable 
influence  on  proprotor  stability  and  on  dynamic 
response.  The  influence  of  disc  loading  on  other 
dynamic  factors  is  small  provided  the  system 
natural  frequency  ratios  are  held  constant. 

(2)  Proprotor  Types  -  A  proprotor  type  which  combines 
stiff  blades  with  the  capability  to  lock  out  flapping 
freedom,  whether  real  or  virtual,  is  desirable. 

Highly  flexible,  or  freely  flapping  or  lightly 
restrained  in  flapping  blades  are  susceptible  to 
instability  at  low  rotor  speeds  and  can  contribute  to 
higher  response  levels  during  feathering. 

(3)  Wing  Stiffness  -  The  wing  chordwise  and  torsional 
stiffness  have  a  strong  influence  on  proprotor/ 
pylon  stability  when  the  blades  are  free  to  flap. 
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Increasing  the  torsional  and  chordwise  stiffness 
increases  the  stability.  The  wing  teamwise  stiff¬ 
ness  has  a  small  influence.  When  flapping  freedom 
is  locked,  the  beamwise  stiffness  is  also  significant. 
However,  with  flapping  locked  out  the  important 
consideration  is  frequency  isotropicity .  Conse¬ 
quently  a  decrease  in  stiffness  can  be  beneficial 
if  the  system  frequency  anisotropy  is  increased. 

The  influence  of  wing  stiffness  on  the  dynamic 
response  is  small  for  a  reasonable  variation  in 
stiffness.  However,  a  reasonable  change  in  stiff¬ 
ness  can  have  a  favorable  influence  if  the  system 
frequencies  are  shifted  toward  a  more  "comfortable" 
range,  from  the  human  factors  standpoint. 

Wing  thickness  ratio,  the  usual  means  of  control¬ 
ling  wing  stiffness  was  found  to  have  a  small 
influence  on  dynamic  response. 

(4)  Wing-Rotor  Spacing  -  Increased  spacing,  obtained 
using  a  longer  mast,  has  a  destabilizing  influence 
when  the  blades  are  free  to  flap,  but  is  stabilizing 
when  flapping  is  locked  out.  A  longer  mast  decreases 
the  dynamic  response . 

Wing  forward  sweep  may  be  used  to  increase  the 
effective  spacing  and  thus  reduce  the  response, 
and  does  not  significantly  influence  stability. 

(5)  Feathering  Rate  -  Based  on  response  studies  for 
25-foot,  38.5-foot,  and  the  D270A  50-foot  proprotor, 
an  optimum  feathering  time  may  be  approximated  by 

J  D 

AT  =  2.75J38.5  seconds  (225) 

where  D  is  the  proprotor  diameter  in  feet. 

This  optimum  is  not  extremely  sensitive  and  may 
be  varied  over  a  considerable  range. 

D.  DESIGN  CRITERIA  IMPLICATIONS 

The  flutter-free  margin  specified  in  MIL-A-8870  and  the 
50-foot  per  second  vertical  gust  encounter  at  Vlf,  specified 
in  MIL-A-8860  appear  to  be  within  the  capability  of  the 
folding  proprotor  concept.  However,  the  2.0  maneuver  load 
factor  at  V^p,  specified  in  MIL-A-8860,  may  be  difficult  to 
meet  if  fatigue  allowable  load  limits  are  applied  during 
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blade  feathering.  It  is  apparent  that  a  criterion  for  allow¬ 
able  loads  during  blade  feathering  is  required. 

As  noted  above,  design  criteria  for  acceptable  crew  station 
acceleration  levels,  and  allowable  loads  during  blade 
feathering  should  be  developed. 

E.  RECOMMENDATIONS  FOR  FURTHER  STUDY 

Experimental  investigations  should  be  extended  to  include 
factors  such  as  wing  flap  and  free-f light  conditions,  which 
may  have  a  significant  impact  on  the  response  during  feathering. 
The  present  theory  can  be  readily  extended  to  include  these 
details . 

Full-scale  tests  are  required  to  identify  the  hardware 
dynamic  problems  associated  with  feathering  and  folding. 

For  example,  the  flapping  lockout  and  blade  folding  mechanism 
could  induce  high  dynamic  loads  depending  on  their 
characteristics . 

Much  of  this  recommended  work  is  presently  being  accomplished 
under  Air  Force  and  NASA  sponsored  research. 
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(1)  Generally  disc  loading  will  be  established  by  a 
gross  weight  consideration  rather  than  dynamics 
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APPENDIX  I 


MODEL  DESCRIPTIVE  DATA 


Transition  Model 


The  experimental  data  on  wing-blade  aerodynamic  interference 
presented  in  this  report  was  obtained  using  the  transition 
model  shown  earlier  in  Figure  24.  The  scaling  of  the  transition 
model  is  the  same  as  the  semispan  aeroelastic  model  discussed 
later  in  this  Appendix.  The  fuselage  is  free  to  pitch  on  the 
wind  tunnel  balance  with  the  pivot  located  at  35  percent  MAC. 

The  fuselage  freedom  may  be  locked  out,  if  desired.  The  pylon 
conversion  angle  relative  to  the  wing  chordline,  and  the 
horizonal  stabilizer  incidence  can  be  remotely  controlled.  A 
three  component  strain  gage  balance  is  used  to  measure  the 
proprotor  hub  normal  forces  and  sideforce  and  pitch  and  yaw 
moments.  A  strain  gage  at  the  blade  root  may  be  used  to  measure 
the  blade  beam  bending  moment. 

Figure  1-1  shows  the  transition  model  layout  and  dimensions. 

The  wing  thickness  ratio  is  0.15.  Two  sets  of  blades  are 
available.  One  set  has  -20  degrees  twist,  the  other  set,  -3 
degrees  twist.  (This  is  the  theoretical  twist,  shaft  center- 
line  to  blade  tip.)  Tests  of  this  model  were  made  in  1968 
under  the  Bell  IR&D  program. 

Semispan  Aeroelastic  Model 

The  model  used  in  the  wind  tunnel  test  conducted  under  this 
program  is  a  0.133  scale,  dynamic  aeroelastic  model  of  a  Bell 
Helicopter  Company  38,000-pound  tilt-proprotor  VTOL  design. 

The  model  was  designed,  fabricated,  and  tested  by  Bell  under 
the  Army  Composite  Aircraft  Program.  At  the  conclusion  of 
that  program  the  model  was  transferred  by  the  Army  to  the 
NASA-Langley  Aeroelasticity  Branch  for  further  research  in 
proprotor  dynamics.  For  a  description  of  the  full-scale  design 
the  reader  is  referred  to  Reference  25. 

A  photograph  of  the  model  installed  in  the  16-foot  transonic 
dynamics  tunnel  was  shown  earlier  in  Figure  29.  The  proprotor. 
pylon,  and  wing  are  Froude  and  Lock  scaled  for  operation  in  air 
at  full-scale  density.  The  model  scale  factors,  based  on  the 
tilt-proprotor  design  are  given  in  Table  I-I. 

The  fuselage  and  mounting  arrangement  may  be  taken  as  effectively 
rigid  with  respect  to  the  proprotor/pylon/wing  system.  While 
the  mounting  is  not  that  normally  employed  for  semispan  models, 
tests  with  a  conventional  reflection  plane  showed  no  significant 
difference  in  stability  or  dynamic  response. 
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MODEL 

TABLE  I- I 

SCALE  FACTORS 

Parameter 

Units 

Scale  Factor 

Length 

ft* 

0.133 

Mass 

lb 

0.00235 

Time 

sec 

0.365 

Froude  Number 

— 

1.0 

Lock  Number 

— 

1.0 

Mach  Number 

— 

0.365 

Reynolds  Number 

— 

0.0487 

Velocity 

ft/sec 

0.365 

Frequency 

1/sec 

2.738 

Inertia 

lb»ft2 

4.21  x  10“5 

Linear  Spring  Rate 

lb/ft 

1.78  x  10~2 

Torsional  Spring  Rate 

lb • ft/deg 

3.16  x  10“4 

Stiffness,  El,  GJ 

lb»in2 

4.21  x  10”5 

♦Inch  units  may  be  used 

as  long  as  units 

are  consistant 

For  the  tests  conducted  under  this  program  several  modifications 
were  made  to  the  model:  The  original  remote  control  collective 
pitch  system  was  replaced  with  a  system  capable  of  feathering 
rates  of  up  to  approximately  150  degrees  per  second.  The  rate 
could  be  varied  (remotely)  from  2  degrees  per  second  to  the 
maximum.  This  allowed  the  proprotor  feather  or  unfeather  time 
to  be  varied  from  30  seconds  to  0.4  seconds.  The  new  collective 
pitch  control  system  also  had  twice  the  range  of  the  original 
system,  approximately  65  degrees  of  collective  pitch.  An  auto¬ 
matic  flapping  lockout  was  installed  in  the  model  to  permit 
investigation  of  the  effect  of  locking  out  flapping  freedom  as 
the  blades  were  feathered.  The  rpm  at  which  the  lockout 
occurred  could  be  preset.  Flapping  could  also  be  locked  out 
over  the  entire  rpm  range.  These  provisions  are  shown  in 
Figure  1-2. 
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gure  1-2  Pylon  Con 


For  the  fold-unfold  tests  the  model  was  modified  as  follows: 

The  hub  gimbal  freedom  and  the  rotational  freedom  of  the 
proprotor  were  locked  out.  A  new  set  of  grips  for  the  blades 
were  installed.  These  contained  elements  which  provided 
flexibility  about  the  feathering  axis.  One  cruciform  was 
located  inboard  of  the  fold  hinge,  the  other  outboard  to 
simulate  different  control  arrangements.  The  fold  hinge  also 
contained  a  cruciform  element  to  sir  ulate  fold  hinge  flexi¬ 
bility.  Several  cruciform  spring  rates  were  available.  In 
the  Bell-NASA  test  a  2000  in-lb/rad  cruciform  was  used  at  the 
fold  hinge,  and  1000  in-lb/rad  cruciforms  as  the  inboard 
feathering  axis  element.  Effectively  rigid  elements  were 
employed  in  the  outboard  feathering  element.  Figure  1-3  shows 
this  arrangement. 

It  should  be  noted  that  while  the  model  tested  is  based  on  a 
tilting-proprotor  design,  the  factors  which  significantly 
influence  dynamic  characteristics  (i.e.  stiffness  and  mass 
distributions  and  geometry)  approximate  Bell  Helicopter  Company's 
folding-proprotor  configuration. 

Model  Dimensions 

Figure  1-4  is  a  three-view  of  the  model  showing  the  principal 
dimensions. 

Proprotor  Parameters 

The  proprotor  blades  consist  of  an  aluminum  spar  covered  with 
fiberglass  skin.  Expanded  foam  is  used  to  maintain  the  blade 
contour.  Lead  weights  are  bonded  to  the  spar  to  achieve  the 
scaled  weight  distribution  and  mass  balance.  The  blade  pitch- 
change  axis  lies  on  the  quarter-chordline  of  the  blade. 

The  proprotor  blades  are  attached  to  a  stainless  steel  yoke 
ring  which  in  turn  is  attached  to  the  mast  with  a  gimbal  (Hooke's 
joint)  •  The  blade  centrifugal  force  is  carried  by  wound  wire 
straps  which  allow  the  blades  to  feather  on  the  yoke  spindles. 
This  semirigid  arrangement  allows  cyclic  flapping  (tip-path- 
plane  tilt)  but  restrains  coning  and  lead-lag  motion.  The  first 
inplane  natural  frequency  of  the  blade  is  above  the  rotor 
rotational  speed  to  preclude  mechanical  instability. 

The  hub  is  lightly  restrained  in  the  cyclic  flapping  sense  by 
a  hub  restraint  of  336  in-lb/radian  about  the  yoke  gimbal 
center. 

The  constant  chord  blades  are  twisted  -27.7  degrees  (hub  to 
tip)  and  the  pitch-change  axis  is  preconed  3  degrees.  The 
blades  have  a  14.6  percent  root  cutout. 
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Table  I-II  gives  the  proprotor  mass  and  stiffness  distributions. 
The  proprotor  natural  frequencies  are  shown  in  Figures  1-5 
and  1-6.  These  are  grouped  in  terms  of  collective  and  cyclic 
modes;  the  collective  modes  are  those  involving  polar  symmetric 
motion  about  the  shaft  axis,  the  cyclic  modes  are  those 
involving  polar  asymmetric  motion. 

Wing  Parameters 

The  wing  spar  is  a  hollow  aluminum  machined  spar  (bonded  top  and 
bottom  hat  sections)  which  scales  the  beamwise,  chordwise, 
and  torsional  stiffnesses  and  duplicates  the  sweep  and  elastic 
axis  location  of  the  full-scale  wing.  Segmented  balsa  fairings 
of  typical  flutter  model  construction  are  used  to  obtain  the 
aerodynamic  contours.  Lead  weights  are  distributed  along  the 
spar  to  obtain  the  correct  weight  distribution. 

The  wing  is  attached  to  the  mount  with  a  yoke  which  clamps 
the  spar  at  each  side  of  the  fuselage.  The  wing  stiffnesses 
and  mass  distribution  is  shown  in  Figures  1-7,  1-8,  and  1-9. 

Pylon 

The  pylon  is  attached  to  the  spar  with  a  tapered  conversion 
spindle  at  the  rear  of  the  spar  and  a  pitch/yaw  stop  at  the 
front  of  the  spar.  The  pylon  stiffness  is  not  scaled  and  may 
be  taken  as  effectively  rigid.  Fiberglass  fairings  (not  shown) 
are  used  to  provide  the  pylon  aerodynamic  contour.  The  pylon 
weight  and  inertia  data  are  given  in  Table  I-III. 
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TABLE  I- I I 


MODEL  BLADE 

MASS  AND  STIFFNESS 

PROPERTIES 

%R 

EIb 

lb  •  in2 

EIC** 
lb  *  in2 

wt/in*** 

lb/in 

Chordwise 
co  %  chord 

2.5 

107,800* 

149,000 

0.110 

— 

7.5 

35,000 

35,000 

0.0430 

— 

12.5 

122,000*** 

163,000 

0.0454 

22.00 

17.5 

75,000 

175,000 

0.04031 

27.00 

22.5 

39,500 

132,000 

0.03093 

27.25 

27.5 

23,000 

120,000 

0.02009 

28.25 

32.5 

13,500 

112,000 

0.01763 

28.50 

37.5 

12,300 

106,400 

0.01645 

28.50 

42.5 

10,450 

101,300 

0.01558 

28.50 

47.5 

8,700 

97,000 

0.01820 

24.50 

52.5 

7,200 

93,500 

0.01445 

28.00 

57.5 

6,000 

89,400 

0.01712 

27.75 

62.5 

5,000 

85,000 

0.01756 

26.00 

67.5 

4,050 

81,400 

0.01784 

24.75 

72.5 

3,200 

78,000 

0.02153 

23.75 

77.5 

2,400 

74,000 

0.02061 

23.25 

82.5 

1,700 

70,000 

0.01963 

22.75 

87.5 

1,200 

66,000 

0.01916 

21.50 

92.5 

850 

63,000 

0.01911 

20.50 

97.5 

600 

58,000 

0.02233 

23.05 

♦Estimated 

♦♦Computed,  except  Station  22.5  through  52.5  percent  radius 
test  data  used  in  this  area 

♦♦♦Calculated 
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TABLE  I-III 


PYLON  WEIGHT  AND  INERTIA  PROPERTIES 

Feather/Stop  Configuration 

Pylon  Weight  8 . 8  pounds 

Pylon  Center  of  Gravity  7.1  inches  ahead 

of  conversion  axis 

Pylon  Inertia  about  eg 

Center  of  Gravity  0.0646  slug-ft2 

Fold/Unfold  Configuration 

Pylon  Weight  9.25  pounds 

Pylon  Center  of  Gravity  7.5  inches  ahead 

of  conversion  axis 

0.0755  slug-ft2 


Pylon  Inertia  about 
Center  of  Gravity 


T.'ifunnn  ■  - .  ^^-v^i^asife«»fjaag^»iaLii  a  e>u  jmpi*  ■<w 


Figure  1-5.  Model  Proprotor  Collective  Modes. 


225 


WEIGHT  DISTRIBUTION  -  LB/IN 


m ms**** 


BUTTLINE  -  INCHES 


Figure  1-8.  Model  Wing  Mass  Distribution 
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CHORDWISE  MASS  LOCATIONS  FROM  1/4  CHORD  -  INCHES 


SWMWfl  A  ifU.  *1**1  V.-W**-**.. 


o  10  20  30  40  50 

WING  STATION  -  INCHES 


Figure  1-9.  Model  Wing  Chordwise  Mass 
Distribution. 


Ins  tr umen ta tion 


The  model  was  instrumented  to  provide  the  data  required  for 
the  correlation  study  and  to  monitor  critical  loads.  Table 
I-IV  lists  the  instrumentation. 

The  test  data  were  recorded  on  CEC  Model  124  direct  write 
oscillographs  and  on  magnetic  tape.  These  data  are  stored  at 
the  NASA-Langley  Aeroelasticity  Branch. 

Model  Data  Summary 

Model  Natural  Frequencies  in  Stop/Start  Configuration 

Prior  to  the  wind  tunnel  test  the  model  natural  frequencies 
and  mode  shapes  were  determined  at  the  contractor's  facility. 
The  wing-tip  spring  rates  and  pylon  inertial  properties  were 
also  measured.  Table  I-V  gives  these  data. 

The  natural  frequencies  were  also  determined  with  the  model 
mounted  in  the  wind  tunnel.  Table  I-VI  gives  the  frequencies 
and  mode  shapes.  The  variance  with  those  given  in  Table.  I-V 
is  attributed  to  the  differences  in  the  model  mounting  arrange¬ 
ment.  This  is  caused  by  a  significant  difference  in  the  wing 
root  constraint. 

The  proprotor  natural  frequencies  were  determined  during  the 
first  tests  of  the  model  during  the  Army  Composite  Aircraft 
Program.  Since  the  blades  were  not  modified  for  this  test, 
the  natural  frequencies  from  that  test  are  still  applicable. 
Table  I-VII  lists  the  zero  rpm  blade  natural  frequencies.  The 
target  frequencies  from  Figures  1-5  and  1-6  are  given  for 
reference. 

Model  Natural  Frequencies  in  Fold/Unfold  Configuration 

The  measured  natural  frequencies  and  mode  shapes  with  the 
blades  removed  and  replaced  by  equivalent  lead  weights  are 
given  in  Table  I-VIII. 

With  the  blades  installed  the  natural  frequencies  are  dependent 
upon  the  blade  fold  angle.  Figure  40,  in  Section  III,  shows 
the  variation  in  frequency  with  fold  angle  for  a  fold  hinge 
stiffness  of  2000  inch-pound/radian  and  an  inboard  feathering 
axis  spring  rate  of  1000  inch-pound/radian.  Th°se  spring 
rates  were  chosen  to  be  tested  since  they  result  in  the  blade 
natural  frequencies  crossing  wing  natural  frequencies  at  some 
point  in  the  folding  cycle.  Consequently,  they  were  expected 
to  have  the  greatest  effect  on  stability. 
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TABLE  I -IV 

MODEL  INSTRUMENTATION 


Parameter  Measured 

Location  of  Pickup 

Transducer  Description 

Wing 

Beam  Bending 

15%  and  50%  Span 

( 

Chord  Bending 

28%  Span 

{350ft,  4Arm  Strain  Gage 

Torsion 

15%  Span 

( 

Pylon 

Forward  Vertical 

Top  of  Transmission 

! 

Acceleration 

Case  | 

Aft  Vertical 

Conversion  Axis  | 

Acceleration 

(Vertical) 

iEndevco  Model  #2264 

Longitudinal 

Mast  Axis  ^ 

(Accelerometers 

Acceleration 

kRange  150g 

Forward  Lateral 

Side  of  Trans-  ; 

^Sensitivity  1.97mv/g 

Acceleration 

mission  Case  j 

Aft  Lateral 

Conversion  Axis  | 

Acceleration 

(Lateral)  ' 

Rotor 

Blade  Beam  | 

f 

>35%,  49%,  and  65% 

350ft,  4Arm  Strain  Gage 

Bending  J 

Blade  Chord 

Radius,  Master 

1  Blade 

350ft,  4Arm  Strain  Gage 

Bending 

Pitch  Link 

Master  Blade 

125ft,  4Arm  Strain  Gage 

(Cruciform) 

Blade  Flapping 

Master  Blade 

125ft,  4Arm  Strain  Gage 

Collective  Pitch 

— 

lKft  Potentiometer 

Rotor  RPM 

Magnetic  Pick-Up 

All  data  recorded  on 

direct  readout  oscillographs. 

Collective  pitch  and 

rpm  displayed  at  the 

model  control  console.) 

High  speed  movie  cameras  mounted  in  the  rotor  plane  and 

upstream  for  photographic  record. 
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TABLE  I-V 

MODEL  PROPERTIES  AND  FREQUENCIES  MEASURED 
PRIOR  TO  WIND  TUNNEL  TEST 


Property 

Wing  Tip  Beamwise  Spring  Rate1 

Wing  Tip  Chordwise  Spring  Rate1 

Wing  Tip  Torsional  Spring  Rate1 

Pylon  Weight2 

Pylon  Center  of  Gravity 

(distance  forward  of  conversion  axis)2 

Pylon  Inertia  about 
Center  of  Gravity3 


Measured  Value 
33.3  lb/in 
94.5  lb/in 
9610  in-lb/rad 
8.8  lb 
7.1  inches 

0.0646  slug-ft2) 


1  References  to  proprotor  shaft 

2  Includes  all  items  outboard  of  wing  tip,  including  proprotor 

3  Does  not  include  blade  flapping  inertia 


Natural  Mode4 

Frequency 

Mode  Shape5 

Wing  Beam 

5.48 

cps 

VZCA  *  °-0196 

rad/in 

Wing  Chord 

9.17 

cps 

VXCA  *  °-054 

rad/in 

Wing  Torsion 

18.6 

cps 

VZCA  =  -°-U7 

rad/in 

Pylon  Yaw 

40 

cps 

^p/XCA  =  “°*30 

rad/in 

4  With  blades  removed  and  replaced  by  equivalent  concentrated 
weights 


s 


Mode  shape  is  expressed  in  terms  of  pylon  pitch  per  inch  of 
vertical  deflection  at  the  conversion  axis  and  pylon  yaw 
per  inch  of  chordwise  deflection  at  the  conversion  axis 


TABLE  I-VI 


MODEL  NATURAL  FREQUENCIES  AND  MODE  SHAPES  IN  WIND  TUNNEL 


Natural  Mode 

Frequency 

Damping 

Mode  Shape 

Wing  Beam 

5.8  cps 

0.005 

ep/zCA 

=  0.0168  rad/in 

Wing  Chord 

9.0  cps 

<0.01 

VXCA 

=  0.040  rad/in 

Wing  Torsion 

17.8  cps 

<0.01 

vzca 

=  -0.130  rad/in 

Pylon  Yaw 

Second  Beam 

36.5  cps 

74.0  cps 

— 

^p/XCA 

=  -0.16  rad/in 

r 

TABLE  I-VII 

PROPROTOR  NATURAL  FREQUENCIES 

ec  =  0  RPM 

Frequency  (cps) 

Mode 

Measured 

Target 

First  Collective 

14.2 

15.2 

Second  Collective 

43.0 

46.8 

First  Cyclic 

13.5 

14.6 

Second  Cyclic 

39.2 

42.7 

Third  Cyclic 

88.6 

99.4 

TABLE  I -VI I I 

MEASURED  NATURAL 

FREQUENCIES  OF  FOLD/UNFOLD  CONFIGURATION 

Mode 

Frequency 

Mode  Shape 

Wing  Beam 

5.7  cps 

Vzca  -  °-o6:!5 

Wing  Chord 

8.8  cps 

=  0.0785 

Wing  Torsion 

18.0  cps 

0p/zCA  =  0.0358 

Pylon  Yaw 

36.5  cps 

— 

0.0358 


Model  Configurations 


For  the  rotating  tests  a  "basic"  configuration  was  defined 
as  follows: 

-  Wing  aerodynamic  panels  installed  with  weights 
simulating  15  percent  fuel 

-  Pylon  fairing  and  spinner  removed 

-  Pylon  rigidly  attached  to  wing  tip  in  pitch  and 
yaw  senses 

-  Pitch-flap  coupling  of  0.414  (63  =  -22*5  degrees) 

-  Isotropic  hub  restraint  of  338  inch-pound/radian, 
zero  degrees  angle-of-attack 

Other  configurations  tested  involved  the  following  changes  to 
the  basic  configuration: 

-  Flapping  locked  out 

-  Variation  in  hub  restraint  from  0  to  1200  inch-pound/ 
radian 

-  Pitch-flap  coupling  varied  from  63  =  0  degrees  to 
63  =  -39  degrees 

For  the  blade  folding  tests  the  model  was  tested  with  an 
effectively  rigid  pitch  change  axis  element  and  an  element 
with  a  spring  rate  of  1200  inch-pound/radian.  The  model  was 
also  tested  with  a  pylon  fairing  with  faired  blade  fold  slots, 
to  study  buffet  at  large  fold  angles. 
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APPENDIX  II 


D270A  DESCRIPTIVE  DATA 


The  D270A  is  a  folding-proprotor  V/STOL  configuration  with  a 
design  gross  weight  of  66,00  pounds.  A  three  view  is  shown  in 
Drawing  D270A-900-110 .  Three  modes  of  operation  are  possible: 
helicopter,  proprotor,  and  airplane.  Four  Lycoming  LTC4V-1 
convertible  fan/shaft  provide  power  in  rotary  wing  and  fixed 
wing  modes.  Two  fifty-foot  diameter  proprotors  provide  a 
design  disc  loading  of  16.8  psf  which  is  optimum  from  a  pro¬ 
ductivity  standpoint.  The  design  wing  loading  is  85  psf,  near 
optimum  for  a  cruise  condition  of  400  knots  at  30,000  feet. 

The  primary  dynamic  system,  consisting  of  the  proprotors,  wing, 
and  engines  and  transmissions,  has  been  designed  to  be  suitable 
for  fuselages  up  to  13-foot  diameter.  The  primary  mission  for 
the  D270A  is  a  rescue  and  a  10-foot  diameter  fuselage  has  been 
selected,  primarily  to  reduce  drag. 

Descriptive  Data  and  Dimensions 

An  inboard  profile  of  the  D270A  is  shown  in  Drawing  D270-470-110 . 
General  descriptive  data  are  summarized  in  Table  II-I. 

Performance 


A  performance  summary  for  the  D270A  is  given  in  Table  I-II. 
Structural  Description 

Proprotor  and  Pylon  -  The  rotors  of  the  D270A  are  mounted  at 
the  wing  tips  on  tiltable  pods.  In  addition  to  the  rotor, 
each  pod  contains  a  rotor  planetary  transmission,  rotating  and 
nonrotating  rotor  controls,  a  bevel  gearbox,  the  pod  tilt 
actuator,  spinner,  cowling  and  accessories.  The  spinner  and 
stationary  cowling  are  recessed  to  contain  the  trailing  edges 
of  the  rotor  blades  when  the  rotor  is  folded.  The  Inboard 
Profile,  D270A-470-110 ,  shows  the  rotor  installation  at  each 
wing  tip.  The  orientation  of  the  rotor  system  on  the  pod  is 
shown  in  Drawing  D270A-410-002 . 

The  proprotor  blades  are  all  metal,  bonded  honeycomb  structure 
with  constant  chord,  linear  twist,  and  a  nonlinear  beam  thick¬ 
ness  taper.  Drawing  D270A-410-001  shows  the  blade  arrangement. 
In  the  primary  structure,  17-7PH  stainless  steel  is  used.  The 
interior  of  the  blade  is  filled  with  an  aluminum  honeycomb 
core,  including  the  spar.  The  honeycomb  in  the  spar  prevents 
buckling,  permitting  the  use  of  thin  gage  material  in  this 
area.  An  aluminum  housing  for  the  pitch  change  bearings  is 
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TABLE  I I- I 

D270A  CHARACTERISTIC 

DATA 

Item 

Units 

Data 

Proprotor 

Blade  Airfoil  Sections 

Root  (CL  Mast  to  R/R) 

— 

64X18 

Tip 

Rotor  Tip  Rotational  Mach  No.  at 

95°F 

■  —  — 

64X08 

Helicopter  Mode 

At  100  percent  rpm 

— 

0.71 

At  90  percent  rpm 

— 

0.64 

f  Conversion  Mode 

At  95  percent  rpm 

— 

0.675 

At  85  percent  rpm 

— 

0.603 

Design  C-p/a 

— 

0.10 

Design  Disc  Loading 

At  Design  Gross  Weight  of  6600 

lb 

lb/ft2 

16.8 

Design  Blade  Twist  (Linear) 

deg 

25 

Solidity 

— 

0.1275 

Number  of  Blades  per  Proprotor 

— 

3 

Diameter 

ft 

50 

Disc  Area  per  Proprotor 

sq  ft 

1963.5 

Blade  Chord 

in 

40.0 

Blade  Area  (3  Blades) 

sq  ft 

250 

Proprotor  Flapping  Clearance 

(Helicopter) 

deg 

9.5 

Hub  Spring  Rate  (Helicopter) 

f t-lb/deg 

1500 

Folding  Hinge  Radius 

in 

26.75 

Hub  Precone  Angle 

deg 

2 

1  Fuselage-Proprotor  Tip  Clearance 

(Airplane) 

in 

32  1 

1  Tip  Speeds 

100  percent  rpm  -  Helicopter 

ft/sec 

822 

95  percent  rpm  -  Helicopter 

ft/sec 

780 

95  percent  rpm  -  Conversion 

ft/sec 

740 

85  percent  rpm  -  Conversion 

ft/sec 

698 

Drive-System  Ratios 

Ratio  of  Engine  to  Proprotor  rpm 

Ratio  of  Engine  to  Interconnect  Shaft 

— 

31.8:1 

1 _ EE2 _ 

— 

1.34:1 

J 

f 


r 


! 
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TABLE  II-I  -  Continued 

Item 

Units 

Data 

Engine 

I  Number  of  Engines 

— 

4 

Designation  (Convertible  Fan/Shaft 

Version) 

— 

LTC4V-1 

Fan  Pressure  Ratio 

— 

1.4:1 

Fan  Bypass  Ratio 

— 

8.0:1 

Military  Rating  (30  Minutes) 

Horsepower  (SLS) 

shp 

4420 

RPM  (Output  Shaft) 

rpm 

10000 

Transient  Torque 

ft-lb 

2960 

Thrust  (30000  ft,  Mach  0.7) 

lb 

1840 

1 

Airframe 

General 

Overall  Length 

ft 

81.25 

Overall  Width  (Proprotors  Turning) 

ft 

114.25 

Overall  Height 

ft 

30.66 

Distance  Between  Proprotor 

ft 

64.25 

Static  Ground  Line  Reference 

Water  Line 

in 

40 

Height  of  Conversion  Pivot 

Above  Reference  Water  Line 

ft 

13.33 

Distance  from  Conversion  Pivot  to 

Horizontal  Tail  1/4  MAC 

ft 

39.88 

Vertical  Tail  1/4  MAC 

ft 

35.66 

Distance  from  Wing  1/4  MAC  to 

Horizontal  Tail  1/4  MAC 

ft 

40.88 

Vertical  Tail  1/4  MAC 

ft 

36.66 

Fuselage 

Fuselage  Station  (0.0) 

Reference  Ahead  of  Nose 

in 

130.0 

Length  (To  Tip  of  Vertical  Fin) 

ft 

81.25 

Maximum  Breadth 

ft 

10.00 

Maximum  Depth 

ft 

10.00 

Empennage 

Vertical  Tail 

Total  Area 

ft2 

205 
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TABLE  I I- I  -  Continued 

Item 

Units 

Data 

Aspect  Ratio 

— 

1.6 

Mean  Aerodynamic  Chord  (MAC) 

in 

139.18 

1/4  MAC  Fuselage  Station 

(WL  291) 

in 

954 

Leading  Edge  Sweep 

deg 

35 

Tip  Chord  (WL  408) 

in 

102.5 

Root  Chord  (WL  192) 

in 

170.0 

Tip  Airfoil  Section 

— 

64008 

Root  Airfoil  Section 

— 

64015 

Rudder  Area 

ft2 

43.9 

Rudder  Chord/Total  Chord 

percent 

25 

Horizontal  Tail 

Total  Area 

ft2 

250 

Aspect  Ratio 

— 

4.5 

Incidence  (Trimmable  Range) 

deg 

20 

Mean  Aerodynamic  Chord  (BL  93.5) 

in 

92 

1/4  MAC  Fuselage  Station 

in 

1004.5 

Leading  Edge  Sweep 

deg 

17 

Tip  Chord  (BL  200) 

in 

75 

Root  Chord  (BL  0.0) 

in 

105 

Tip  Airfoil  Section 

— 

64008 

Root  Airfoil  Section 

— 

64012 

Elevator  Area 

sq  ft 

51.25 

Elevator  Chord/Total  Chord 

percent 

26.5 

Wing 

Geometric  (Between  Rotor  ) 

Span 

ft 

64.25 

Area 

sq  ft 

706 

Aspect  Ratio 

— 

5.85 

Root  Chord  (BL  0.0) 

in 

157 

Tip  Chord  (BL  385.0) 

in 

107 

Mean  Aerodynamic  Chord  (BL  188.6) 

in 

133.4 

1/4  MAC  Fuselage  Station 

— 

514 

Foot  Airfoil  (BL  0.0) 

— 

64X18 

Tip  Airfoil  (BL  385.0) 

- — 

64X15 

Leading  Edge  Sweep 

deg 

6 

Dihedral 

deg 

2 

Root  Incidence 

deg 

3.7 

* 


■mm 
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TABLE  II-I  -  Continued 


Item 

Units 

Data 

Flaperon 

Area/Side 

sq  ft 

23.3 

Span 

in 

103 

Chord 

in 

32.6 

Flap 

Area/Side 

sq  ft 

34.94 

Span 

in 

148.0 

Chord 

in 

34.0 

Pylon 

Conversion  Axis/Rotor  Shaft  Axis 

Intersection 

Fuselage  Station 

— 

526.0 

Waterline 

— 

200.0 

Buttline 

— 

385.0 

Conversion  Axis  Wing  MAC  Location 

percent 

35.75 

Conversion  Axis  Sweep  (Forward) 

deg 

6 

Conversion  Axis  to  Rotor  Flap  Axis 

in 

122.0 

Detail  Strength  Requirement 

Design  Gross  Weight 

lb 

66000 

Alternate  Gross  Weight 

lb 

82500 

Positive  Limit  Load  Factors  at 

Design  Gross  Weight 

Helicopter 

g 

2.5 

Conversion 

g 

.> ,  0 

Airplane 

g 

.} .  0 

Negative  Limit  Load  Factors  at 

Design  Gross  Weight 

Helicopter 

g 

-0.5 

Conversion 

g 

-1.0 

Airplane 

g 

-1.0 

Design  Sinking  Speed  at  Design 

Gross  Weight 

ft/sec 

12 

1 - 

TABLE  II-II 

D270A  PERFORMANCE  SUMMARY 

Desian  Gross  Weight 

! 

66,000  lb 

Hover  Ceilings  at  Design  Gross  Weight 

Standard  Day,  OGE 

13,900  ft 

Standard  Day,  IGE 

14,800  ft 

95°F  Day,  OGE 

7,950  ft 

95°F  Day,  IGE 

8,900  ft 

Rate  of  Climb,  Design  Gross  Weight,  SLS 

Helicopter  Mode 

4,150  ft/min 
at  70  kt 

Proprotor  Mode 

5,050  ft/min 
at  135  kt 

Airplane  Mode 

4,200  ft/min 
at  240  kt 

Speed  at  Design  Gross  Weight 

Helicopter  Mode ,  Maximum 

140  kt 

Proprotor  Mode,  Maximum 

324  kt 

Airplane  Mode ,  Maximum 

452  kt 

Airplane  Mode,  Maximum  Cruise 

415  kt  per 
3000  ft 

Airplane  Mode,  Best  Range 

340  kt  per 
3000  ft 

Alternate  Gross  Weight 

82,400  lb 

Maximum  Gross  Weight  (VTOL,  IGE,  SLS) 

93,500  lb 

Ferry  Range  (30  minute  reserve) 

2,600  nm 

bonded  into  the  blade  at  the  root  end.  The  integral  housing 
eliminates  the  weight  and  softness  of  a  bolted  joint  in  this 
area  and  provides  a  clean  aerodynamic  configuration. 

Rotor  blade  pitch  change  occurs  about  a  spindle  located  out¬ 
board  of  the  folding  hinge.  Oil  lubricated  needle  bearings 
carry  bending  loads  from  the  blade  into  the  spindle. 

The  proprotor  hub  is  designed  for  a  three-bladed ,  hub- 
restrained  ,  gimbal-mounted  semirigid  rotor.  Centrifugal  force 
is  transmitted  from  the  blade  to  the  foldable  hub  spindle 
through  a  wire  tension-torsion  strap.  The  strap  is  pinned  to 
the  blade,  and  has  a  fitting  with  an  integral  shoulder  bearing 
against  a  counter  bore  in  the  spindle  for  attachment  to  the 
hub. 

The  principal  structural  element  of  the  hub  consists  of  a 
titanium  ring  with  lugs  which  attach  the  folding  pitch  change 
spindles.  The  yoke  is  attached  to  the  mast  with  a  gimbal 
joint.  Oil  lubricated  needle  bearings  accommodate  rotor 
flapping;  as  an  alternate,  laminated  elastomeric  bearings  can 
be  used  in  this  area  to  eliminate  the  lubrication. 

In  the  helicopter  and  proprotor  modes  of  flight,  hub  restraint 
is  provided  by  a  nonrotating  elastomeric  spring  mounted  below 
the  rotor.  This  is  shown  in  Drawing  D270A-410-002 .  A  ball 
bearing  transmits  rotor  flapping  to  the  spring.  The  hub  spring 
rate  is  3840  foot-pound/degree .  Before  rotor  rpm  is  decreased, 
hydraulic  cylinders  below  the  hub  spring  lock  out  the  flapping 
freedom. 

A  blade  folding  hinge  is  provided  between  the  yoke  and  the 
blade.  Drawing  D270A-410-003  shows  the  configuration.  Folding 
is  accomplished  inboard  of  the  pitch  change  axis  to  permit 
tracking  adjustment  to  be  made  to  the  blade  without  affecting 
the  alignment  of  the  blade  when  folded,  and  so  that  aerodynamic 
and  dynamic  moments  generated  during  folding  are  isolated  from 
the  control  system.  Also,  folding  inboard  of  the  feathering 
axis  facilitates  driving  of  the  folding  mechanism.  The  hinge 
consists  of  a  multiple  clevis  attachment  below  the  feathering 
axis,  actuated  by  an  overcenter  toggle  linkage  above  the 
feathering  axis.  A  planetary  power  hinge  drives  the  toggle 
linkage.  Dual  hydraulic  motors,  attached  to  the  rotor  hub, 
power  the  system.  A  scissor  arrangement  between  the  pitch  horn 
and  the  blade  transmits  control  motion  across  the  folding  hinge. 
The  pitch  horn  is  split  for  assembly  purposes  and  is  supported 
by  split  angular  contact  bearings  around  the  yoke  spindle. 

Figure  II-l  shows  the  proprotor  stiffness  and  mass  distributions. 
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Figure  II-l.  D270A  Proprotor  Stiffness  and  Mass 
Distributions . 
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WEIGHT  DISTRIBUTION  -  LB/IN 


The  pylon  conversion  mechanism  consists  of  a  screw  jack  which 
attaches  the  transmission  case  to  the  forward  part  of  the 
wing.  By  retracting  this  screw,  the  pod  assembly  may  be  con¬ 
verted  from  the  helicopter  to  the  proprotor  mode.  Interconnect 
shafting  between  the  rotors  on  the  conversion  actuator  assures 
that  the  right  and  left  pylons  convert  simultaneously.  The 
interconnect  shaft  also  drives  the  range  shift  actuator  in 
the  fuselage  which  rephases  the  control  system  during  con¬ 
version. 

Wing  -  The  wing  is  a  two-spar,  single-cell  structure  using 
aluminum  honeycomb  sandwich  in  the  upper  and  lower  surfaces 
and  in  the  rear  spar.  The  wing  sections  and  structures  are 
indicated  in  Drawing  D270A-470-110 .  The  ailerons  and  flaps 
are  rib-stiffened  aluminum  sheet.  Concentrated  loads  from 
the  engine  installations  and  tip-mounted  transmission-rotor 
assemblies  are  carried  entirely  by  the  torque  box.  All  other 
contour  panels  and  ribs  carry  airloads  only. 

The  main  structural  box  of  the  wing  consists  of  a  constant 
section  64X17  airfoil.  The  front  spar  is  located  at  5  percent 
chord,  and  the  rear  spar  at  50  percent  chord. 

Forward  of  the  front  spar,  the  wing  contains  a  shaft  extending 
from  wing  tip  to  wing  tip  to  interconnect  the  pod  tilt 
actuators.  Access  to  the  couplings  for  maintenance  of  this 
shafting  is  through  removable  panels  on  the  lower  surface  of 
the  wing.  The  trailing  edge  of  the  wing  is  made  up  almost 
entirely  of  flap  and  aileron  assemblies  which  are  supported 
by  hinge  fittings  attached  to  the  torque  box,  rotor  drive 
interconnect  shafting  and  supports,  engine-to-interconnect 
shaft  gearboxes,  rotor  cyclic  and  collective  control  tubing, 
wing  flap  and  aileron  controls  and  actuators,  and  hydraulic, 
electrical,  and  fuel  lines. 

The  wing  strength  requirements  were  dictated  by  a  2.5g  takeoff 
in  helicopter  mode.  (This  is  equivalent  to  5.0g  in  airplane 
mode.)  This  applies  to  beam  bending  and  torsional  strength. 

The  chordwise  strength  is  in  excess  of  the  requirement.  The 
wing  stiffness  was  obtained  by  using  the  required  bending 
material  so  as  to  achieve  a  high  torsional  stiffness.  Figures 
II-2  and  II-3  show  the  wing  stiffness  and  weight  distributions. 

Empennage  -  The  empennage  is  of  conventional  construction. 

The  horizontal  tail  is  located  approximately  at  the  midspan  of 
the  vertical  tail,  and  the  entire  stabilizer  area  is  controllable. 

The  fin  is  a  two-spar  structure  with  all  the  beam  shear  and 
bending  carried  by  the  spars.  The  interspar  coverings, 
together  with  the  spar  webs,  provide  the  torsional  shear¬ 
carrying  structure.  The  fin  spars  are  continuous  from  the  tip 
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Figure  II-2.  D270  Wing  Stiffness  Distribution 


WEIGHT  DISTRIBUTION  -  LB/IN 


BUTTLINE  -  INCHES 

Figure  II-3.  D270A  Wing  Mass  Distribution. 
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to  their  attachment  to  fuselage  bulkheads.  The  entire  fin 
structure  is  attached  to  the  fuselage  by  fittings  joining  the 
spars  and  the  fuselage  bulkheads.  The  horizontal  tail  single 
cell  structural  box  is  a  typical  skin-stringer  construction 
utilizing  7075-T6  skins,  stringers  and  spar  caps.  Ribs  are 
located  at  the  elevator  hinge  points  and  at  the  horizontal  tail- 
vertical  fin  intersection. 

Both  the  horizontal  and  vertical  tail  sizes  are  determined  by 
aircraft  dynamic  stability  requirements.  The  structural 
requirements  are  expected  to  be  influenced  by  aeroelastic 
considerations  in  high  speed  mode.  However  the  aeroelastic 
aspects  have  not  been  considered  to  date.  The  empennage 
aeroelastic  characteristics  are  not  expected  to  impact  on  the 
feather/fold  sequence. 

Fuselage 

Description  -  The  fuselage  is  a  conventional  structure,  semi- 
monocoque  configuration,  72  feet  in  length,  with  a  circular 
cross  section  10  feet  in  diameter.  The  structure  design  is 
based  on  7075-T6  aluminum  alloy  materials  consisting  of  formed 
skins,  rolled  frames,  extruded  stringers,  and  machined  forgings. 
Inboard  Profile  Drawing  D270A-470-110 ,  indicates  fuselage 
layout. 

The  frames  are  of  three  basic  types .  Splice  frames  are  used 
to  accommodate  production  and  compartment  breaks.  Forged 
bulkheads  are  used  on  each  side  of  the  main  landing  gear,  and 
one  provides  a  tie  between  the  aft  wing  spars.  Intermediate 
rolled  zee-shaped  frames  are  installed  between  the  main  frames 
on  18-inch  centers.  The  extruded  zee-shaped  stringers  are 
spliced  through  the  major  bulkheads  using  forged  splice  fittings 
which  are  the  same  for  all  locations.  The  stringers  are  riveted 
to  the  skin  and  intermediate  frames  every  9  degrees  around 
the  fuselage. 

The  fuselage  strength  requirements  are  determined  by  in  flight 
gust  loads  and  landing  loads.  A  preliminary  estimate  of  the 
fuselage  stiffness  and  mass  distribution  is  given  in  Figure 
II-4 . 
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FUSELAGE  STATION  -  INCHES 

Figure  II-4.  D270A  Fuselage  Stiffness 

Distribution. 
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WEIGHT  DISTRIBUTION  -  LB/IN 


CHARACTERISTICS 


WEIGHTS 

DESIGN  GROSS  wFIGhT 

WEIGHT  EMPTY 

OIL  ANO  TPAPPED  FUEL 

CREW 

STOL  GROSS  WEIGHT 


66,000  LB 
43,505  LB 
624  LB 
1,200  LB 
•2,500  lB 


POWER  PLANT 

MANUFACTURER  and  model 
MILITARY  RATED  POWER 
MILITARY  RATED  THRUST 
TRANSMISSION  POWER  LIMIT 


LYCOMING 

LTC4V-I 

DERIVATIV  E 

4  X  4,420 

(SL  STD) 

1^  680 

HP 

4  X  1.854 

(400  KTS 

3Q000  FEET) 

7,4  6 

LB 

c  [2  X  3,980] 

15  9?0 

HP 

RCTCR 


diameter 

500 

FT 

number  cf  blades  /  rotor 

3 

oi$:  area/ rotor 

1,963.5 

SG-  OT 

BLADE  AREA /ROTOR 

250  C 

SC  FT 

DISC  LOADING 

DESIGN  GR-'SS  A  EIGHT 

1  6.8 

LO/SQ  -  FT 

BLADE  AIRFOIL 

theoretical  root  l  rcj 

64  X  16 

TIP 

64  xOb 

blade  chord 

CONSTANT 

4  :  o 

IN 

SOLIDITY 

0.1275 

BLADE  TWIST 

t  TO  TIP  LINEAR 

2  5.0 

DEGREES 

TIP  SPEED 

HELICOPTER  MCCE  (283*314  RPM) 

740  822 

FT/SEw 

CO  NVE  Rr  iCN  MOOF  (207-298  RPM 

698  -  780 

FT/St 

WING 

SPAN 

04  25 

F  T 

AREA 

706  0 

SO  -f  T 

WING  LOADING 

DESIGN  GROSS  WEIGHT 

93.  5 

l  B/SQ-F  T 

ASPECT  RATIO 

DL  <188.6 

5.85 

MEAN  AERODYNAMIC  CHORD 

133.4 

IN. 

AIRFOIL 

6  L  CO 

64  X  18 

B  L  (385  5, 

64  X  15 

FLAP  AkEA/SlDE 

AFT  OF  HINGE 

34.94 

SC.  FT 

flaperon  area  /side 

AFT  OF  HINGE 

233 

CG-FT 

EMPENNAGE 

hCRIZCNTAl  Tail  AREA 

25C.0 

SG  FT 

ASPECT  RATIO 

4  5 

AlRFC’L 

BL  ^0.0) 

NACA  64-01? 

5  L  (200.0; 

NACA  64  8 

ELEVATOR  AREA  TOTAL 

AFT  Of  HINGE 

51.25 

SQ-ft 

vertical  tail  area 

205.0 

CQ  FT 

ACP£  CT  RATIO 

1.6 

AIRFO'L 

W  L  092 . 0)  ROOT 

NACA  64-015 

W  L  (408.0)  IIP 

NACA  64-008 

RUDDER  AREA 

Af  T  Of  HINGE 

43.9 

SG  FT 

30  66  FT 


4 
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TYPICAL  WING  SECTION 
SCALE  I/IO 


RES< 

24  l 


2" 


5.56mm  MINIGUN 
TURRET 


/—RETRIEVAL  HOOKS 


L-  ELECTRONIC  EQUIP 
TYP-  UN  PRESSURIZED 


5  56mm  MINIOUN 
TURRET 
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FOLDING  BENCHS 

.-STORAGE  EQUIPMENT 
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